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Preface 


The Space Programs Summary is a six-volume, bimonthly publication that documents 
the current project activities and supporting research and advanced development 
efforts conducted or managed by JPL for the NASA space exploration programs* The 
titles of all volumes of the Space Programs Summary are: 

Vol. I. The Lunar Program (Confidential) 

VoL II. The Planetary^lnterplanetary Program (Confidential) 

Vol. Ml. The Deep Space Network (Unclassified) 

Vol. IV. Supporting Research and Advanced Development (Unclassified) 
Vol. V. Supporting Research and Advanced Development (Confidential) 
Vol. VI. Space Exploration Programs and Space Sciences (Unclassified) 

The Space Programs Summary, Vol, VI consists of an unclassified digest of appro- 
priate material from Vols. I, II, and III; an original presentation of technical supporting 
activities, including engineering development of environmental-test facilities, and qual- 
ity assurance and reliability; and a reprint of the space science instrumentation studies 
of Vols. I and II. 
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I. Introduction 


The Deep Space Network (DSN1 >^ established by the 
NASA OJRce of Tracking and Data Acquisition, is under 
the system management and technical direction of JPL. 
The DSN is responsible for two-way communications 
with unmanned spacecraft travelling from approximately 
10,000 miles from Earth to interplanetary distances. 
Tracking and data-handling equipment to support these 
missions is provided. Present facilities permit simultane- 
ous control of a newly launched spacecraft and a second 
one already in flight. In preparation for the increased 
number of U.S. activities in space, a capability is being 
developed for simultaneous control of either two newly 
launched spacecraft plus two in flight, or four spacecraft 
in flight. Advanced communications techniques are being 
implemented to make possible obtaining data from, and 
tracking spacecraft to, planets as far out in space as 
Jupiter. 

The DSN is distinct from other NASA networks such 
as the Scientific Satellite Tracking and Data Acquisition 
Network (STADAN), which tracks Earth-orbiting scien- 
tific and communication satellites, and the Manned Space 


Flight Network (MSFN), which tracks the manned space- 
craft of the Gemini and Apollo programs. 

The DSN supports, or has supported, the following 
NASA space exploration projects; (1) Ranger, Suroet/or, 
Mariner^ and Voyager Projects of JPL; (2) Lunar Orbiter 
Project of the Langley Research Center; (3) Pioneer Proj- 
ect of the Ames Research Center, and (4) Apollo Project 
of the Manned Spacecraft Center (as backup to the 
Manned Space Flight Network). The main elements of 
the network are: the Deep Space Instrumentation Fa- 
cility (DSIF), with space communications and tracking 
stations located around the world; the Ground Communi- 
cations System (CCS), which provides communications 
between all elements of the DSN; and the JPL Space 
Flight Operations Facility (SFOF), the command and 
control center. 

The DSIF tracking stations are situated such that three 
stations may be selected approximately 120 deg apart in 
longitude in order that a spacecraft in or near the ecliptic 
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plane is always within the field of view of at least one 
of the selected ground antennas. The DSIF stations are: 


0—p spot# 
Communication 
Complox (DSCC) 

Doop Space 
Station (DSS) 

DSS serial 
designation' 

Goldstono 

Pioneer 

11 


Echo 

12 


Venus 

13 


Mart 

14 

Canborro* 

Woomera 

41 


Tidbinbilla 

42 


Booroomba’ 

43 


Johannesburg 

51 

Madrid* 

Robledo 

61 


Cebreros^ 

62 


Rio Cofio* 

63 


Cope Kennedy 
(Spacecraft Monitoring) 

71 


AKension Islond^ 
(Spocecraft Guidance 
ond Command) 

72 

'At of June 1966, DSIF serial designations (e.g., DSIF-11) ore no longer used. 
^Ptonned. 

’Station not yet oudtorized. 


^Station not yet operational. 



JPL operates the U.S. stations, and will operate the 
Ascension Island Station. The overseas stations are nor- 
mally stafiFed and operated by government agencies of 
the respective countries, with the assistance of U.S. sup- 
port personnel. 

The Cape Kennedy Station supports spacecraft final 
checkout prior to launch, verifies compatibility between 
the DSN and the flight spacecraft, measures spacecraft 
frequencies during countdown, and provides telemetry 
reception from lift-off to local horizon. The other DSIF 
stations obtain angular position, velocity (doppler), and 
distance (range) data for the spacecraft, and provide 
command control to (up-link), and data reception from 
(down-link), the spacecraft. Large antennas, low noise 
phase-lock receiving systems, and high-power transmit- 
ters are utilized. The 85-ft diameter antennas have gains 
of 53 db at 2300 MHz, with a system temperature of 
55 *K, making possible the receipt of significant data 
rates at distances as far as the planet Mars. To improve 
the data rate and distance capability, a 210-ft diameter 
antenna has been built at the Goldstone Mars Station, 
and two additional antennas of this size are planned for 
installation at overseas stations. 

In their present configuration, all stations with the 
exception of Johannesburg, are full S-band stations. The 
Johannesburg receiver has the capability for L- to S-band 


conversion. The Ascension Island Station will be basically 
full S-band when it becomes operational. 

It is the policy of the DSN to continuously conduct 
research and development of new components and sys- 
tems and to engineer them into the network to maintain 
a state-of-the-art capability. Therefore, the Goldstone sta- 
tions are also used for extensive investigation of space 
tracking and telecommunications techniques, establish- 
ment of DSIF/spacecraft compatibility, and development 
of new DSIF hardware and software. New DSIF-system 
equipment is installed and tested at the Goldstone facili- 
ties before being accepted for system-wide integration 
into the DSIF. After acceptance for general use, it is 
classed as Goldstone Duplicate Standard (GSDS) equip- 
ment, thus standardizing the design and operation of 
identical items throughout the system. 

The GCS consists of voice, teletype, and high-speed 
data circuits provided by the NASA World-Wide Com- 
munications Network between each overseas station, the 
Cape Kennedy Station, and the SFOF. Voice, teletype, 
high-speed data, and video circuits between the SFOF 
and the Goldstone stations are provided by a DSN micro- 
wave link. The NASA Communications Network is a 
global network consisting of more than 100,000 route mi 
and 450,000 circuit mi, interconnecting 89 stations of 
which 34 are overseas in 18 foreign countries. It is en- 
tirely operationally oriented and comprises those cir- 
cuits, terminals, and switching equipments interconnect- 
ing tracking and data acquisition stations with, for ex- 
ample, mission control, project control, and computing 
centers. Circuits used exclusively for administrative pur- 
poses are not included. 

During the support of a spacecraft, the entire DSN 
operation is controlled by the SFOF. All spacecraft com- 
mand, data processing, and data analysis can be accom- 
plished within this facility. The SFOF, located in a 
three-story building at JPL, utilizes operations control 
consoles, status and operations displays, computers, and 
data-processing equipment for the analysis of spacecraft 
performance and space science experiments, and com- 
munications facilities to control space flight operations. 
This control is accomplished by generating trajectories 
and orbits, command and control data from tracking 
and telemetry data received from the DSIF in near-real 
time. The SFOF also reduces the telemetry, tracking,! 
command, and station performance data recorded by the/ 
DSIF into engineering and scientific information for 
analysis and use by scientific experimenters and space- 
craft engineers. 
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II. Tracking and Navigational Accuracy Analysis 

N67 


A. DSIFTwo-Way Doppler Inherent 
Accuracy Limitations: III. 
Charged Particles 

D. W. Trask and L Efron 

1. Introduction 

This is the third in a series of articles (SPS 37-38, -39, 
VoL III, pp. 8~12, 7—15, respectively) on the inherent 
accuracy limitations to the DSIF data type commonly 
referred to as two-way doppler. In practice the doppler 
is continuously counted and sampled at fixed intervals 
such that the observable used by the orbit determination 
program is the range diflFerence between the observer 
and spacecraft, which is accumulated over a fixed count 
interval. Although the accuracy of these data types 
has consistently exceeded the quality necessary for the 
successful fulfillment of the current Earth-based radio 
guidance role to which it has been committed, an effort 
is continually in progress to improve the quality of these 
data so that needs of future missions may be met. This 
article discusses the effect of the charged particles which 
exist in the ionosphere and the space plasma between 
the spacecraft and the tracking station. 


A set of standard data-taking conditions has been 
selected for this series of articles to illustrate the inherent 
accuracy limitations on two-way doppler (f 2 ). These 
errors^ are expressed in terms of their contribution to 
the range difference derived from counting the doppler 
tone over 10^ sec. This relatively long count time 
(Tc ^ 3 hr) is chosen to show the advantage of contin- 
uously accumulating doppler in the presence of bounded 
errors such as spacecraft limit cycle motion, doppler 
counter quantization, and system phase jitter. These 
errors would look deceptively large (percentagewise) for 
small Tc. In addition, Tc is also made large to reduce the 
correlations of successive data points. Examples of the 
occurrence of relatively long correlation widths include 
the deviation between the real universe and the fitters 
model of atmospheric corrections and of the relationships 
among station time tags, universal time, and ephemeris 
time. In addition, it is assumed that the count interval 
is centered about the zenith, i.e., maximum elevation 
angle. This has the effect of minimizing the effects of 
atmospheric refraction. 

Both lunar and planetary missions are considered. For 
the lunar case the spacecraft is assumed to be at a lunar 
distance which results in a radio signal round-trip transit 


*A11 statistics are quoted as l-o- values unless otherwise specified. 
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time (t) of 3 sec while the spacecraft distance is 1.5 AU 
(t = 1,500 sec) for the planetary case. This influences, 
for example, the data sensitivity to transmitter oscillator 
stability and, in the case of the space plasma, the number 
of charged particles the radio signal encounters. 


2 . Discussion 

a. Electromagnetic wave propagation. When an elec- 
tromagnetic wave propagates through a region contain- 
ing charged particles there is an interaction between the 
ionized particles and the electric and magnetic fields 
associated with the wave. Any external magnetic field 
will complicate the interaction. The motion of an ion of 
mass m and charge e* is described by 

mr = e*E H rXH — mvf 

c 


where 


V is the effective collision frequency 
H is the magnetic field 
E is the electric field vector 


The conduction current, due to the presence of freely 
moving charged particles of density N, is given by 

I = Ne*^ = - 

at rtuo 

The conductivity of the medium <j, is the ratio 1/E and is 

_ 

^ — 

rruii 


Since a is a purely imaginary quantity, it implies that 
the current and field waves are 90 deg out of phase with 
one another, and no power is consumed. Although, to 
these approximations, the wave is not damped, its veloc- 
ity of propagation will be altered by the presence of 
the ions. 


For a wave described by 

E — a cos {iot — kx) 

the velocity of propagation, in a dispersive medium, is 
defined by the spatial variation of its phase. Surfaces of 
constant phase propagate at a velocity v, given by 


c is the speed of light 
r is the displacement of the ion 

Magnetometer measurements aboard spacecraft have 
to date, shown that the magnetic fields encountered in 
the interplanetary medium are negligible with respect to 
their effects on radio tracking signals. The collision fre- 
quency may also be assumed to be near zero. Neglecting 
terms of order (1/c) we have simply (Ref. 1). 

mr = e*E 

For a plane E M wave of amplitude Eq and angular fre- 
quency <0 we will have E = Eq and 

mr = e*Eo 

the solution of which is 



Since B merely specifies an origin and A implies a con- 
stant component of velocity, they are usually both taken 
to be zero. This result is equally valid for electrons or 
heavy ionic particles. The mass of a proton or ion will be 
eighteen hundred times that of an electron. Thus, their 
amplitudes of motion will be very much smaller than 
those of electrons. 


0 > 

We define the index of refraction n of a medium by the 
relation 

c 

n = — 

V 

Taking the dielectric constant of free space to be unity, 
it can be shown that 



At this point we note that n < 1 for both electrons or 
positive ions, and the phase velocity through the ionized 
region is greater than in free space. For any given elec^ 
tron density, it is also noted that the phase velocity 
decreases with the square of the signal frequency. 

Since monochromatic waves exist in textbooks only, 
let us look at a signal composed of two waves of equal 
amplitude, but slightly different frequencies and propa- 
gation constants. If the two waves are given by 

El = cos [<ot — kx] 

E 2 = cos [(<D + h<o)t - (k + 8k)x] 
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where 8<u and Sk are small compared with <a and 8, the sum E = Ei -!- is given by 

E = 2cos|^^t - ^xjcosj^^<o + ^y~(k + 


The result is a wave exhibiting a beat frequency. The 
carrier is essentially the original frequency Ei, but the 
envelope is characterized by a wave for which planes of 
constant amplitude have a velocity of propagation u, 
given by 

8(0 

The quantity u is referred to as the group velocity. 
Further analysis gives 

uv = 

and we see that while the phase velocity may exceed c 
in an ionized medium, the group velocity will be less 
than c. 

For N expressed in electrons per cubic meter and the 
signal frequency f in Hz we have 



The effect of a nonunity index of refraction on the dop- 
pier and ranging methods utilized in the radio tracking 
of a spacecraft is that the observed doppler will be 
advanced (the range difference accumulated by counting 
the doppler tone is decreased) by 

: 40.3 d 

and the ranging signal, which propagates at the group 
velocity, will be retarded (the extended travel time indi- 
cating a fictitious increase in the measured range) by an 
amount proportional to 



The integration is along the entire signal ray path. 

The integral may be broken down into two parts. One 
through the Earth s ionosphere and the other through the 
interplanetary medium. Additional planetary ionospheres 
encountered along the ray path would also bear con- 
sideration. 

b. Charged particle density and flux in the interplan- 
etary medium. Plasma probes aboard Mariners II and /V 
and aboard Pioneer VI have measured or are measuring 


local values of charged particle density and flux in the 
interplanetary medium. The information gathered to date 
indicates that nearly all the charged particles detected 
are in the so-called solar wind which consists of a flow 
of electrons, protons and alpha particles moving radially 
outward from the Sun. The plasma flow seems to have 
three source mechanisms. The first is a perpetual general 
outflow from the solar corona while the second is the 
emission of intermittant jets of varying intensity. The 
duration of an individual jet may be several solar rota- 
tions. These latter give the solar wind a pinwheel 
streamer configuration. The final mechanism is anomaly 
effects such as plasma emissions due to solar flares. When 
viewed from along the Suns axis of rotation, we see that 
the streamer regions of higher density take on a garden 
sprinkler appearance. 

Under the assumption of local charge neutrality the 
electron density Ne will be given by 

K, = N,(l + 2^) 

where Na/Np is the ratio of alpha particles to protons. 
From Mariner II (Ref. 2) results this ratio appears to be 
about 0.05. For a spherical outflow, conservation of mass 
requires Nvr^ — constant. This requires particle density 
flux to follow a r~2 relation with respect to distance from 
the Sun (Fig. 1) (obtained from Ref. 3). 

Mariner II measured proton densities in the interplane- 
tary medium between the orbits of the Earth and Venus. 
Proton velocities varying from 300 to 800 km/sec (3-hr 
averages) have been recorded (Fig. 2; obtained from Ref. 
3). At these velocities the gravitational deceleration be- 
tween the orbits of Venus and the Earth would be negli- 
gible, and this is indicated to be so in that Mariner II 
detected no dependence of particle velocity on helio- 
centric distance. At a velocity of 500 km/sec a particle 
leaving the Sun would reach the vicinity of the Earth’s 
orbit in 3^ days. During this time the Sun would have 
rotated approximately 45 deg. Fig. 3 illustrates how par- 
ticles flowing radially outward near the Sun’s equatorial 
plane would generate spiral streamers in the solar wind. 
The Sun’s equator is inclined by 7 deg 15 min to the 
ecliptic, and hence all the interplanetary probes to date 
would have experienced the solar wind as described. 
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Fig. 1 . Logarithmic plot of proton density and flux vs 
distance from the Sun. The measurements have 
been averaged over a solar rotation. 

Mariner // (1962) (Ref. 3) 


The flux density in the solar wind is probably related 
to the solar sunspot cycle. Pioneer VI, which was launched 
into a heliocentric orbit between the Earth and Venus 
late in 1965, after the last minimum in the sunspot cycle 
(Fig. 4), has aboard two plasma probes from MIT and 
Ames. Recent telephone conversations with members of 
the experiment teams indicate it has detected the pres- 
ence of only four or five dense streamers, while Mariner 
II, launched between the previous maximum and this 
minimum, detected the presence of many more and 
broader regions of high density. However, both spacecraft 
seem to be encountering very nearly the same average 
values of local electron density over a solar revolution. 
Local electron densities deduced from the proton density 
measurements, as shown in Fig. 2, range from 1--30 
electrons/cm^. An average near 10 electrons/cm^ is 
indicated. 

The Stanford dual frequency experiment aboard 
Pioneer VI consists of the transmission of two coherent 
signals at widely separated frequencies from Stanford to 
the probe (Ref. 4). Both signals are simultaneously mod- 
ulated by a low frequency signal. The inverse square 
frequency dependence of phase and group velocities on 
electron density makes possible the detection of relative 
carrier phase and group envelope phase differences. From 
these measurements the total integrated electron density 
and its time rate of change along the entire signal ray 
path may be deduced. Pioneer V7 is in an orbit having 
an aphelion and perihelion of 1 and 0.8 AU, respectively. 
In recent telephone conversations, Dr. Bruce Lusignan 


of Stanford indicated that preliminary results from the 
Stanford dual frequency experiment show average values 
over the entire interplanetary portion of the Stanford- 
Pioneer VI signal ray path that vary from 5-25 
electrons/cm® with an average value over a solar revolu- 
tion near 10 electrons/cm^ This wide range may be 
accounted for by the presence of only a few dense 
streamers coupled with a geometry such that the ray 
path may at times lie entirely within or outside a streamer. 

Subsequent Pioneer missions will explore the inter- 
planetary medium beyond the orbit of the Earth out to 
a heliocentric distance of 1.2 AU. Referring to Fig. 1 we 
see that extrapolation of available data indicates that the 
effects on radio tracking signals due to the presence of 
the solar wind would be greater for a typical Venus mis- 
sion than a typical Mars mission. 

For a typical Mars mission near the minimum of the 
solar cycle we may expect an average electron density of 
5 electrons/cm^ along an S-band ray path through the 
interplanetary medium from Earth to Mars. Over 1.5 AU, 
this would produce a range error of approximately 11 m at 
S-band. The time rate of change of the integrated electron 
density would depend greatly on the number of streamers 
in the solar wind and the spacecraft-Sun-Earth geometry. 
A worst condition would occur when the geometry was 
such that the ray path went from almost entirely outside 
to almost entirely inside (or vice versa) a streamer. As- 
suming a linear variation over 3 days between extremes, 
2 and 15 electrons/cm^ could result in a phase advance 
of approximately 0.011 cm/sec (===^0.0085 Hz) at 1.5 AU 
for an S-band signal. This would represent an error of 
1.1 m in the change in range over 10^ sec. 

References 2, 3, 5-8 are a source for further readings 
on solar wind research. 

The Earth's ionosphere is also dependent on the state 
of solar activity, but the major source of ions seems to be 
primarily dependent on the ionizing radiation of the Sun 
in the upper atmosphere (Ref. 9 and 10). Particle flow 
from the Sun into the planetary atmosphere is negligible, 
if not nonexistent. 

c. Ionospheric effect on range measurements. The 
ionospheric effect on range measurements is illustrated 
in Fig. 5 for the range nighttime (low sunspot index) to 
daytime (high sunspot index). This figure is based on an 
internal report by A. Liu and D. Cain which assumes a 
Chapman atmosphere and studies the variation of 8/p 
with maximum electron density, height of maximum elec- 
tron density, scale height, and temperature. The range 
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Fig. 2. Calculated proton velocity and density (3-hr averages). Mariner II (1962) (Ref. 3) 
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Fig. 3. The pinwheel solar wind pattern illustrating two 
streamers formed by charged particles flowing 
everywhere radially outward from the Sun 



1962 1964 }966 1966 1970 1972 


Fig. 4. Solar activity cycle 

correction 8/p varies linearly with the maximum electron 
density. This parameter causes the greatest variation in 
8jP and, in fact, diflFers by an order of magnitude for the 
conditions presented in Fig. 5. However, this eflFect is 
largely predictable from the sunspot index disseminated 
by the Bureau of Standards, Central Radio Propagation 
Laboratory, Boulder, Colorado, (Ref. 11 and 12) and a 
knowledge of the tracking geometry. In fact, the chief 
parameters determining the state of the ionosphere above 
any point on the Earth s surface appear to be the state 
of solar activity, geomagnetic latitude, time of day and 
time of year. The latter two are sometimes accounted for 



Fig. 5. Effact of ionosphere on rouge measurements 
DSIF $-bond system (f = 2.3 X 10» Hz) 


by considering the solar zenith angle. Present ability to 
model the ionosphere for the purpose of predicting elec- 
tron density profiles is very limited. Although it may be 
possible to predict the slope and level of the 8/p vs ele- 
vation angle curve to about 10% on a day-to-day basis, 
at a particular location based on accumulated data,^ the 
accuracy of long-range or world-wide predictions is very 
poor. Tlie correlation between stations several hundreds 
of kilometers distant is often poor. 

The correction to counted doppler accumulated over 
Tc sec is essentially the difference between the two range 
corrections Tc sec apart. In particular, the formula for 
this correction is 

~ {[^cpp{t) + ^cpp(t r)] 

[8ci»p(^ Tc) “b 8cppif T c) 

where 

€cp^^ — effect of charged particles on the counted 
doppler range difference 

8cpP = effect of charged particles on range 

i = time that accumulation of Ap is completed at 
the receiving tracking station (in practice t is 
the data time label) 

Tc — count (accumulation) time for Ap 
T = round trip signal transit time 

£ I j == expected value operator 

The above formula applies to the effect on accumu- 
lated doppler of charged particles in the ionosphere 

2This represents extrapolation from a perfect knowledge. Otherwise 
it is probably no better than 20% . 


8 





JPL SPACE PROGRAMS SUMMARY NO. 37-41, VOL. Ill 


and/or solar plasma. For the nominal conditions con- 
sidered in this article, t = 1500 sec for 1.5 AU and 
Tc = 10* sec, and it is assumed that the elevation angle 
at the end points of the doppler accumulation interval 
is greater than 30 deg. A typical time history of elevation 
angle during the same lunar mission is shown in Fig. 6 
for a deep space station (DSS) at Goldstone, California 
(latitude = +35.4 deg) and at Woomera, Australia 
(latitude = —31.4 deg). The maximum elevation angle 
(ym<w) is a function of the DSS latitude (</>) and the local 
declination (8^,) of the probe. In particular 

ymar ~ 90 deg ~ I 8l “ I 

The majority of the DSS are located at either + 35 
deg or — 35 deg. Therefore, for a distant probe 
y > 30 deg for some of the DSS at the end points of the 
Tc = 10* sec doppler accumulation period. To minimize 
the effects of variations in elevation angle, it is assumed 
that the count interval is centered about the maximum 
elevation angle. This means that the range corrections 


which are being differenced will be at nearly identical 
elevation angles, and, therefore, errors in the predicted 
slope of the correction vs elevation angle will be mini- 
mized. It is assumed that with a reasonable model for the 
ionosphere will be good to 10% (although Sjp is 
good only to 20%). The principal limitation is due to 
the fluctuations of n with time. 

In some instances the spacecraft will carry equipment 
which, along with ground based equipment, is sufiBdent 
to calibrate the total (ionosphere and space plasma) 
charged particle effect. Two approaches are available. 
The first is the dual frequency experiment of Stanford 
University (Ref. 4) which is aboard the Pioneer VI space- 
craft and has been approved for the Mariner 67 mission. 

The second approach also makes use of the fact that 
the charged particles retard the group velocity and ad- 
vance the phase velocity. The DSIF two-way doppler 
experiences the phase velocity advance while the pulse 



-7 ~6 -5 -4 -3 -2 -I 0 I 2 3 4 5 6 7 

TIME FROM ZENITH, hr 

Fig. 6. Typical elevation angle time history for a lunar mission 
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code of the DSIF ranging system suflFers the group 
velocity retardation. These two eflFects are nearly equal 
in magnitude and opposite in sign. Therefore, when the 
spacecraft carries a DSIF doppler/ranging transponder, 
the influence of the charged particles can be measured 
during the period when both range and fz are obtained. 
Missions (per present plans) which will be able to employ 
the phase-group velocity calibration technique include the 
Lunar Orbiter, Apollo, Mariner 67, Mariner 69, and 
the Voyager missions. 

This method of calibration is limited by those errors 
(to the extent they cannot be measured or predicted) 
which do not have an equal eflEect on the range difference 
obtained by integrating the two-way doppler and the 
range increment obtained by differencing the range 
measurement at the start from the one at the end of the 
doppler accumulation interval. 

The size of these combined errors is shown in Table 1 
for the Lunar Orbiter, Mariner 67 Venus, and post 1970 
missions. Time lags, tracking errors, and the resolution 
of the present ground handling system are the main 
sources of error. The measurement of charged particle 
effects was envisioned relatively recently. Specifications 
resulting from this requirement in several cases greatly 
exceed those imposed during the design of the spacecraft 
transponder and ground data handling system. Although 
this calibration technique will be exercised dxuing the 
Lunar Orbiter mission the errors of the system will be 
large compared to the effect being measured. It should 
be noted that although bounds on the error sources con- 
sidered in constructing Table 1 are reasonably well 
known, their actual variation over 10* sec has not been 
rigorously tested. The Lunar Orbiter data will aid in 
analyzing some of these characteristics. 


Table 1. Limitations to range measurements: 
field performance 


Mittlon 

Uncoitointy In pCts) — p(tv} 
(ti ~ ti = 10* soc), m 

Spocacrafl 

Ground tytiom 

lunar Orbifor 

1 - 6 

1 -7 

Morinar 67 



Standard ground tyttem 

0.5-5 

1 -7 

Rasearch ground systom 

0.5-5 

0.1 - 1 

Pott 1970 



(Noxt genorafion 

0.05 - 0.5 

0.01 - 0.1 

spacocraft transponder) 




The Mariner 67 mission will be the first chance to 
measure a sizeable solar plasma effect. In addition, the 
calibration of the spacecraft transponder will be superior 
to that performed for the Lunar Orbiter and an improved 
ground system will be in operation at the 210-ft dish at 
Goldstone on an experimental basis. This system will 
provide data of sufficient quality to make meaningful 
comparisons with the measurements obtained by the 
Stanford University dual frequency experiment which is 
part of the Mariner 67 scientific payload. 

The group/phase velocity requirements will be incor- 
porated in the specifications of the next generation space- 
craft transponders and ground data handling systems. It 
is anticipated that the charged particle effects will be 
measurable from 0.05 to 0.5 m during post 1970 missions. 

Although such plans are being made for the future, 
charged particle effects are currently ignored at JPL 
when processing the two-way doppler data. The past 
plus the predicted history of the data processing is de- 
lineated in Table 2 for the case of the ionosphere and 
space plasma. This information is shown in a different 
form in Fig. 7 starting with the Ranger III mission and 
projecting into the 1970’s. Event markers for the Ranger 
III-IX and Surveyor 1 (lunar missions) as well as for the 
launch and target encounter of Mariner II (Venus), 
Mariner IV (Mars), and Mariner 67 (Venus) are super- 
imposed on this figure. Both planetary (spacecraft at 1.5 


Table 2. Effect of charged particles on two-way 
doppler accumulated over 10* sec" 


Past and prodlctod history 
of data procosting 

Source of 
chorged 


Timo poriod 

Conditions 

particles 


1962-1964 

Uncalibroted’L band 

ionosphere 
Space plasma 

1.0 to 6.0 

3.0 to 9.0 

1964-1967 

Uncalibratod-S band 

Ionosphere 
Space plasma 

0.17 to 1.5 
0.5 to 1.5 

1967 

Modolod ionosphero 
uncolibrated 

Ionosphere 

0.1 to 1.0 

1967 

Phato-group velocity | 
calibration 

Totol 

0.5 to 1.8 

1970 

Modeled ionosphere 

Ionosphere 

0.05 to 0.5 

1970 

Phase-group velocity 
calibration 

Total 

0.05 to 0.5 


•p — 1.5 AU [Mar* MiiiJonJ and »h« \0* fc oc<umti lotion porlod it contorod 
about fho maximum olovation anglo. 
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AU s) and lunar conditions are illustrated. The planetary 
missions include the ionospheric and space plasma effect 
while the lunar conditions illustrate only the effect of the 
ionosphere. (The space plasma effect is negligible.) Under 
favorable conditions a modeled but uncalibrated iono- 
spheric correction will be accurate to 0.1 m. The calibra- 
tion of the total charged particle effect for planetary 
distances to 0.1 m or less will not be available until the 
post Mariner 67 missions. These techniques, such as the 
phase-group velocity calibration, do not separate the 
contributions due to the ionosphere and space plasma 
but measure the total charged particle effect. 

3. Summary 

By 1970 the effects of charged particles will be mea- 
surable to less than 0,1 m by the group/phase velocity 
technique. That is, the charged particle effect on the 
phase velocity of DSIF doppler is nearly equal and op- 
posite to its group velocity effect on the DSIF range 
code. Therefore, a comparison of the range difference 


obtained by accumulating doppler with the range differ- 
ence obtained from direct range measurements at the 
ends of the doppler accumulation period will provide a 
measurement of the influence of the charged particles. 
This approach will be tried during the Lunar Orbiter 
mission. But the first chance to obtain meaningful space 
plasma measurements will come during Mariner 67. The 
Mariner 67 mission will also carry a Stanford University 
dual frequency experiment for comparing the two mea- 
suring techniques. 
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B. AStudy of Planetary-Approach 
Orbit-Determination Policies Using 
Mariner IV Doppler Data 

G. W. Null 

I. Introduction 

The Mariner TV spacecraft was launched November 28, 
1964 and made its closest approach to Mars July 15, 1965 
at 01:00:57 GMT. At the time of Mars encounter (E) 
the distance from Mars was 13201 km. 

During the mission, highly accurate two-way doppler 
tracking data were taken by the DSN. During the period 
of E~5 days to E, four DSN tracking stations provided 
continuous coverage: Goldstone, California; Canberra, 
Australia; Woomera, Australia, and Johannesburg, South 
Africa. Goldstone and Canberra had regular S-band sta- 
tions while Woomera and Johannesburg operated with 
L- to S-band conversion kits, which gave a factor of 3 
poorer resolution. The Goldstone station was equipped 
with a "times 8” multiplier which multiplied the doppler 
frequency by 8 before counting, thus increasing the 
resolution. The limiting error source for all stations was 
instability in the rubidium standard. (See SPS 37-39, Vol. 
Ill, pp. 7-15 for a detailed description of the Mariner TV 
data noise and its causes.) This error source contributed 
a l-<r noise of 0.0005 m/sec for 60-sec sampled data or 
0.00025 m/sec for 600-sec sampled data. 

The accuracy of the tracking data for 60-sec data (ex- 
pressed in range rate) at Mars distance was: 

For Goldstone — 0.0005 m/sec 

For Canberra a; = 0.001 m/sec 

For Woomera and Johannesburg af = 0.003 m/sec 

When data are processed by counting over 600 sec, 
the error for all stations is <rf < 0.0005 m/sec. 

For the time span of interest, 600-sec data were taken 
from £—5 days to E~4 hr (H) and 60-sec Goldstone 
data from E— 4hrtoE+4hr. Thus, each data point had 
an accuracy of 0.0005 m/sec la. When assigning a priori 
standard deviations to the tracking, as used in orbit 
determination, the assumption was made (which will be 
substantiated by the results to follow) that low frequency, 
long correlation width noise external to the data was the 
predominant orbit error source. Therefore, for orbits, the 
600-sec data are assigned a standard deviation which is 
a factor of (lO)'^ smaller than the 60-sec data. The 
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weights (standard deviations) for the results to be dis- 
cussed are all quoted for a 60-sec sample rate. 

This article will present the orbit determination target 
error results obtained by analyzing the tracking data from 
£—5 days to the following times: 

(E-72H, E-48H, E-24H, E-17H, 

E-5H, E-IH, E-0.3H, E+0.25H) 

A weighted least squares single precision orbit deter- 
mination program was used for the analysis. The quan- 
tities of interest are the predictions of the following 
closest approach quantities 

B = a vector from the center of Mars normal 
to the incoming asymptote 

B • Te, B • Rc 

where Tc is a unit vector in the plane of the ecliptic 
normal to the incoming asymptote S. Rc makes up a 
right-handed RST coordinate system. For Mariner IV, Rc 
is very nearly normal to the elliptic plane. Finally, the 
errors in time of closest approach (TCA) are investigated. 
The evaluation of the target point prediction error, as a 
function of time, is accomplished by assuming that the 
"true” target point may be obtained by fitting data from 
£—5 days to £+5 days. A full discussion of this orbit 
and its accuracy is given in SPS 37-39, Vol. Ill, pp. 
24-30. Summarizing these results, it was shown that the 
target point was accurate to at least 20 km. Since most 
of the errors to be discussed here are much greater than 
20 km and since we are mainly interested in the devia- 
tions from the final predicted target point, it may be 
considered that the "true” target point is exact. The astro- 
nomical unit (At/) was determined to within ±500 km 
of la and the gravitational constant of Mars (fi) to within 
±8 km^/sec^. 

The reason for fitting only the near-encounter tracking 
data for encounter operations, as was done for Mariner 11 
and Mariner TV, has been the difficulty in obtaining good 
solutions over long tracking arcs. One particularly trou- 
blesome error source is small unpredictable disturbance 
accelerations from the attitude control system (^2X10“^ 
cm/sec^ for Mariner IV and 25 to 50 X 10*^ cm/sec^ for 
Mariner II). Other important error sources are the plane- 
tary ephemerides and the A17. Although the use of longer 
tracking arcs in the presence of a determinate model 
should always give better and more reliable results, the 
actual nondeterminate errors negate this conclusion for 
current probes. Certain future missions including Voyager 
have such requirements for high accuracy as to necessi- 
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tate fitting over long arcs. This requires that the ephemeris 
errors and small force errors be held to a minimum. The 
work to control these errors is being planned such that 
by 1970 an accuracy of 100 km at £ —30 days may be 
achieved. 

Although it is hoped that such error sources will be 
reduced to a negligible amount by 1970, this is by no 
means certain. Therefore, it is important to examine 
carefully the accuracy attainable using only short arc 
encounter data, both for current spacecraft and for Voy- 
ager era spacecraft. The purposes of this study are: 

(1) To compare, using the current single precision orbit 
determination program (SPODP), the actual real- 
time Mariner IV orbit determination errors with 
the errors which would have existed if alternative 
orbit determination policies had been adopted. 

(2) To obtain a better indication of the accuracy ob- 
tainable using current software and present track- 
ing data accuracies. This information is needed for 
planning for future missions in the following ways: 

(a) Demonstrating typical error patterns for the 
approach phase. 

(b) Giving a better background for “committing” 
orbit accuracy to the flight projects. Normally 
the committed accuracy will be at least as good 
as the current “state of the art” accuracy. 

(c) Showing which a priori statistical assumptions 
for data weighting and initial conditions are 
consistent with the actual errors. These as- 
sumptions may then be applied to study a wide 
range of trajectories for future missions. 

(d) Aiding the selection of an optimum orbit deter- 
mination policy for future missions. 

(3) To answer certain questions of general interest 
such as: 

(a) What accuracy in determining AU and /x is 
possible as a function of time before encounter? 

(b) At what time does the smallest B plane uncer- 
tainty (SMIA) shift from being along B • T to 
being along the B vector, as predicted by orbit 
theory? 

2 . A Priori Stathfics Used in Computer Runs 

a. Real time runs (plot symbol X), 

(1) cr; = 0.03 m/sec for 60-sec sample 

(2) '^position ~ 500 km spherical distribution at E — 5 

days 


(3) ^veiocitv = 0.1 m/sec spherical distribution at E— 5 

days 

(4) ^Av = 2000 km 

(5) = 860 kmVsec^ = 2.0% 

The errors in the a priori positions were ^ 1000 km 
and the errors in the a priori velocities were ^0.1 m/sec. 
The error in the a priori AU was ^ 1000 km and the error 
in fx was ^ 100 km^/sec^. 

b. Postflight runs, 

(1) Common characteristics 

(a) = 0.003 m/sec for 60-sec sampled 

data 

(b) (^position 00 st E 5 days 

(c) o^veiocity 00 3it E 3 days 

(2) Runs solving for position and velocity only 

(a) Use true AU, fx (plot symbol 0) 

(b) Use true AU but 

^fx = 150 km^/sec^ (no plot since results were 
oflF scale ^ 3000 km) 

(c) Use true fx but 

aAU = + 1500 km (plot symbol A) 

(3) Solve for position, velocity, /x, and AU with no 

a priori information (plot symbol ■) 

3 . Explanation of Plots 

The various types of runs shown in the plots are dis- 
cussed in Section 2. Also displayed on most plots are 
the 1-(T limits which come from a no a priori solution for 
position and velocity with ak = 0.003 m/sec. The excep- 
tion was Fig. 8 for | B | since SPODP does not compute 
(Tjs. For most parameters the error is shown in two plots; 
one for the time span E— 5 days to E— 5 hrs and the 
other for the time span E— 5 hrs to E+0.25 days. 

Among the cases not shown was a set of runs defined 
by Section 2 (b) in which A/x = 150 kmVsec^ The target 
error for this case was greater than 3000 km for the entire 
period E— 5 days to E+0.25 hrs. Other cases not shown 
include a set using simulated data with no noise, solving 
for position and velocity with AAU = 500 km and A 
Earth-Mars ephemeris of 100 km spherical error. No 
significant differences were seen between this case and 
those displayed here. 
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80 70 60 50 40 30 20 10 0 

HOURS BEFORE ENCOUNTER 

Fig. 8. Vector magnitude error vs hours from Niariner IV 
encounter E — 72 H fo E — 5H 


The figures and their contents are listed as follows: 


Fig. 

Quantities 

Time span, hrs 

8 

AB 

E-72to£-5 

9 

ab -r 

E-72toE-5 

10 

ab -r 

E— 5toE 

11 

V 

E-72toE 

12 

aAU 

E-72toE 

13 

Orientation 

E-72toE 


angle of SMI A 
from B 


14 

AB • R 

E-72toE-5 

15 

AB • R 

E- 5toE 

16 

AB 

E— 5toE 

17 

aTCA 

E-72toE-5 

18 

MCA 

E — 5 to E 

19 

SMAA and SMZA 

E-72toE-5 

20 

SMAA and SMIA 

E- 5toE 


SMAA and SMI A are the semimajor and semiminor axes 
of the 40% dispersion ellipse of the predicted target 
point in the B plane. 

4. Significant Features Seen in the Plots 

Comparing the 1-a limits with the actual errors, it may 
be seen that an a priori data weight of a; — 0.005 m/sec 
for a 60-sec sample is consistent with the observed errors. 
Since simulated data showed similar errors, it is likely 
that software error is the major limiting error source and 
that a dramatic improvement in accuracy by perhaps a 
factor of five can be made with improved software. (High 
frequency noise on the data ^ 0.0005 m/sec.) Limiting 
error sources such as ephemeris errors become propor- 
tionately more important for these higher accuracies. 

The real-time runs tend to have larger errors and to 
show different patterns than the postffight runs because 
the a priori statistics used on position were too small in 
proportion to the errors in the a priori position. This was 



Fig. 9. B • Tc error vs hours from Manner IV 
encounter £ — 72HtoF— 5H 
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especially true for B*T (Figs. 9 and 10) since the Y 
a priori position had the largest a priori error and lay 
nearly along B • T. 

Based on the results of runs where A/a = 150 km Vsec^ 
(not plotted), which resulted in target prediction errors 
of over 3000 km for the entire approach phase, it is man- 
datory to solve for the gravitational constant when a large 
uncertainty exists. 

Strong AU and /a solutions are obtained by E— 1 day 
with the AU being known to within a few hundred kilo- 
meters and /A to within 10 km®/sec^ See Figs. 11 and 12. 

As expected, the largest uncertainty in the B plane 
(SMAA) fell along B • R prior to £ — 1 day and then 
gradually shifted to be normal to the B (Fig. 13). Most 
of the shift occurred between £—10 hrs and £—5 hrs. 
The alignment of the SMAA normal to B can be shown 
to follow from the fact that for two-body motion with 
Earth and Mars fixed, there is an unobservable coordi- 
nate of the two-body orbit. This coordinate is the right 
ascension of the ascending node, where the reference 
plane for the orbit is the “plane of the sky” (a plane nor- 
mal to the Earth-Mars line). A change of the ascending 



Fig. 1 1 . GMisfars errors vs hours from 
Mariner IV encounter 


node does not change the range or range rate as seen 
from the Earth, and thus the corresponding coordinate is 
unobservable with DSN tracking data. For Mariner TV 
the incoming and outgoing asymptotes lay very nearly 
along the Earth-Mars line, and thus the “unobservable” 
coordinate was essentially normal to the B and in the 
B plane. This may be seen in Figs. 8-10 and 14-18 
where B and the time of closest approach are deter- 
mined much faster than B • T and B • R which are com- 
ponents of B. The movement of the Earth-probe line 
in inertial space and the gravitational attraction of the 
Sun combine to make the “unobservable” coordinate 
slightly observable. The original work at JPL on this 
relationship was done by Carl Pfeiffer (Ref. 13). 

Comparing the cases where /a and A17 are estimated 
(the ■ cases) with those where only position and velocity 
are estimated, there is a large difference in the SMAA 
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versus time for the time span £ — 30 hr to E — 15 hrs and 
again at £—0.5 hr to £ (Figs. 19 and 20). 

These results indicate the accuracy of the orbit may be 
improved by several hundred kilometers (from 1000 to 
600 km) at £ — 20 hrs if accurate and reliable values for 
the gravitational constant and AU are available before 
the flight. Also, there is a constant diflFerence of 25 to 
50 km between the two SM/A’s from £ — 50 hrs to £ — 1 
hr, which amounts to a large percentage error at £ — 1 
hr (1 km vs 25 km). The £ — 5 H to £ SMI A is essentially 
uncertainty in B so that these results may be compared 
with the B errors in Fig. 16. This figure indicates that 
AU and ft errors are important to the determination of B 
in the £ — 5 to £ period. 

The errors in the components of B (Figs. 9, 10, 14 
and 15) show that the results are insensitive to A17 errors 



Fig. 13. Orientation angle of semiminor axis (SMfA) vs 
hours from Mariner IV encounter (STD 6 
solution with no a priori) 



Fig. 14. B * Rc error vs hours from Moriner IV 
encounter £ — 72Hto£ — 5H 


and to /X errors, provided fi is included in the solution, 
thus confirming the insensitivity of the SMAA to these 
factors. 
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Fig. 15. B • error vs hours from Mariner IV 
encounter £ — 5 H to E + 0.25 H 



Fig. 16. B vector magnitude error vs hours from 
encounter E — 5 H to E + 0.25 H 



Fig. 17. Time of closest approach error vs hours from 
Mariner IV encounter E — 72HtoE — 5H 


5. Implication of the Results for Future Ptanefory 
Encounter Phases Before 1970 

Data weights used for studies and real-time work 
should be no greater than 0.005 m/sec and could be 
considerably smaller depending on what software is avail- 
able and how serious the limiting error sources are. All 
stations will be high resolution S-band for the future 
missions. 

If it appears possible to achieve the necessary accu- 
racy using encounter data with loose a priori, then infor- 
mation from preencounter should be used with extreme 
discretion to avoid the type of errors which occurred on 
Mariner IV, The single precision orbit program appears 
to give stable solutions for the approach phase. 

If there is any serious doubt of the accuracy of the 
planetary gravitational constant, then this constant should 
definitely be included in the solution with fairly loose 
a priori statistics. 
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Fig. 18. Time of closest approach error vs hours from 
fAartner IV encounter E — 5 H to £ + 0.25 H 

The question of whether to estimate AU and perhaps 
the Earth-Mars ephemeris appears to be an open one 
since roughly equivalent answers are obtained in these 
cases or in the case where there is no AU error. This 
conclusion might change if better software were avail- 
able, such as the double precision orbit program sched- 
uled for completion in April 1967, 

This study does not replace a thorough investigation 
of the encounter trajectories for each mission, especially 
with regard to the “unobservable coordinate" (normal 
to B). The variation of accuracy in this coordinate with 
trajectory geometry is not well understood at this time 
and deserves careful study. However, barring any un- 
usual surprises, it may be expected that achieved accu- 
racy on future missions, both to Mars and Venus, will be 
better than that shown here. Improved software will 
be available, ranging data will be taken starting with 
Maritier Venus 67, and for Venus missions the gravita- 
tional bending will be proj)ortionally greater than Mars 
missions. 



Fig. 19. SMA4 and SMI4 v$ hours from Mariner IV 
encounter E — 72HtoE --5H 

N67 U424‘ 

C. Timing Data and the Orbit 
Determination Process at JPL 

P. M. Muller 

1. Introduction 

JPL produces bimonthly summaries of time-data nec- 
essary to the real-time orbit determination program 
process. Samples of the published information are illus- 
trated herein. Fairly adequate predictions of future 
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Fig. 20. SM4A and SMIA vs hours from Mariner IV 
encounter £ — 5 H to E + 0.25 H 


conditions are generated. Those faced with time prob- 
lems will find the published summaries convenient and 
useful. We also recommend the references for background 
on this increasingly important geophysical and astro- 
dynamical topic. 

Ref. 14 documents the role of time in accurate orbit 
determination based on Earth-based radio guidance. 
Necessary time corrections require polynomial approxi- 
mations to real UTl, UTC, and ET. The author has a 
computer program, TPOLY, which accepts U. S. Naval 
Observatory time-bulletin data. The program produces 
a convenient listing of raw data and auxiliary quantities, 
a set of polynomials which fit the data, and a punched 
deck which initializes the orbit determination program 
(ODP) with respect to timehandling. 

The purpose of this article is to publicize this program 
and its regularly published output, which is available to 
all outside users of the JPL orbit determination program 
and any other interested people. A second purpose of 
this article is to outline the current methods of handling 
these timing problems, and to consider the direction of 
future effort. 


2. Program Description 

Input to the TPOLY program is obtained from the 
U.S. Naval Observatory, Washington, D.C. The Observa- 
tory takes no responsibility for our use of their data, 
and this fact must be emphasized. Readers who are not 
familiar with its service should write the Observatory 
for information. Ref. 15 documents its production of data 
bulletins and reprints of important journal articles. The 
N.B.S. in Boulder, Colorado, also participates in these 
activities (Ref. 16). The author takes responsibility for 
the material presented in this article, and for the TPOLY 
program outputs. 

Refs. 14, 17, and 18 give substantial theoretical back- 
ground for this material. Any of them will assist the 
reader. The input data includes the date of observation, 
UTl-UTC, UT2-UTC, and A.l-UTC on that date, and 
ET-UTl for 1958.0. The program assumes that the differ- 
ence ET— atomic time is a constant. 

The program computes various auxiliary quantities 
based upon the input [Fig. 21 (a)]. The title of Fig. 21 
gives the dates identifying the run. Predictions are gen- 
erated for approximately 6 mo ahead, as in Fig. 21 (b). 
Figs. 21 (c)-(e) illustrate the individual polynomial fits 
(least-squares quadratic to UTl, linear for UTC). Fig. 
21 (f) is a summary of the polynomials produced, and 
their equations will be found at the top of the listing. 
Fig. 21 (g) is a daily print of the polynomial time-values 
and their residuals against the observed data, and against 
each other at the interval ends. The sum total of this 
output is duplicated and distributed by the author bi- 
monthly, or on special order at other times. The January 
output contains listings for 1961 through the current date. 
Subsequent output listings are sufficient to make the data 
current. The use of recent polynomials and data is im- 
portant to the more critical applications. Those interested 
in copies of these listings, or of Refs. 14 and 18, should 
write the author. Please ask separately for yearly TPOLY 
outputs, updates of same, copies of the punched poly- 
nomials, and the technical memoranda. Ref. 18 gives a 
detailed user s guide to the program, its printed outputs, 
and generally elaborates on the material presented herein. 

The predictions of future UTl are of special interest. 
Refs. 19 and 20, and their bibliography, are all con- 
cerned with finding a model for the irregularities of the 
Earth s rotation. We are also involved in this problem 
because our ODP must be able to use real-time data. 
It must have some overlap of polynomial coverage to 
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PKOGI^AM Rim TFOLY-1620 VERSION 7» DATED MARCH 1966 
UPDATE UF TPOLY MAY 5, 1966 PREDICTS TO JAN I, 1967 


(a) COMPUTATION OF TPOtY FOR SP(»OP 
INPUT DATA AND AUXILIARY QUANTITIES 

DATE^a/19SQ J«D«ND* I S/50) *E-6 ET-UTl . .ET-UTC A. l-UT 1 _A • 1-UT2 A. l-UTC UT 1-UTC UT2-t)TC UT2-UT1 D IAl^ UTC ) OU> l-UTl) 


*660105 5848 "^2639130.5 505.26720 36.6790 36.6110 4.4290 4.4230 4.3610 -.0680 -.0620 .0060 .0150 .0150 

^660101 5844" 2439126.5 504.92160 OPFSET • 300 STEP « 0 IFIT TO TIME) 

660112 5855 2439137.5 505.87200 36.6900 36.6290 4.4400 4 .4420 4 . 3790 -.06 10 -.0630 -.0020 

.660119 ,.5062_ 2439144.5 506.47600 36.7020 36.6470 4.4520 4.4550 .4,3970 -.0550 -.0580_ -.0030 

660126 5069 2439151.5 507.08160 36.7160 36.6660 4.4660 4.4680 4.4160 -.0500 -.0520 '-.0020 

.660202 .. 5076 2439158.5 507.68640 36,7310 36.6840 4.4810 4.4020 4.4340 -.0470 -.0480 -.0010 

660209 5883 2439165.5 508.29120 36.7475 36.7020 4,4975 4.4970 4.4520 -.0455 ' -.0450 ,0005 

660216 5890_2439172.5 508.89600 36.7650 36,7200 4,5150 4.5 160 _ 4.4700 -.0450 -.0460 ^.0010 

660223 5897 2439179.5 509.50080 36.7830 36.7380 4.5330 4.5330 4.4880 -.0450 ‘-.0450 0.0000 

.660302_.3904 2439186,5 510. 10560 36.8015, 36.7560 4.5515 4.5480 _.4.5060 , -.0455 •0035 

..660301__5903_2439185.5 510.01920, OFFSET ■ 300 STEP - . 0 IFIT TO TIME)_ 

660309 ^5911 2439193. 5_510. 71040 36.0200 36.7740 4.5700 4. 5640_ 4.5240 _-,0460 ^-.0400 

"660316 5910 2439200.5 511.31520 36.0370 36.7930 4.5870 4.5790 4.5430 -.0440 -.0360 

660323 5925 2439207,5 ^11.92000 36,8520 36.8110 4.6020 4.5910 4.5610 _-.0410 -.0300 

660330 ’~5932 ~ 2439214.5 512.52480 36.8690 36.8290 4,6190 4.6050 4.5790 -.0400* -.0260 

660406 5939 2439221.5 513.12960 36.8860 36.0470 4.6360 , 4 .6200 _ 4 .5970 _-.0390 _-.0230 

"660413 * 5946"' 2439228.5 513.73440 36. 9070 36. 8650 4.6570 4.6370 4.6150 -.0420 *"-.0220 

_660420 5953 2439235 . 5..5 14. 33920. 36 . 9270 36.8830 4.6770_4.6540_4,6330 _-.0440 TjOllO 

"660427 5960 2439242.5 514.94400 36.9455 36.9010 4.6955 4.6710 4.6510 -.0445 -.0206 

5.967_2439245.5_5J5..5408O^ 36.9>5Q.3.6.9200 4.7150 ^4..6880_ .4.670p_,-.p450 “•Pl_8p 

.iuS0.501„5P6^.-.2^mA*>_5I5,#2896O..OFESeL.f,...>_O0.STeP -^.,... 0 .„ 

Fig. 21. Program run TPOLY-1620 version 7, dated March 1966 update of TPOLY May 5, 1966 predicts to January 1, 1967 
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(b>. COMPUTATION OF TPOLY FOR SPOOP 


PREDICTS LISTING 


DATE JO/1950 J.O.NO. 

(S/50»*E-6 

6T-UT1 

ET-UTC 

A.l-UTl 

A.l-UTC 

UTl-UTC U(ET-UT1» 

01 ET-UTC » 

0(UT1-UTCI 

666501 5 96A'"2A 39246. 5 

515.28960 

OFFSET 

■ 300 

STEP • 

0 PREDICTS INTERVAL TIME 

— 


5967 2439249.5 

515.54880 “■ 

36.9650 

■ 36.9200 

4.7150 

4.6700 

-.0450 .0195 

.0190 

-.0005 

. 5974.. 2439256.5 

516.15360 

36.9819 

36.9381 

4.7319 

. 4.6881 

-.0438 .0169 

.0181 

• 0011 

5961 2439263.5 

516.75840 

36.9980 

36.9562 

4.7480 

4.7062 

-.0417 ,0160 

.0181 

• P020 

_5988 _ 2439270.5 

517.36320 

37.0135 

36.9744 

4.7635 

4,7244 

-.0390 .0154 

• 0181 

•0026 

5995 2439277.5 

517.96800 

37.0288 

36.9925 

4.7788 

4.7425 

-.0363 .0153 

.0181 

' .002 7' 

6002 2439284.5 

518.57280 

37.0442 

37.0107 

4.7942 

4.7607 

_-.0335 .0153 

.0181 

• 0027 

6009 2439291.5 

519.17760 

37.aS96 

37.0288 

4.8096 

4.7788 

-.0307 ,0153" 

.0181 

»Q02 7 

_6016__2439298.5 

519.78240 

37.0746 

37.0470 

4.6246 

4.7970 

-.0276 .0149 

.0181 

• 0031 

6023 2439305.5 

520.38720 

37.0882 

37.0651 

4.8382 

4.6151 

-.0231 .0136 

• 0181 

•0044 

_6030 2439312.5 

520.99200_ 

37.0991 

37.0832 


4.8332 

-.0158 .0109 

.0181 

• 0072 

660701 6025 2439307.5 

520.56000 

-OilFSET 

■ 300 

STEP • 

0 PREDICTS INTERVAL TIME 
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.. (cl COMfPUTATION . OF. . TPOLY. FOR SPOOP. 



DATA USED IN 

POLY FIT. 




DATE JO/1950" 

(S/50l«E-6 

ET-UTl " 

ET-UTC OFFSET 

STEP 

660101 5844 

'2439126,5 504.92160 


300 

6 

5841 

504.66240 

36.6640 

36.5960 


5848 

505.26720 

36.6790 

36.6110 


5855 

505.87200 

36.6900 

36.6290 


5862 

506.47680 

36.7020 

. 36.6470 


5869 

507.08160 

36.7160 

36.6660 


5876 

507.68640 

36.7310 

36.6840 


5883 

508.29120" 

36.7475 

36.7020 


5890 

508.89600 

36.7650 

36.7200 


5897 

509.50080" 

36.7830 

36.7380 


5904 

510.10560 

36.8015 

36.7560 


660301 5903 

2439185.5 510.01920 


300 

0 


n 


Fig. 21. Cont'd 
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(d) POLYNOMIAL CURVE FIT USING LEAST SQUARES METHOD 
FIT OVER DATA FROM 1 660101 1 58A^ ) TO (660301,5903) FOR INTERVAL (5841 TO 5904) ET-UTl 


TERM COEFFICIENT 

0 •596513845830E^03 

1 ' “ -.2231b53154fl6E-05 

^ j,222368316125E-l4 

STANDARD DEVIATION* .0026 


SEC/50 

ET-UTl 


FIT 

RESIOKAL 

5O46624U0. 

36.6640E-00 


36.6697E-00 

-.0057 

.505267200..... 

36.6790E-00 


36.6782E-00 

.0007 

505872000. 

36.6900E-00 


36.68856-00 

.0014 

506476800. 

36.7020E-00 


36.70036-00 

.0016 

507081600. 

36.7160E-00 


36.7l38€-00 

.0021 

507686400. 

36.7310E-00 


36.7289E-00 

.0020 

508291200. 

36.7475E-00 


36.7456E-00 

.0018 

508096000. 

36.7650E-00 


36.7640E-00 

.0009 

509500600. 

36.7830E-00 


36.78406-00 

-.0010 

510105600. 

36.80156-00 


36.6056E-00 

-.0041 


t^Li^OLYNOMIAL CURVE .FIT. USING LEAST SQUARES METHOD .. 


,F1T OVER OATA FROM 

(660101.5844) TO (660301.5905) 

FOR INTERVAL (5841 TO 

5904) ET-UTC 




TERM 

0 . 

1 

SJANDARD.DEViAIKia* 

COEFFICIENT 
.214480341136E-*-02 
,3000992847756-07 
.0002 

- •- 

- 



. FIT 


SEC/50 

ET-UTC 

RESIDUAL 

505267200. 

36.6110E-00 

36.61106-00 

0.0000 

505872000. 

36.6290E-00 

36.62926-00 

-.0002 

506676800. 

36.6470E-00 

36.6473E-00 

-.0003 

507081600. 

36.6660E-00 

36.6655E-00 

.0004 

507686400. 

36.6840E-00 

36.6836F-00 

.0003 

508291200. 

36.70206-00 

36.70186-00 

.0001 

. . 508896000. 

36.7200E-00 

36.7199E-00 

0.0000 

509500600. 

36.7380E-00 

36.7381E-00 

-.0001 

51010560.0.,. 

36 ♦ 75 60 E -00. 

36.7562E-00. 

-•0002 


Fig. 21. 

permit practical operation. The consensus of these refer- 
ences is that we do not have a satisfactory model at the 
present time. We have, therefore, been forced to use an 
empirical method of prediction. It uses year-old UTl 
data as a priori for the current predictions. Six-month 
averages in UTl slope are used to estimate the expected 
average slope. Ref. 18 gives the details and equations, 
but the model is simply a linear polynomial added to 
year-old data. 


Cont’d. 

The success of this method is due to the regular sea- 
sonal variations in UTl, and to the fact (Ref. 14, p. 109) 
that large, rapid (6-mo to 5-yr) changes in UTl slope are 
fairly rare (generally decades apart). The author, and 
others, are actively investigating these problems from 
our own particular points of view. The plot, Fig. 22, gives 
some typical residuals between predictions and observa- 
tions. It also includes some approximate statistical limits 
on this method. There is no certainty that our published 
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Fig. 22. Typical prediction accuracies of UT 1 


(o) an' example of the best to be expected^ 




MONTHS AFTER LAST DATA 


predictions will fall within these limits. That they have 
done so over the entire period 1962 to the present gives 
some slight confidence in their future value. But Ref. 14 
particularly points out how dangerous these assumptions 
can be. 

N67 \hkl5 

D. Theoretical Basis for the Double 
Precision Orbit Determination 
Program (DPODP); III. Probe 
Trajectory 

T. D. Moyer 


7 . Introduction 

This third article in the series describes the differential 
equations of motion of the probe which are integrated 
numerically in a rectangular coordinate system to give 
the probe ephemeris with ephemeris time (ET) as the in- 
dependent variable. The x axis is directed along the mean 
equinox of 1950.0, the z axis is normal to the mean Earth 
equator of 1950.0, directed north, and the y axis com- 
pletes the right-handed system. The center of integration 
is located at the center of mass of the Sun, the Moon, or 
one of the nine planets. It may be specified as one of these 
bodies, or it may be allowed to change as the probe 
passes through the sphere of influence of a planet (rela- 
tive to the Sun) or the Moon (relative to the Earth). For 
this case, the center of integration will be that body 
within whose sphere of influence the probe lies. At a 
change in center of integration, the position and velocity 
of the probe relative to the old center of integration are 


incremented by the position and velocity, respectively, 
of the old center relative to the new center. The injection 
position and velocity components may be referenced to 
any body (not necessarily the center of integration). The 
injection epoch is specified in the UTC time scale and 
must be transformed to ephemeris time. 

The acceleration of the probe consists of: 

(1) The Newtonian point mass acceleration relative to 
the center of integration. 

(2) The direct oblate acceleration due to a near planet 
or the near Earth and Moon. 

(3) Acceleration due to solar radiation pressure. 

(4) Acceleration due to small forces, such as operation 
of the attitude control system and gas leaks. 

(5) Acceleration due to motor bums. 

The equations for computing each of these acceleration 
terms are given in Section 3. A motor bum of short dura- 
tion or a spring separation may be represented alterna- 
tively as a discontinuity to the probe trajectory. The 
rectangular components of the velocity increment Ar and 
the bum time tb are solved for parameters. The increment 
to the probe position is computed as Ar = % Artb. 

The DPODP has avaflable precomputed position and 
velocity ephemerides for the celestial bodies within the 
solar system. They consist of the heliocentric ephemerides 
of eight planets and the Earth-Moon barycenter, and the 
geocentric lunar ephemeris. The heliocentric ephemerides 
(with the exception of Mercury and Neptune) are ob- 
tained from a separate numerical integration for each 
body with epoch values chosen to obtain a least squares 
fit to source positions. The source positions, which repre- 
sent astronomical observations, are obtained from an 
evaluation of certain general perturbation theories for 
the four iimer planets, from a simultaneous integration 
of the motion of the five outer planets (corrected for 
motion of the iimer planets) by W. J. Eckert, D. Brouwer, 
and G. M. Clemence, and from an evaluation of the 
Brown Improved Lunar Theory. The acceleration due to 
each perturbing body is computed from the source posi- 
tion or position from the fitted ephemeris for the perturb- 
ing body if previously computed. The fitted ephemerides 
are generated one at a time in the following order: 
preliminary Jupiter, Saturn, final Jupiter, Uranus, Pluto, 
Mars, Earth-Moon barycenter, and Venus. The ephem- 
erides for Mercury, Neptune, and the Moon are obtained 
directly from the source positions with velocity obtained 
by numerical differentiation. A description of the JPL 
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(f) POLYNOMIAL SOMMAKY FOR 5POOP 
ET-UTC‘* A+a*T UTl-FT * C^O«T4F*T»T 


OATE 

JO/1950 

A 

B 

D6V. 

c 

0 

6 

OEV. 

FROM 

TO 

FROM TO 








640101““ 

640401 

"(5113»5;?04> 

"“.284440206+02 

. 149558026-07 

.0002 

‘ <-.55H735A9E + 02 

.116424566-06 

~ .157313936-15 

".0004 “ 

640401 

640710 

(5204,5304) 

.285354756+02 

.149750586-07 

0.0000 

. lA5051H2E+03 

-.773140126-06 

.827272316-15 

.0027 

640710 

640901 

(5304,5357) 

.285407336+02 

. 149636306-07 

0.0000 

-.211380996+03 

.777932376-06 

-.860134506-15 

.0006 

640901 

641101 

(5357,5418) 

.285346666+02 

.15192921E-07 

.0002 

-.236831746+03 

.889063766-06 

-.981435526-15 

• 0021 

641101 

650101 

(5416,5479) 

.285900956+02 

. 150740386-07 

.0002 

-.915242566+02 

.26298UAE-06 

-,30 7090eoE-l!> 

.0005 

650101 

650301 

(5479,5538) 

.206636226+02 

.15131A50E-07 

.0004 

.129556566+03 

-.671701546-06 

.680820416-15 

.0009 

650301 

650501 

(5530,5599 ) 

.287966^ 7fc-*-02 

• 15061 31AF-07 

.000 2 

342462106+03 

.130375866-05 

13fl6072/>6-14 

.0014*^ 

650501 

650701 

(5599,5660) 

.288841606+02 

.148609156-07 

0.0000 

.209064416+03 

-.984533566-06 

.987360376-15 

.0012 

650701 

650901 

(5660,5722) 

.288740816+02 

.151060466-07 

.0002 

-.331879496+03 

.122094376-05 

-.126053046-14 

.0012 

650901 

651101 

(5722,5783) 

.290479036+02 

.149560926-07 

.0002 

-.357600546+03 

.132290786-05 

-.136153326-14 

.0011 

"“651101“ 

660101 

(5783,5044) 

■ .290407926+02 

. 149711026-07. 

.0002 

' .176284746+03 

-.818121456-06 

.785002176-15 

.0011 

660101 

660301 

(5844,5903) 

.214480346+02 

.300099286-07 

.0002 

-.596513846+03 

.223155316-05 

-.222368316-14 

.0026 

660301 

‘660501" 

(5903,5964) 

.214441236+02 ' 

.300174646-07 

.0003 

196676796 + 03 

‘ “ .653383826-06 

"“-.666465516-15 

' .0014 ■' 

660501 

660701 

(5964,6025) 

.214201386+02 

.300643276-07 

.0002 

.216112506+03 

-.951989376-06 

.894391596-15 

.0010 

“ “660701~* 

660901 

(6025,6087) 

.214535496+02 

.299999756-07 

0.0000 

-.292514586+03 

.996156916-06 

-.971049306-15 

.0008 ■ 

660901 

661101 

(60&7,6148) 

.214535376+02 

.299999966-07 

0.0000 

-.459732196+03 

.163056116-05 

-.157276006-14 

• 0012 

661101 ‘ 

670101 

(6148,6209) 

".214535436+02 

.299999056-07 

0.0000 

.108124246+03 

-.514352616-06 

.452677646-15 

.0005 

670101 750101 (6209. 9 1,3U._ 

_.214535356+02 

.300000006-07 

0.0000 

-.214346466+02 

-.300000046-07 

. .320906866-23 

0.0000 


(g) FIT OF POLYNOMIAL OVER PERIOD 660101 TO 660301 INCLUSIVE 
PRtVI^Tjrii'FTsEin^ 150~ THIS OFFSET « 300 NEXT OFFSET » 300 “(PARTS 7 n KF + 1()) 


TYPE 

DATE 

JD/50 

. 0 . NO • 

(S/50) *6-6 

6T-UT1 

ET-UTC 

A.l-UTl 

A. 1-UTC 

UTl-IITC R(ET-UTn RIET-UTC) 

' POLYl 

660101 

56AA 

2439126.5 

" 504.92160 

36.6731 

36.6006 

4.4231 

4.3506 

-.0724 

POLYL. 


5845- 

_ 2439127.5 

505. 00800. 

-.36.6744 

... 36.6032 

. 4.4244 

.4.3532 

-.0711 

POLYl 


5846 

2439128.5 

505.09440 

36.6756 

36.6058 

4.4256 

4.3550 

-.0697 

^-POLYl. 


...5847 .. 

..2439129.5 

- .505.180B0 

36.6769 

. 36.6084 

..4.4269 

4. 3584 , 

-.0684 

POLYl 


5848 

2439130.5 

505.26720 

36.6782 

36.6110 

4.4282 

4.3610 

-.0672 

DATA 


5848 

2439130.5 

505,26720 

36.6790 

36.6110 

4.4290 

4.3610 

-.0600 .0007 0.0000 

POLYl 


5849 

2439131.5 

505.35360 

36.6796 

36.6136 

4.4296 

4.3636 

-.0659 

POLYl 


.5850 

..2439132.5 

505.44000 

36.6810 

36.6162 

4.4310 

4,3662 

-.0648 

POLYl 


5851 

2439133.5 

505.52640 

36.6824 

36.6188 

4.4324 

4.3608 

-.0636 

POLYl 


5852 

2439134.5 

505.61280 

36.6839 

36.6214 

4.4339 

4.3714 

-.0625 

POLYl 


5853 

2439135.5 

505.69920 

36.6854 

36.6240 

4.4354 

4.3740 

-.0614 

POLYL 


5854_ 

_2439136.5_ 

505.78560 

_ 36.6869 

36.6266 

4.4369 

4.3766 

-.0603 

POLYl 


5855 

2439137.5 

505.87200 

36.6885 

36.6292 

4.4385 

4.3792 

-.0593 

DATA 


5855 " 

~ 2439137.5" 

”505.07200 

36.6900 

36.6290“ 

“ 4 .4400 ^ 

" 4.3790 ’ 

.06 10 ".“0014 .0002 

"pOLYl 


5898 

2439180.5 

509.58720 

36.7869 

36.7407 

4.5369 

4.4907 

-.0462 

POLYl 


5899 

.. 2439181.5 

^ . 509.67360 

36.7900 

_ 36.7433 

4.5400 . 

, 4.4933 . 

-.0467 

POLYl 


5900 

2439182.5 

509.76000 

36.7930 

36.7458 

4.5430 

4.4958 

-.0471 

..... POLYl 


.5901 . 

2439183.5 

. 509.84640 

36.7961 

36.7484 

4.5461 

4.4984 

-.0476 

POLYl 


5902 

2439184.5 

509.93280 

36.7992 

36,7510 

4.5492 

4.5010 

-.0481 

....POLYL 


. 5903 . 

2439185.5 

510.01920 

. 36.8024 
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ephemeris tapes and the method by which they were 
generated is given in Refs. 21 and 22. 

For three of these precomputed ephemerides, the 
orbital elements of an osculating conic at some epoch are 
solve for parameters. The ephemeris correction consists 
of the diflFerence between the solve for conic and the 
osculating conic. Section 2 gives the formulas for cor- 
recting the basic ephemerides and for computing the 
corrected position and velocity of the Sun, Moon, or a 
planet relative to any other of these bodies. These quan- 
tities are required (1) to compute the above mentioned 
acceleration terms, (2) at phase changes, and (3) at the 
injection epoch, if the injection position and velocity are 
referred to a body other than the center of integration. 

The probe ephemeris is represented by the three nu- 
merically integrated sum and diflFerence (SAD) arrays. 
These arrays may be interpolated at any value of ET 
to give r, r, r, and r of the probe relative to the center 
of integration. 


(E) 

r (l.u.) = position vector of Moon relative to Earth 
in dimensionless units (l.u,), obtained directly 
from the lunar ephemeris 

Ae = conversion factor from astronomical units (a.u.) 
to kilometers 

Rfi = conversion factor from dimensionless lunar units 
(l.u.), to kilometers 

E = osculating two-body orbital elements for helio- 
centric ephemeris of planet or Earth-Moon bary- 
center, or geocentric lunar ephemeris 

0r/0E — partial derivatives of position on osculating 
conic with respect to orbital elements (Ref. 23, 
p.241) 

aE = solve for corrections to osculating orbital ele- 
ments 

The position of the Moon relative to the Earth-Moon 
barycenter B and of B relative to the Earth is given by 


2. Planetary and Lunar Ephemerides 

The corrected heliocentric position of a planet or 
Earth-Moon barycenter is given by 

0 r (S) 

*• (km) = (a.u.)AB + ■ AE„ (1) 

Similarly the corrected geocentric position of the Moon 
is given by 

0r 

Each of these vectors has rectangular components re- 
ferred to the mean Earth equator and equinox of 1950.0. 
The X axis is along the mean equinox of 1950.0, the z 
axis is normal to the mean Earth equator of 1950.0 
directed north, and the y axis completes the right-handed 
system. In Eqs. (1) and (2) above, 

(S) 

r ^ ' (km) = corrected position vector of planet or 
^ Earth-Moon barycenter relative to Sun in units 
of kilometers 


r 


(S) 

P 


(a.u.) = position vector of planet or Earth-Moon 
barycenter relative to Sim in astronomical units 
(a.u.), obtained directly from the ephemeris 


r 


(E) 

M 


(km) = corrected position vector of Moon rela- 
tive to Earth in units of kilometers 


where 


(B)_ (£) 

M 1 + M 

(3) 

(E) _ 1 (E) 

B 1 + M 

(4) 

11 

(5) 


where 

fiE = gravitational constant of Earth, km^/sec^ 
fiM = gravitational constant of Moon, km^/sec^ 


The following table lists sums of the position vectors 
Eqs. (1) to (4) to give the position vectors of the Earth 
(E), Moon (M), Sun (S), or planet (p) (other than the 
Earth) relative to the center of integration, with each of 
these bodies, respectively, as center. 


Earth = Center 


r 


(E) 

M 


= r 


(E) 

M 


r (£) _ r (^) _ r 

S B B 



B 


(S) 

B 


+ r 


(S) 

P 


f 
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Moon = Center 


where 


.m ^_(E) 

M 

(B) 

M 

(B) 

M 


.(M) 

S 

.(M) _ 


= — r 


— r 


(S) 

B 


= — r 




Sun = Center 

E B B 

,(S) ^^(S) (B) 

' M B M 

r(S) 

P P 


Planet = Center 

r(B) =-,(S) +,(S) _r(E) 

E 'p^ B B 

r (B) - - + r + r 

^ M - ^ p B M 


(P) = _ r (S) 


r(P) 

Pi P Pi 


All of the formulas in this section also apply when r is 
replaced by f, where the dot denotes diflFerentiation with 
respect to ephemeris time. 

The solve for parameters which affect the relative posi- 
tion and velocity between two celestial bodies are the: 
scaling factor Ae for the heliocentric ephemerides; the 
scaling factor Re for the lunar ephemeris; corrections to 
osculating orbital elements Ae for any three of the eph- 
emerides; and the gravitational constants of the Earth 
and Moon, fXE and These are known as reference 
parameters. 


/i = ftij + fip for the heliocentric ephemeris of a planet 

= /AS H- At£? + /am for the heliocentric ephemeris of the 
Earth-Moon barycenter 

= /ae + /am for the geocentric lunar ephemeris 


r = f - j!iL 


(7) 


where 


r = 


r • r 


(8) 


The values of r and r are computed from Eqs. (3) 
M D 

. JE) ..(S) ..(S) ..(E) ..(B) ..(E) 

and(4)usingry. Given >■ b ’ M ’ ’’ M ’ ’’ B ’ * 

relative r, *r between two celestial bodies is computed 
from the table of vector sums where r is replaced by r,'r. 


3. Probe Acceleration 

The equations for computing each term of the total 
probe acceleration relative to the center of integration 
are given in the following subsections. 

{ a. Newtonian point mass acceleration, 

r= ^;^i[h(r-r,) + h(rO] (9) 

i 

where i is summed over each planet, the Sun, and the 
Moon, including that body which is the center of integra- 
tion, and 


fii = gravitational constant of body t, kmVsec^ 

r = position of probe relative to center of integra- 
tion in 1950.0 Earth equatorial rectangular 
coordinates 

r . = position of body i relative to center of integra- 
tion in 1950.0 rectangular coordinates 

h(r) = vector function of r = — r/r^ 


In order to compute doppler observables, the relative 
acceleration r and jerk r between two celestial bodies 
are required. Given corrected r, f of a planet or Earth- 
Moon barycenter relative to the Sun or Moon relative to 
the Earth, r and r are computed from two-body formulas: 

'^= (6) 


i b. Oblate body acceleration. The acceleration of the 
probe due to the oblateness of a planet or the Moon is 
the direct acceleration of the probe due to the oblate- 
ness. The indirect acceleration of the center of integration 
due to the oblateness of a perturbing body is ignored in 
general. However, for the case where the Earth is the 
center of integration and the Moon is a perturbing body 
or vice-versa, an expression has been derived by F. Sturms 
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(SFS 37-29, Vol. IV, pp. 1-6) for the indirect acceleration, 
which accounts approximately for the oblateness of each 
body. 

The direct acceleration of the probe due to oblateness 
of a body is derived from the generalized potential func- 
tion (Ref. 24, pp. 173-174) for that body: 

X {Cn,m COS m\ + Sn,m sin mX)^ (10) 

where 

ft == gravitational constant of body, 
km^/sec^ 

radius, latitude, and longitude (posi- 
tive east of prime meridian) of probe 
relative to body 

dp = mean equatorial radius of body (an 
adopted constant used for 17) 

(sin <f>) = associated Legendre function of the 
first kind. The argument sin <f> will be 
omitted here 

Cn,m and m = numerical coefficients (tesseral har- 
monic coefficients) 


The inertial acceleration of the probe is computed in 
a rectangular coordinate system {x'y'z') with the x' axis 
directed outward along the instantaneous radius to the 
probe, the y' axis directed east, and the z' axis directed 
north. Fig. 23 shows these axes relative to body-fixed 
axes Xbt/ftZb, where Xf, is along the intersection of the prime 
meridian and equator of the body, z& is directed north 
along the axis of rotation of the body, and {/& completes 
the right-handed system. The transformation from body- 
fixed coordinates r& = (x&, {/&, z&) to r' = (x', t/', z') co- 
ordinates is given by 


where 
R = 



r' = Rrft 


(15) 

cos <fi cos A 

cos <l> sin A 

sin </> 


— sin A 

cos A 

0 

(16) 

— sin <l> cos A 

-'sin<|>sin A 

cos (f> 



The position of the probe relative to the body in rec- 
tangular coordinates referred to the mean Earth equator 
and equinox of 1950.0 is given by r — where r and ri 
are defined in Section a. Transformation from these iner- 
tial coordinates to body-fixed coordinates is defined as 


n = T^(r - u) (17) 

The overall transformation from (r — ri) to r' is thus 

r' = G(r - rO = RT^{r - n) (18) 


The associated Legendre function is defined by 


where 


pm = ^ 


d{sm <f>y 


Pn 


( 11 ) 


P„ = Legendre polynomial of degree n in sin 

The zonal harmonic coefficient Jn is defined as 

7n=-C,.o (12) 

Eq. (10) may be written in three terms corresponding to 
the potential of a point mass, zonal harmonics, and tes- 
seral harmonics 



X P {Cn,m cos mA + Sn,m sin m\) (13) 

with each term defined as 


The inverse transformation is 

(r - n) = GV = TRV (19) 



f/ = ^ + U{J) + [7(C, S) 


(14) 


Fig. 23. x^ y' and z' axes relative to fixed 
body axes x^y^Zh 
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The transformation T is specified in the DPODP for the 
Earth, the Moon, and Mars. Using body-fixed coordinates 
from Eq. (17), trigonometric functions of and A are 
given by 

sin<t> = ^ (20) 


cos = 


r 


( 21 ) 


where Ui has a maximum value of 15, and tiz has a maxi- 
mum value of 8. 


The Legendre polynomial is computed reciu*sively 
from 


P„ = ^ sin < 1 , P„-x - (-^^IT^) (29) 

starting with 


sin A — 


V*6 + vl 


(22) 


+ vl 

An expression will be developed for the inertial accel- 
eration of the probe due to oblateness for any body with 
rectangular components along the instantaneous direc- 
tions of x' i/', and zf axes, denoted 'r'. This acceleration 
can be broken down into r' (/) due to the zonal harmonics 
and r' (C, S) due to the tesseral harmonics. With these 
terms, the contribution to the probe acceleration r rela- 
tive to the center of integration in Earth equatorial rec- 
tangular coordinates due to the oblateness of any body is 


Po=l 


(30) 


Fi = sin 

The first derivative of Fn with respect to (wrt) sin <^, 
denoted F^ , is given by 

F; = sin^F;_, +nF„_i (31) 

starting with 

f; = 1 (32) 

The function sec <#> F^ is computed by first generating 

sec </) P” = (2m - 1) cos <#> [sec <f, P^\ ] (33) 


r = = QT (j) + -if (c, S)} 


The components of if' are computed from 


x' = 


dU 

dr 


1 dU 

r cos <f> Sa 

1 dU 

r d<l> 


(23) 

(24) 

(25) 

(26) 


where the point mass term of U has been accounted for 
in section a., and is ignored here. Carrying out these 
derivatives gives 





- (n + 1)P„ - 
0 

-cos<^p;_ 


(27) 


starting with 


sec </. P[ = 1 (34) 

and then generating 

‘ 4> P: = (^"J^^ ) sin <!> [sec P” , ] 

/n + m— l\n «« 1 


sec < 


(35) 


For each value of m between 1 and na, n is varied from 
m + 1 to fig. The general term F^ is zero if b > a. 


The function cos P^\ where P^' is the derivative of 
P^ wrt sin is computed from 


■^(c.s) = ^i:t(T)” 

' n=l m=i \ / 


X 


" — (n + 1)P” {Cn,m COS m\ + S„,„ sin mx) 
m sec </> P“ { — sin mx + cos mX} 
cos <f> Pf {C„,m cos mX + S„,m sin mX} 


(28) 


cos(^P™' = - nsin.^ [sec<^P“] + (n + m) [sec P“_,] 

(36) 


\ c. Acceleration due to solar radiation pressure and 
attitude control. The acceleration of the probe due to 
solar radiation pressure and small forces such as gas 
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leaks from the attitude control system, noncoupled attitude control jets, etc., is represented by 

r = { [flr + hr(t — T^ci) + Cr(t — Taci)*] [«(# ~ T aC\) ~ U(t ~ T AC2)] 

+ ^Gr] U*(t - VsP 

{[^X + bx(t — TaCi) + Cx(t — TaCiY^ [u{t — TaCi) u{t — TaC 2 )] 

+ Aa, + [G, + G; (EPS) + aG.] u*(t - Tsbp)) X* 

+ {[^y T^ac^ *i~ Cy(t Txci)^] ” Taci) ~ ^,402)] 

+ [G, + G; (EPS) + aG^] tt*(f - Ts^p)) Y* 


(37) 


The terms in this equation are defined as 


m = instantaneous mass of probe, kg 


U^-p = unit vector from Sun to probe 

X*,Y^ = unit vectors along spacecraft x and 
y axes (X* X ¥♦ = Vgp) 

ai,bi,Ci where i = r,x, or y = solve for coeffi- 
cients for acceleration polynomials, 
km/sec^ km/sec®, km/sec^ 

t = ephemeris time 

Taci,Tac 2 = epochs at which the attitude control 
(AC) acceleration polynomials are 
turned on and off, respectively. The 
epochs may be specified in the UTC, 
S.T., or A.1 time scales 

u(t Taci) ~ 1 for t ^ T^in» 0 for t <C Taci 
TaCI TaC2 

AarAo^xAciy = input acceleration (not solve for), 

km/sec^. The value for each A«i will 
be obtained by linear interpolation 
between input points specified in any 
time scale 


where 


C, 


JAg2 ^ 

c 


1 km^ 


= 1.031 X 10« 


km^kg 

sec^m^ 


J — solar radiation constant 

= 1.383 X 10^ watts/m^ 

= 1.383 X lO-’ kg/sec^ 

Afi = 1.496 X 10« km 

c = 2.997,925 X 10« km/sec 

Ap = nominal area of spacecraft projected 
onto plane normal to Sun-probe line, 


fsp = distance from Sun to probe, km 

Tsrp = epoch at which acceleration due to 
solar radiation pressure is turned on 
The epoch may be specified in the 
UTC, S.T., or A.l time scales 

u*{t—TsRp) = 1 for f > Tsrp if spacecraft in sunlight, 
0 for ^ < Tsrp ot if spacecraft in 
shadow 

Gr = solve for effective area for acceleration 
of spacecraft in radial direction due to 
solar radiation pressure, divided by 
nominal area Ap 

Gx — solve for effective area for acceleration 
of spacecraft in the direction of its 
positive X axis (along X* vector) 
divided by Ap 

Gy = solve for effective area for acceleration 
of spacecraft in the direction of its 
positive y axis (along Y* vector) 
divided by Ap 

G' G' G' = solve for derivatives of Gr, Gx, Gy with 
respect to Earth-probe Sun angle, EPS 

EPS = Earth-probe Sun angle, radians 

AGr,AG;r,AGy = increments to Gr, Gx, and Gy obtained 
by linear interpolation of input points 
specified in any time scale 


The unit Sun-probe vector Vsp is computed from 


VsP 


r - rf > 

l|r-rr II 


(38) 
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where 

r = rectangular coordinates of probe relative to center 
of integration, referred to mean Earth equator 
and equinox of 1950.0 

~ rectangular coordinates of Sun relative to center 
of integration C, referred to mean Earth equator 
and equinox of 1950.0 

The spacecraft X* and Y* unit vectors are obtained as a 
rotation of the tangential T and normal N vectors through 
the angle K: 



The angle K is an input (nonsolve for) constant. Com- 
putation of the unit vectors T and N requires the unit 
vector U/j 


U/? — unit vector from spacecraft to reference body 
which orients the spacecraft about the roll axis 
(Sun-probe line). The reference body may be 
a star, a planet, or the Moon. 

If the reference body is a star 


cos 8 

cos 

a 

cos 8 

cos 

a 

sin 

8 



(40) 




rf» -r 

-••I 


(41) 


where 


j.^c) ~ rectangular coordinates of reference body B 
relative to center of integration C, referred to 
mean Earth equator and equinox of 1950.0 

The unit normal vector N (normal to Sun-probe 
reference body plane) is computed from 

Ur X Us/> 


N 


U« X UsP I 


(42) 


The unit tangential vector T (tangent to Sun-probe- 
reference body plane) is 

T = N X Vsp (43) 


The EPS angle may be computed from 


cos EPS = — U^p • U' 


0 < EPS < 180° 


(44) 


where 

U^ is computed from Eq. (41) using B = Earth, 

\ d. Acceleration due to motor burn. The acceleration 
of the probe in 1950.0 Earth equatorial rectangular coor- 
dinates due to a motor burn is given by 

r - - To) - /x(f - Tf)] (45) 

where 

a = magnitude of r vs time 
U = unit vector in direction of r vs time 

To — effective start time of motor, the ET value 
of the solve for UTC epoch, To (UTC) 

Tf = effective stop time of motor, ET 
t = ephemeris time 

The efiFective stop time Tf is given by 

T, = To + T (46) 

where 

T = solve for burn time of motor, ET 


where the right ascension a and declination 8 are referred 
to the mean Earth equator and equinox of 1950.0. If the 
reference body B is a planet or the Moon (normally the 
Earth), 


The acceleration magnitude a is given by 

F{t) ^ Eo + Fit + F^t^ + F,t^ -f F,t^ ^ 

“ ~ mo - Mot - V2Mrr- - - ViUr 

(47) 
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where 

F{t) = magnitude of thrust at time t. The 
polynomial coefficients of F{t) are solve 
for parameters 

m{t) = spacecraft mass at time f 

C — 0.001 for F in newtons and m in kg. 
For F in lb and m in Ibm, 

C = 0.009,806,65 

trio — spacecraft mass at To (ET) 

= polynomial coefficients for propellant 
mass flow rate (positive) at time t, 
iH{t) = Uo + Mst" 

(not solve for parameters) 

t — E.T. — To (ET), seconds 
The unit vector U in the direction of thrust is given by 

(cos 8 cos 

> = ^ cos S sin a > (48) 

U,} \ sins ) 

where 

a — right ascension of U 
S — declination of U 
given by 

a = tto + Ui# — 02 ^“ + a^t'^ -I- ait* (49) 

S = So + Sif - Sst- + Ss^" -1- Sit* (50) 

where the polynomial coefficients of Eqs. (49) and (50) 
are solve for parameters. 


N67 14426 

E. Theoretical Basis for the Double 
Precision Orbit Determination 
Program (DPODP): IV. 

Light Time Solution 

T. D. Moyer 

7. Introduction 

This fourth article in the series gives the formulation 
for the solution of the light time problem, which is the 
first step in the computation of all observable quantities. 
An electromagnetic signal is transmitted from a tracking 
station on Earth at time ti. This signal is received by the 


probe and retransmitted at time arriving at the same 
or different tracking station on Earth at time h. Alterna- 
tively, the signal may be transmitted directly by the 
probe at time tz. All observables are related to charac- 
teristics of this electromagnetic radiation, i.e., the angle 
of the incoming ray, the ratio of received to transmitted 
frequency, or the round trip transit time. The transmit- 
ting and receiving station(s) and the probe are referred 
to as direct participants, and ti, tz, and ts are the epochs 
of participation of each participant. The solution to the 
light time problem consists of the heliocentric position, 
velocity, acceleration, and jerk of each direct participant 
evaluated at its epoch of participation. The next articles 
in the series will give the formulation for computing 
doppler, range, and angular observables, starting with 
the solution to the light time problem. 

The solution to the light time problem is obtained by 
solving the light time equation for each leg of the path 
of electromagnetic radiation from the transmitting to 
the receiving station. The light time equation relates the 
light time between two points to the heliocentric posi- 
tions of each of the two participants evaluated at their 
epochs of participation. Starting with the known recep- 
tion time fa, the light time equation is solved by an 
iterative technique for the down leg of the light path 
to give the epoch of participation for the probe, tz. 
Given tz, the light time equation is solved iteratively for 
the up leg of the light path to give the transmission 
time fi. 

The light time equation is derived in Section 2 and 
accounts for the effects of general relativity. The method 
of computing the heliocentric position vector of each 
participant is given in Section 3. If a participant is a 
station, the position of the station relative to the center 
of the celestial body on which it is located must be com- 
puted. The formulation for doing so is given in Section 4. 

2. Light Time Equation 

Light moves along a geodesic path in the four- 
dimensional non-Euclidean geometry of space time. The 
geometry is embodied in the coefficients of the expres- 
sion for the distance or interval ds between two points 
(Ref. 25, p. 93) 

ds^:^(l - -^+ -^)c-dt^ 

\ c-r c*r^ j 

- (l + (df= + r- deF + sin* e d4>^) (1) 
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where 

r,e,4> - Schwarzschild isotropic spherical 

coordinates of a point in a heliocentric 
nonrotating frame of reference 

t = uniform coordinate time in heliocentric 
frame of reference (ephemeris time) 

fji = gravitational constant of Sun, 

1.32715 X 10“ kmVsec^ 

c = speed of light at infinite distance from 
Sun, 299,792.5 km/sec 

ds = interval between two points with 

differentials dr, dd, dij>, and dt in the 
space-time coordinates 

The departure from pseudo-Euclidean geometry in Eq. (1) 
is due to the mass of the Sun. The effects of the mass of 
each planet and the Moon have been ignored. 

The square of the coordinate speed of light s is 

For the motion of light, the interval ds between two 
points on the light path is zero (Ref. 1, p. 93). Setting 
d« = 0 in Eq. (1), substituting Eq. (2), and ignoring 
1/c^ terms gives 



If light moved along a straight line at speed s relative to 
the Sun, the light time equation would be of the form 

) ( 4 ) 

where ti is the time of emission of light from point 1; 
t 2 is the time of reception of light at point 2; and fi 2 is 
the distance between point 1 at time and point 2 at 
time # 2 , measured in a heliocentric nonrotating frame of 
reference. The 1/c® term is due to the departure of the 
speed of light from the constant c. 

For light passing by the Sun on a curved geodesic path 
with minimum radius R, the angle between the incoming 
and outgoing asymptotes is given by 4 ^/c^R, which has 
the well known maximum value of 1'75 when R is set 
equal to the radius of the Sun, 700,000 km. The maxi- 
mum angle between the straight line path between two 
points and the actual curved path is 4fx/c^R, If the 
nominal length of the light path between two points is L, 


the difference in length aL between the straight line 
path and the cmved geodesic path will be less than 



which is of order 1/c^. Thus, the bending of light pro- 
duces a term of order 1/c® in the light time equation 
[Eq. (4)]. If all terms of order 1/c® and greater are 
ignored in the derivation of the light time equation, the 
assumption of straight line motion at speed s given by 
Eq. (3) is valid. 

Fig. 24 shows the geometry for straight line motion of 
light relative to the Sun S. Light is emitted from point 1 
at ephemeris time fi and arrives at point 2 at ephemeris 
time ^ 2 . The angles and P 2 are computed from 


cos 81 = — 

. £11 

( 6 ) 

ri 

fl2 

cos 82 = — 

.£ii 

( 7 ) 

f2 

fl2 

where [r^ = (ti)], [rg = 

rf' (#.)] 

= heliocentric 


position vectors of point 2 at ephemeris time ^2 and point 
1 at ephemeris time ti, with components referred to a 
nonrotating coordinate system. In the DPODP this coor- 
dinate system is referred to the mean Earth equator and 



Fig. 24. Geometry for straight line motion 
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equinox of 1950.0. The x axis is along the mean equinox 
of 1950.0, the z axis is normal to the mean Earth equator 
of 1950.0 directed North, and the y axis completes the 
right-handed rectangular coordinate system. 

Ti2 = T2 Ti (8) 

fl2 ~ I I **12 I I (9) 

Noting the definitions of Eqs. (6) and (7), the angles 
and /?2 are in quadrants I and II and can be determined 
from their cosines. The time for light to travel from 
point 1 to point 2 at speed s given by Eq. (3) is 


^2 ti — 

From Fig. 24, 




da = 


dr 


r dr 


cosfi - R2 

Substituting Eq. (11) into Eq. (10) gives 


^2 ti 


r dr 


^ I yjr^ ” R\ 


: + 


^ r dr 


(11) 


( 12 ) 


Evaluating the integrals gives 


^2 tl — 


i r^^prRT _ ^ In 

C L'' 2 ^ ^ “ j + Vr" - Ro . 


The 1/c term is simply fn/c and is r cos /3. Thus, 


But fj sin J02 = Tl sin = R„ 


ti tx 


Txi , ^ . r r2_(lH2_cos_^“ 

c ” fi (1 + cos /3i)_ 


(13) 


(14) 

(15) 


Substituting raAi from Eq. (15) into Eq. (14) gives the final expression for the light time equation 


ti - tx = lH. + %ln 
c 


r tanl^/3i ~l 
[_ tan % j32 J 


(16) 


where Txi is evaluated from Eqs. (8) and (9), and ^x and ^2, which are in quadrants I and II, are evaluated from Eqs. 
(6) and (7). Note that px is always greater than ^ 2 - Since the right side of Eq. (16) is a function of #2 and tx, the equa- 
tion must be solved by an iterative procedure. 


The 1/c® term in Eq. (16) is indeterminate for light traveling radially to or from the Sun. For light traveling 
radially from the Sun, jSx = ^2 = 0° and 


tan % _ sin fix . 1 + cos ^2 _ 1 + cos ^2 f2 _ ^2 

tan % ^82 1 + cos / 3 x sin ^2 1 -I- cos /Ji tx r, 

where f2 > ri. For light traveling radially toward the Sun, j8i = J82 = 180° and 


(17) 


tan%y3i _ 1 — cos )8i ^ sinj32 _ 1 — cos fix fi _ fi 
tan Vifi 2 ~ sin ^x 1 — cos /?2 1 — cos jSj f2 


(18) 


where Tx > T 2 . Hence, for light traveling radially to or from the Sun with r larger and Temaiur equal to the larger and 
smaller values of Tx and f2, 


^2 


C Cr "^STnaller 


(19) 


In Eq. (16), if = 90°, fj = 0.7 X 10* km (radius of Sun), r2 = 1 a.u. = 1.5 X 10* km, ^2 16' and 



sec) In 


1 

1/430 


= 6,1 X 10-* sec 
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In Eq. (19), if fi = 0.7 X 10® km (radius of Sun) and = 1 a.u. = 1.5 X 10® km 

^ In ^ 215 = 5.4 X 10-® sec 

c® Tsmaller 


3. Heliocentric Position Vector of Each Participant 

The light time equation is solved for each leg of the 
path of electromagnetic radiation from the transmitter 
to the receiver in reverse order starting with the known 
reception time at the receiving station. The solution 
consists of the epoch of participation ti and the helio- 
centric position, velocity, acceleration, and jerk of each 
participant evaluated at its epoch of participation. The 
components of each vector are referred to the 1950.0 
Earth equatorial rectangular coordinate system. 


r, r, and Ir. The formulation for computing the helio- 
centric position, velocity, acceleration, and jerk of a 
celestial body is given on pp.25 and 26 of this summary. 
The position, velocity, acceleration and jerk of the probe 
relative to the center of integration is obtained by in- 
terpolation of the probe ephemeris sum and difiFerence 
arrays. The formulation for computing the position, 
velocity, acceleration, and jerk of a station in 1950.0 
Earth equatorial rectangular coordinates centered at the 
body on which the station is located is given in Section 4 
of this article. 


The light time equation (Eqs. 16 or 19) relates the 
light time in ephemeris time (ET) for a given leg of light 
path to the heliocentric position vectors of the two par- 
ticipants evaluated at their ET epochs of participation. 
However, for a receiving station on Earth, the known 
epoch of participation is in station time (ST) and must be 
converted to ET. Furthermore, to compute the helio- 
centric position vectors of certain participants requires 
the epoch of participation in time scales other than ET, 
i.e., UTl for an Earth station. Thus, the time transforma- 
tions given in the second article [SPS 37-39, Vol. Ill, 
pp. 36-n38, Eqs, (3), (7), (8), and (9)] must be used exten- 
sively in the solution of the light time problem. 


For a station on the Earth (E), the position of the 
station relative to the Sun (S) at its epoch of participa- 
tion ti is computed from 

r (ti, ET) = r g (ti, ET, UTl) + r (|) (ti, ET) 

( 20 ) 

For a station or landed spacecraft on another celestial 
body B, 

( 21 ) 

For the probe with center of integration C, 

r (ti, ET) = r (ti, ET) + r (ti, ET) 

( 22 ) 

Each of these vector sums in 1950.0 Earth equatorial 
rectangular coordinates applies also for r replaced by 


4. Body-Centered Station Position, Velocity, 
Acceleration, and Jerk in 1950.0 
Rectangular Coordinates 

a. Introduction. This section gives the formulation for 
computing the position, velocity, acceleration, and jerk 
of a station in a 1950,0 Earth equatorial rectangular 
coordinate system centered at the body on which the 
station is located. 

A fixed station may be located on the Earth, the Moon, 
or Mars with spherical or cylindrical coordinates. Al- 
though not discussed here, the DPODP can also account 
for the motion of a moving ship station on Earth, 

For a fixed Earth station, two models are available 
to represent the station position. The simplest is to repre- 
sent the station position by constant values of the co- 
ordinates relative to the true axis of rotation, equator, 
and prime meridian. However, due to wandering of the 
pole (axis of rotation) relative to the Earth, the actual 
position of a fixed station relative to the true pole, equa- 
tor, and prime meridian is a variable quantity. Therefore, 
the alternate method of representing the station position 
is to use constant values of station coordinates relative to 
the mean pole, equator, and prime meridian of the pe- 
riod 1900-05. These constant coordinates are transformed 
to variable coordinates relative to the true pole, equator, 
and prime meridian of the date using the time-varying 
coordinates of the true pole of date relative to the mean 
pole of 1900-05 supplied by the U. S. Naval Observatory. 
This model is relatively simple, should provide a better 
fit to the data, and has the added advantage that solved 
for station locations obtained from missions covering a 
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long period of time should be the same, regardless of 
the wandering of the pole relative to the Earth which 
has occurred during this period of time. 

Section b. gives the formulation for computing body- 
fixed rectangular coordinates of a station where Xt is 
along the intersection of the prime meridian (passing 
through the true axis of rotation) and true equator of 
date; Zb is along the instantaneous axis of rotation di- 
rected north; and t/b completes the right-handed rec- 
tangular system of coordinates. Section c. gives the 
formulation for transforming body-fixed rectangular co- 
ordinates to 1950.0 Earth equatorial position, velocity, 
acceleration, and jerk. The details of the transformations 
are given only for an Earth station. 

b. Body-fixed rectangular coordinates. Fixed stations 
on any body, using constant coordinates relative to true 
pole, equator, and prime meridian: 

A fixed station may be represented by spherical or 
cylindrical coordinates. The spherical coordinates are 

r = radius from center of body to station, km 

<f> = body-centered latitude measured from true equa- 
tor (plane normal to instantaneous axis of rotation 
and containing center of mass) 

A = longitude measured east from prime meridian 
(passing through instantaneous axis of rotation) 

The cylindrical coordinates are: w = (r cos <j>), v = 
(f sin <^), A. 

For spherical coordinates, the body -fixed coordinates Vb 
are 



" Xb" 


“ r cos sin A “ 

rb = 

yb 


r cos </> sin A 




r sin 


(23) 


For cylindrical coordinates, 



■ Xt 


“ M cos A ~ 

n = 

Vb 

~ 

Msin A 


Xb 


V 


(24) 


Fixed Earth station using constant coordinates relative to 
mean pole, equator, and prime meridian of 1900-05: 


The solve for coordinates are spherical or cylindrical 
relative to the mean pole, equator, and prime meridian 
of 1900-05. The spherical coordinates are r, <^o, Ao; the 
cylindrical coordinates are Uo, t?o, Ao. The coordinates 


relative to the true pole are denoted by the same quan- 
tities without the subscript. The transformations are 


<^ — + A<^ (25) 

A = Ao + AA (26) 

M = Mo + Al/ (27) 

t) = Oo + An (28) 


With these quantities, body-centered coordinates r*, are 
computed from Eqs. (23) or (24). 

Fig. 25 shows the latitude <j>o and longitude A<> of a 
station S relative to the mean pole of 1900-05 (Fo), and 
the instantaneous latitude 4> and longitude A relative to 
the true pole of date (F). The pole Po and associated 
grid of equator and meridians are rotated through the 
angle G> carrying Fo to F. The angular coordinates of F 
relative to F« are x measured south along the Greenwich 
meridian of 1900-05 (strictly the 1900-05 meridian of zero 
longitude) and y measured south along the 90 W meridian 
of 1900-05. Values of x and y are supplied weekly by the 
U. S. Naval Observatory. They are represented by linear 
polynomials in the DPODP with coeflBcients specified 
by time block. Since x and y are only a few meters in 
linear dimension, an approximate expression for 
A<^ ~ <f>o is 

A<^ — 3C cos A .0 — 1 / sin Xo (29) 



Fig. 25. Latitude and longitude relative to mean pole 
of 1 900-05 and to true pole of date 
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Noting ao and a on Fig. 24, 


A<, = Uo + tan-* \”^) 

(30) 

A = a + tan-^ ^ ^ 

(31) 

Thus, 


AA *“ A Ao tx Uo 

(32) 

From the spherical triangle P, Po, S, 


sin Oo _ sin a 
cos <\> cos 4>o 

(33) 


Cross multiplying and using Eqs. (25) and (32) gives 

sin Oc cos <l>o = sin (co + AA.) cos {<f>o + A<^) (34) 

Expanding, noting that AA and A</> are very small angles, 
and ignoring the higher order term containing AA 
gives 

AA = tan tto tan <f>o A<j> (35) 

From Eq. (30), 

tan Ao 4- — 

tan tto 

1 — — tan Ao 

X 

Substituting Eqs. (36) and (29) into Eq. (35) gives 

AA = tan <j)o (x sin Ao 4 t/ cos Ao) (37) 

The cylindrical coordinates relative to the pole of 
1900-05 and to the true pole of date are 

Uo = r cos <l>o u = r cos <!> = r cos (<^o 4 A<^) 

(38) 

Vo = r sin<t>o v — r sin <j> ~ r sin (<f>o 4 A</>) 

(39) 

Solving for Am = u “ i/o and Av ~ v — Vo gives 

Am = — i?o A</> (40) 

Ao =: Mo A<t> (41) 

where A<j> is given by Eq. (29). Using cylindrical coor- 
dinates, AA is computed from 

AA == (x sin Ao 4 y cos A«) (42) 


c. Transformation of body-fixed rectangular coordinates 
to 1950.0 position, velocity, acceleration, and jerk. The 
position, velocity, acceleration, and jerk of a station in a 
1950.0 Earth equatorial rectangular coordinate system 
centered at the body on which the station is located are 
denoted by rgo, Fso, Tso, and r'so. The transformations from 
body-fixed coordinates r*, for a fixed station to Tso, Tso, 


Tso and r’so are 

Tso = Tvb (43) 

rso = fvb (44) 

T5o = fn (45) 

r'so = fvb (46) 

where, for the Earth 

Te = {BN AY (47) 

Te - {BN A 4 BN A 4 BNAf (48) 

« {BN AY (49) 

Te - (bn AY (59) 


where the dots indicate differentiation with respect to 
ET. The matrices A, N, and B are 

A = precession matrix, transforming from coordinates 
referred to the mean Earth equator and equinox 
of 1950.0 to coordinates referred to the mean 
Earth equator and equinox of date 

N — nutation matrix, transforming from coordinates 
referred to the mean Earth equator and equinox 
of date to coordinates referred to the true Earth 
equator and equinox of date 

B ~ rotation from coordinates referred to the true 
Earth equator and equinox of date to body-fixed 
coordinates = (X&, t/&, z^Yy where x^, t/&, is a 
right-handed rectangular coordinate system with 
Xft and yb in the Earth's true equatorial plane, x^ 
passing through the prime meridian, and Zb di- 
rected North 

The matrices B, B, B, and *B are computed from 


COS B 

sin B 

0 ■ 



— sin B 

cosB 

0 


(51) 

0 

0 

1 _ 



” — sin ^ 

cos B 

0" 



— cos B 

— sin B 

0 

B 

(52) 

0 

0 

0 
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—COS B 

—sin B 

0" 



T = Julian centuries of 36,525 ephemeris 

B = 


sin B 

— COS^ 

0 


(53) 

since 1900, January 0, 12*^ ET 


- 

0 

0 

0_ 



JED = Julian ephemeris date 


“ 

sin B 

— COS^ 

0" 



_ JED - 241,5020 

B = 


COS B 

sin B 

0 


(54) 

36,525 



0 

0 

0 





(61) 


where 

B = apparent (true) sidereal time = Greenwich hour 
angle of true equinox of date 


The parameters used in Eqs. (55) through (61) are de- 
fined as 

Bu = mean sidereal time = Greenwich hour angle of 
mean equinox of date 


The quantities of B and B are 

computed from the fol- 

= nutation in longitude 


lowing formulation 



c = true obliquity of ecliptic 


B — Bm + cos c 
B = Bm cos € 

(rad) 

— €8^ sin € 

(55) 

T — mean obliquity of ecliptic 
8c = nutation in obliquity 


e = Se 

(rad/ephem sec) 
(rad) 

(56) 

(57) 

The mean sidereal time is computed as a function of UTl 
(See SPS 37-39, Vol. Ill, p. 36) 

Om = UTl + / + KTv + LTl (62) 

€ = €■+ S€ 
(Ref. 26, p. 98) 

(rad/ephem sec) 

(58) 

where 

J = 6''38"'45f836 = 23,925f836 


__ A + BT + CP + DT» , j, 
' 206,264.80625 

(59) 

K = 8,640, 184f542 
L = 0?0929 



A = 23^27'8'/26 = 84,428'/26 
B = -46'/845 
C = -(K'0059 
D = (K'OOISI 

B + 2CT + 3DT2 


Tu = number of Julian centuries of 36,525 days of 
UTl elapsed since 1900, January 0,12*" UTl 


€ — 


^ _ JD (UTl) - 241,5020 
36,525 

JD (UTl) = Julian date computed from UTl 


(63) 


86,400 X 36,525 X 206,264.80625 

(rad/ephem sec) (60) 


The true sidereal time B may be computed from Eqs. (55) 
and (62) as 


_ f/rTri -i_ ^ , H COS C 


= (UTI + - 


86,400 


+ 


27T 


-) V 


(rad) (64) 


where 

UTi = [JD (UTl) - 0.5UcUn.i (65) 


DifFerentiating Eq. (62) with respect to ET gives 


$u — 


dUTl / K + 2I.Tv \ 7T 

dET \ 36,525 X 86,400 ) 43,200 


(rad/sec) 


( 66 ) 


From SPS 37-39, Vol. Ill, p. 37, 

UTl = ET - (ET - A.1) - (A.1 - UTl) 


and 


_^UT^ =14. cesium 

dET 9,192,631,770 


- g - 2ht 


(67) 

(68) 
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where 


t = seconds past start of time block containing f, g, and h 


Thus, 


=( 


— ( 1 + 


K + 2LTu 


y){ 


36,525 X 86: 
and 0 is computed from Eqs. (56) and (69). 


1 + 




9,192,631,770 


(^) 


( 69 ) 
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F. Ranged Angle Corrections 
Due to the Ionosphere 

A. Liu 

1. Introduction 

Radio waves passing through an ionosphere are not 
only deflected from a straight line, (refracted), but also 
the speed at which the wave is propagated is altered. The 


refraction due to the ionosphere for S-band frequencies 
can be as large as 12 sec of arc, although a more typical 
value is 6 sec of arc. The range correction can be as large 
as 20 m at low elevation angles. The ionospheric effects 
turn out to be more severe for missions to the inner 
planets than for missions to the outer planets. The reason 
is that for missions to Venus, for example, the observa- 
tions are always taken more nearly in the direction of 
the Sun, and hence, through a much thicker ionosphere 
than for observations taken for a Mars mission, where for 
about half the duration of the mission, the measurements 
are taken at night (away from the Sun). Fig. 26 shows 
the ion density as a function of direction of Sun and 
altitude. 


5 AM 


3 AM 



II AM 


A'# < 10,000 



II PM 


7 PM 


I PM 3 PM 5 PM 

Fig. 26. Ion density vs Sun direction and altitude 


A program was written to ray-trace through a Chap- 
man ionosphere, and curves are presented to illustrate 
the effects of variations in the ionospheric model upon 
the range and angular measurements as a function of ele- 
vation angle. 


2. Ray Tracing Formulae 

(SPS 37^28, Vol. IV, pp. 1 and 2) 

The assumptions made are 

(1) The ionosphere is spherically layered, the density 
being given by a Chapman distribution. 

(2) The ray path takes the shortest time route between 
two given points (Fermat’s Principle). 

Define the following 

c = speed of light in vacuum 
V = speed of light in medium 
(f, = polar angle 
n = index of refraction = c/v 

The problem is to find the curve C for which the travel 
time is at a minimum. 
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Consider the situation illustrated in Fig. 27. A light ray 
is required to leave point A, pass through a medium 
with index of refraction n, and arrive at point B in the 
least amount of time, T, i.e., 


= f iii 

Jc r 


is to be an extremal by finding the 
proper path C, or (1) 


where 


Cl = To rio cos Eo 


(8) 


a. Angular refraction. To evaluate the integral given 
by Eq. (7) make a change of variable to w = m. It fol- 
lows that, 


du = (m' -h n)dr 


where 


n' = 


dn 

dr 


(9) 


( 10 ) 


Integration by parts of Eq. (7) gives 



- (12) 




where n is the index of refraction. 


The path C is determined by finding = </>(r), A 
necessary condition for rendering T an extremal is that 
the integral of Eq. (2) satisfy the Euler-Lagrange equa- 
tion. 


where 


and 



F = n(<^, r) r, </>i) (5) 


We note that will have a singularity whenever u = Ci, 
and to remove this, singularity transform again the var- 
iable w to X by 

cosX = -^, (13) 

yielding for 8<^, 

** = r 

J X ~ 

b. '"Range” or time of flight correction. The path is 
now known, given by Eq. (14) and the transit time, T 
or R, the measured range, is given by 



At this point, we note that for the case where spherical 
symmetry exists. 


dF 

d(f> 


0 


The total bending to infinity can be found as 

dr 


s= r 

(rVn='r^-q) 


Changing variable r to u and integrating Eq. (15) by 
parts as before gives 


(6) AR = r ~ ^ \ du 


(7) 


+ foTio sin Eo 
= SR — fo sin (Eo 


To sin (Eo — 8<t>) 
(16) 


+ fo n„ sin E„ (no = 1) 
(17) 
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We note again, that 8H suffers the same singularities 
as for u — C^ and to remove this singularity, a trans- 
formation to X as before will give 

The actual numerical method for integrating Eqs. (14) 
and (18) depends upon the function in the integrand, 
which in turn depends upon the model atmosphere. The 
model atmosphere chosen is given by 

"Ne = 2V(max)exp [1 ~ u — exp( — m)]}, (19) 

where 

"Ne = number of electrons/cm^ at height h 

N(max) — maximum number of electrons/cm^ occurring 
at h(max) 

h = altitude above Earth 
B =: scale height and, 
u — [h — h(max)]/B 

The index of refraction is related to Ne by ; 

n = 1 -^yNe, (20) 

where 

f = transmitter frequency = 2.3 X 10® Hz 
— 80.6 for Ne in units of electrons/m"^ 

3 . Numerical Results 

Eqs. (14) and (18), which give the angular and range 
refraction corrections^, respectively, for a given Chapman 
ionosphere, have yielded numerical estimates for various 
assumed model atmospheres as a function of elevation 
angle. These estimates are summarized in Figs. 28-33. 

Fig. 28 shows the variations in the range correction due 
to fluctuations of one order of magnitude in the maximum 
electron density. These corrections vary directly and 
linearly with the electron density and are the greatest 
variations in the corrections. The tropospheric range 
correction is superimposed for comparison purposes. 

Fig. 29 shows the variations in the range correction due 
to changes in scale height, i.e., variations in temperature. 

“The range corrections are in meters, and angular corrections are 
given in microradians ( l/xrad ^ 0.2 sec) 



ELEVATION, deg 


Fig. 28. Range correction as a function 
of electron density 



ELEVATION ANGLE, deg 

Fig. 29. Range correction as a function of 
scale height 



ELEVATION ANGLE, deg 

Fig. 30. Range correction as a function of HMAX 
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Fig. 31. Elevation angle correction os a 
function of electron density 


These changes produce the next greatest variations in the 
range corrections. 

Fig. 30 shows the changes in the range correction due 
to variations in the height of maximum electron density 
(HMAX). HMAX will vary from 200 to 400 km during a 
day, and the curves illustrate the fact that the range 
corrections are relatively insensitive to changes in HMAX. 

Figs. 31-33 show the angular corrections for the same 
set of parametric fluctuations. They may be summarized 
by saying that for the most part the corrections are on the 
order of 5" of arc. The largest correction of this set is 
12" of arc (Fig. 33). From these curves, it appears that 
angular corrections are not important enough to warrant 
further study, other than the manner in which they inter- 
act with the range corrections. 



Fig. 32. Elevation angle correction as a 
function of scale height 
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Fig. 33. Elevation angle correction as a function of HMAX 
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G. Apollo Project Support 

D. W. Curkendall 

The task reported here has as its objectives to provide 
the systems analysis effort to support the development of 
Unified S-band to be used in the Apollo program, to 
develop the capability of participating in the ApoUo 
S-band qualification tests, and to assist the Manned 
Spaceflight Center (MSC) and Goddard Spaceflight Cen- 
ter (GSFC) in the planning and execution of the naviga- 
tion aspects of the Apollo lunar landing mission. Specific 
items under this task will be selected each report period 
for further discussion. 

The Apollo task includes, but only partially funds, items 
reported under Tracking and Navigational Accuracy 
found elsewhere in Sect. II of this volume. In Sects. H 
and I, D. Curkendall contributes two articles concerning 
the doppler data noise model and its relation to orbit 
accuracy. In the first of these, the influence of the oscil- 
lator stability on the data is analyzed assuming that the 
oscillator generates white frequency noise. Both two- and 
three-way doppler are investigated and it is shown that 
while the oscillator frequency noise is impressed directly 
onto the three-way data, a transformation to colored 
phase noise takes place when considering the two-way 
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configuration. The importance of this difference is 
discussed. 

In the second article, an uncorrelated phase noise dop- 
pler model is applied to the three-dimensional navigation 
problem first analyzed by Hamilton and Melbourne using 
the white frequency noise model in SPS 37-39, Vol. Ill, 
pp. 18-23. The results obtained, which indicate a signifi- 
cant increase in accuracy over those obtained by Hamil- 
ton and Melbourne, are related to the oscillator error 
source equations and it is shown that the information 
content in a single pass of doppler data is inversely pro- 
portional to the round-trip light time to the probe when 
oscillator noise is the dominant error. 

Also found in Sect. J is a paper entitled ‘‘Orbit Deter- 
mination of Two Spacecraft Through the Use of Multiple 
Link Tracking Data.” As the title implies, the article is 
a study of the navigational accuracy resulting when 
multiple-link doppler is obtained by coherently trans- 
mitting from the Earth to the first spacecraft, to the 
second, and finally back to the Earth. Three examples 
of performance are discussed. 

The first example analyzes a direct LEM ascent to a 
rendezvous with the CSM. It is found that the addition 
of the multiple link data results in an orbital accuracy 
improvement of a factor of five over that using two- and 
three-way doppler alone. This is thought to be a marginal 
improvement considering the difficulty involved. How- 
ever, if the three-way data is deleted, the improvement 
factor is approximately 50, suggesting that multiple-link 
data can supplant the function of three-way data for a 
rendezvous at interplanetary distances. 

The second example analyzes the case of tracking a 
spacecraft which is behind the Moon using a lunar com- 
munications satellite as a coherent relay for the tracking 
signal. It is seen that entirely adequate orbital estimates 
can be obtained before direct contact is made with 
the spacecraft. 

The final case examines the tracking of a released 
capsule-lander and its parent spacecraft for an inter- 
planetary mission. It is found that the multiple link data 
cannot establish the orbits of the two spacecraft much 
better than two-way data alone at least in the time span 
9 days to 1 day from encounter when the probes are 
experiencing nearly rectilinear motion. The relative orbits, 
however, are found to be precisely determined if the 
interspacecraft link is added. 
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H. The Influence^of Oscillator 
Instability on Orbit Accuracy 

D. W. Curkendall 

7 . Introduction 

In previous issues of this Summary, the ability of 
doppler tracking data corrupted by stationary-phase 
(bounded-range) errors to determine the radial velocity 
of a spacecraft was investigated. In SPS 37-33, Vol. Ill, 
pp. 20-24, the assumption was made that the noise con- 
sisted of white phase jitter; in SPS 37-39, Vol. Ill, 
pp. 30-34, this assumption was relaxed and the case of 
exponentially correlated noise was treated. Many of the 
error sources present in a doppler tracking system con- 
form very well to this simple model and the results 
derived can be very useful in assessing the importance 
of such errors as counter truncation (SPS 37-38, Vol. Ill, 
pp. 8-13), atmospheric scintillation, and RF loop phase 
jitter. An error source of particular importance and whose 
noise characteristics contain significant departures from 
the simple models treated before is the frequency drift 
inherent in the master oscillator in the tracking station. 

This article derives the statistical properties of the 
noise impressed on the data assuming that the frequency 
spectrum of the oscillator is flat. Both two-way (coherent) 
and three-way (non-coherent) tracking modes are con- 
sidered. As before, application of the resulting noise 
model will be made to a simple one-dimensional tracking 
problem where it is desired to determine the velocity of 
a probe receding radially from the tracking station. 

2 . Noise Model for Two-Way (Coherent) Doppler 

Figure 34(a) depicts a convenient visualization of a 
two-way doppler tracking system. The frequency trans- 
mitted is the sum from two oscillators; one whose output 
is a perfectly stable fo, and a second noise generator 
whose output is time-dependent, or /e(t). For this analy- 
sis, it is assumed that the spectrum of fe is flat; i.e., 

E S(u - v) ( 1 ) 

where E [ ] denotes the expected value of [ ] and 8 is 
the Dirac delta function. This sum frequency is trans- 
mitted to the spacecraft and returns r seconds later, 
translated an amount proportional to the range-rate, p, 
of the spacecraft. Upon return, the signal is beat with 
the current transmitter frequency as shown. DiflFerence 
frequencies are then integrated producing range change 
data, Ap, and noise. 
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Suppose we are taking a counter reading every M sec- 
onds and treating the sequence of readings so obtained 
as the data vector (this is the total-count data-taking 
scheme discussed in SFS 37-38, VoL III, pp. 20-24). 
Then, from Fig. 34(b) the error in the ith measurement is: 

= 

c, (2) 

r%At fiM-T 

^ I - j fe{u)du 

then with the restriction that iAt > t 



Note that the left-hand term is not dependent on t and 
will be present in all measurements on a single pass of 
data. It is thus possible to treat this term as a constant, 
subject to estimation; thus, we will re-define the tth error: 



Fig. 35. Normalized autocorrelation function for 
two-way coherent doppler 


3. Noise Model for Three-Way (Non-Coherent) 
Doppler 

In the three-way configuration, the essential change 
from Fig. 34(a) is that a second oscillator is used to dif- 
ference the incoming doppler signal. Thus, assuming 
that there is no unknown bias between the two oscilla- 
tors, the error in the ith measurement is written: 


(, = C,~K= fe(t)dt (4) 

JiAt-r 

We may also compute the auto-covariance between 
the ith and the (i + /)th errors as: 

( fiAt f(i+J)At ) 

E Uiiui] = EJ fe (u)du / fe {v)dv [ (5) 

iJiAt-T J{i+j)At-T j 

Interchanging the ordering of integration and the ex- 
pectation operation and employing Eq. (1) 

= <t|t 1^1 for I / 1 < T 

E 

= 0 for I / 1 At > T 

In particular when i = 0 

E[q] -cjr (7) 

Fig, 35 illustrates the normalized auto-correction function 
obtained from Eq, (6), and assuming / a continuous 
variable. Note that Eq. (7) shows the variance of the 
error to be stationary and not an increasing function 
with i. 


/ iAt 

K(t)-U^{t-r)dt (8) 

where the 1 and 2 subscripts refer to the independent 
oscillators. Then: 

E [Ci = 

( fiAt r{i+})At ) 

E jj J (fiju) - f,Ju)) (U^iv) - U^(v)) du dv j 

= 2aHM i > 0 (9) 

= 2al(i + i)M i < 0 

It is readily seen that the variance of the error grows 
linearly with time as contrasted with the stationary error 
found to be present in the previous case. For subsequent 
(6) analysis, it will prove convenient to also obtain the statis- 
tics assuming that the data will be processed after differ- 
encing each counter reading with the previous readout 
and dividing the result by At, i.e.: 

^ Api Api~\ ~h Cj (10) 

where the subscript i refers to the ith measurement. 

The statistics on the error 

€i = Ci- Ci-jAt 
are easily computed. 
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1 

A#" 


E 





- fi(v) f2(u-r) - fi^iu) fi^iv-r) 
+ hjs>—r) dudv > 

r 2<r? 

ifi = i 

= 0 if i ^ / 


( 11 ) 


The reader will now readily see the motivation for ob- 
taining two equivalent descriptions of the three-way 
noise characteristics. By reverting to the diflFerenced data 
type, we have obtained an uncorrelated model which will 
prove of distinct convenience when application is made 
to the tracking problem in Sect. 5. 


4 . Application of the Two-Way Model to a 
One-Dimensional Problem 

Imagine that we are tracking a spacecraft moving radi- 
ally from the station with constant but unknown speed v. 
The probe is tracked for T seconds, and the integrated 
doppler is read-out every At seconds for a total of N 
measurements (NAt = T). For two-way coherent track- 
ing, the errors, in the measurements are described in 
Eq. (4). Relating the measurements, t/i, to the unknown 
velocity v and the constant K defined again by Eq. (4), 
we have: 


yi 


O = A 


V 

K 




( 12 ) 


where A is the partial derivative matrix of the y’s with 
respect to v and K, i.e., 


A^ = 


1 , 


2Af, 

1 , 


If At > T, then 


NAt 

1 


E = af rl = Ae 

where I is the N X N identity 


(13) 


(14) 


If the data is then used to determine v and K and a min- 
imum variance (MV) estimate is formed, the covariance 


(2X2) of the errors in the estimate is given by 


(15) 


This description of the problem is identical to that 
presented in SPS 37-38, Vol. Ill, pp. 20-24, which assumed 
that the noise on the data consisted of white phase errors 
in the read-out of the range-change. The one difference 
is that the phase error having variance is replaced 
here by the term <tJt, We may quote the results of this 
reference directly then, and it is found that the variance 
of the error in the estimate of v (we are not directly inter- 
ested in how well we may determine K) is given as 

<?(N) = (16) 

which for large N 


Q(N) = 


12aj T 

T^N 


These results do not include the a-priori knowledge 
on K (which is easily computed as <r^r); if this is in- 
cluded, we have: 

— 12o-f tN 

Q(N) T‘(N + 1) (N+2) 

which has the same asymptotic behavior as Eq. (14). 

These results underscore the first important conclusion: 
If the interval between samples is greater than the round- 
trip light time, the error variance in the estimate goes 
inversely with the number of samples. Moreover, the 
sample errors are uncorrelated producing a diagonal noise 
matrix (Eq. 14) and the MV estimate is easily formed. 

Increasing the sampling density. If, in the above, 
Af < T it is no longer true that the noise co-variance 
matrix is uncorrelated and we have no previous analysis 
to quote in order to specify the error variance. In order 
to see the effect of sampling at a higher rate than once 
per round-trip light time, it is convenient to pick a single 
interval of t and vary the number of samples taken in 


TRANSMIT 

RECEIVE 
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^’"1 I I ^ 
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I /*o 
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Fig. 36. Observation schedule for high-density sampling 
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that one-interval. Fig. 36 depicts the observation sched- 
ule which calls for measurements at 


' T, T 


n’ 


T 


n ’ 


r (n-l)r 
n 


or a total of n measurements within this single t interval. 
Here the small n is to be differentiated from the total 
of N samples for the entire data span. 


Formulating the problem in this matter, the correlation 
between and is easily computed 





(18) 


Using (18), the n X n covariance matrix of observation 
noise is given by: 


A.= 




n 

\ 

\ n 
\ 

Symmetri^ 


n — 1 n — 2 
n — 1 
n 


1 

2 

3 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ n — 1 
n 


(19) 


•• r ^ ~ (/~o ~ 

u 


In order to form a MV estimate with these data, we have 
need of the inverse of this matrix. It is easily seen to be: 


(n-l-2) 
(n -1-1)2 

_ 1 _ 

2 


0 

0 


0 

1 

2(n-hl) 


X 0 0 ■ • ■ 0 2(„+i) 

1 -X 0 " 


2 \ \ 2 \ 

\ \ \ 

\ \ \ 
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\ 
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\ 
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\ \ \ 

\ \ 

\ \ 

\ \ 

\ \j_ 

\j_ (n+2) 


\ 


\ 


2 2(n-l-l) 


( 20 ) 


This expression is valid for n > 3, the inverse for the >Theauthor wishes to thank C. B. SoUoway for obtaining the inverse 

2X2 and 1X1 are of course easily computed directly. of Eq. ( 19 ). 


! 
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For simplicity, let us now assume that the constant K 
has been sufficiently well-determined (from the rest of 
the data, perhaps) so that we may drop it from the esti- 
mate list and form a one-dimensional estimate of v. The 
A matrix in Eq. (14) is now given by: 

^ _ (n-2)r y _ (n-lM 

n ’ ” J 

Inserting these expressions into (13) and performing the 
indicated summations, the final result becomes 




= f7’ T —— T — — 

L ’ n’ 




for n = 1 


= Q(n) = 


9<Ty T 


4 ( 37’2 _ 


f or n = 2 


(n+1) 

n [2T^ - 2Tt + T^] + 2Tt - 


f or n ^ 3 

( 21 ) 


For T > > T, the general form of (21) reduces to 

<r| T(n+ 1) 


Q{n) = 


2T^n 


for r > > r (22) 


and for large n 


Q(n) ^ for T > > T and large n 


(23) 


Comparing (23) and (21a), we see that as n-^oo, the 
error variance reduces by only a factor of 2 over that 
obtained for n = 1. Thus, we reach the 2nd important 
conclusion: Sampling faster than the round-trip light time 
brings but marginal improvement to the orbit estimate. 
Moreover, the noise covariance matrix is no longer 
diagonal and the MV estimate is difficult to form. 


5. Application of the Three-Way Model 

In applying the three-way noise model to the same one- 
dimensional problem, it is convenient to assume the data 
type used is the differenced doppler (11) whose noise 
model is specified by (12). Thus, 

y = Av -h £ 


and 


E 




I 


Application of (15) results in the mean-squared un- 
certainty in v: 


“ T' 


(24) 


Eq. (24) is to be contrasted with (17) which shows the 
variance goes as r/T^N. Thus it is seen that the lack of 
coherence in the three-way data dramatically increases 
the effects of the oscillator instability. 


6. Concluding Remarks 

The characteristics of the noise impressed on the 
doppler data have been obtained assuming a flat fre- 
quency spectrum present in the station s master oscillator. 
The important results of this analysis have been under- 
scored in the main body of text and will not be repeated 
here. Although the models used are simple ones, it is 
argued that these conclusions apply to more complex 
circumstances. Found elsewhere in this volume is a 
second article contributed by the author which applies 
the bounded-range error model to a three-dimensional 
tracking problem first analyzed by Hamilton and Mel- 
bourne in SPS 37-39, Vol. Ill, pp. 18-23. 


fi(>7 U430 

I. Application of a Bounded Range 
Error Model to the Evaluation of 
the Information Content of a 
Single Pass of Doppler Data 

D. W. Curkendall 


7. Introduction 

In SPS 37-39, Vol. Ill, pp. 18-23, Hamilton and Mel- 
bourne presented an analysis determining the information 
content of a single pass of doppler data from a distant 
spacecraft. It was shown that the range rate of the probe 
with respect to the observing tracking station could be 
approximated by: 


where = [1, 1, 1 1] 


p = f + f s<t> cos 8 sin<ot + <i)tors<o cos 8 cos <at (1) 
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where 

r = average geocentric range-rate during the pass 

fs = distance of the station from the Earth s spin axis 

<0 = Earth’s rotation rate 

8 = declination of the probe 

t = time measured from the expected time that the 
probe will cross the meridian of the station 

to = actual time of meridian crossing. 

By making successive observations of p and performing 
a regression analysis using Eq. (1), estimates of f, cos 8, 
and to (which is related to the probe’s right ascension) are 
obtained. The original analysis assumed that the noise 
in the successive observations of p was white, gaussian, 
and of zero mean. Using this assumption, results con- 
cerning the accuracy obtainable in the estimated param- 
eters were obtained without explicit reference to the 
system which secured the observations. 

It is our purpose here to make an explicit focus on the 
measurement system used and perform a similar analysis 
employing the noise model used in SPS 37-38, Vol. Ill, 
pp. 20-24. It was argued in this reference that when 
tracking a spacecraft in the two-way coherent mode, the 
primary error sources are better characterized as white 
range (as opposed to range rate) errors. This argument is 
pursued further in a separate article in this Summary, and 
here it is shown that the frequency drift inherent in the 
stations master oscillator produces doppler errors con- 
forming to this model if the time between successive 
samples is greater than the round-trip light time to the 
probe. 

As in the original article, calculations of estimation 
accuracy are made and it is shown that the white range 
error model produces significantly better estimates of 
perfonnance. 

2 . Regression Analysis 

To facilitate the application of the bounded range 
error model, it is convenient to use the total counted 
doppler as the basic data type. That is, the doppler fre- 
quency is integrated by means of an electronic counter, 
and its current reading minus the reading at the start of 
the tracking period is the observable. This quantity is 
interpreted as the change in the station-probe range. 


Thus using the definitions: 
a ~ r 

b = r^ cos 8 
c = (iitorg cos 8 
and integrating Eq. (1) 

Ap(fi) = p{t) dt = a{ti — ts) — b{cos iiiti — cos <o 

+ c(sin <D t\ — sin <o (2) 

where ts is the start of pass, and ti the time of read-out 
at the ith data time. The ith data point then will be the 
range change plus noise: 

di = Ap(fi) -f Ui — Ho (3) 

where the no is an error common to all measurements and 
is interpreted as the initial phase error at the start of the 
pass (or alternately, as a combination of the initial phase 
error and the oscillator induced error common to all 
measurement, that was discussed in the preceding article, 
Sect. H). 



Fig. 37. RMS estimation accuracy, (Tc versos 

half-width of pass, \j/, for crAp = 0.06 m, 
1*min samples 
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Our objective is to determine the accuracy to which 
the constants a, b, and c may be estimated if Ui is as- 
sumed uncorrelated with n, and has stationary variance, 
. If a minimum variance estimate is performed, the 
covariance of uncertainty in the estimates is written 

A=(A^A;M)-‘ (4) 

where A is the partial derivative of the data with respect 
to t parameters to be estimated, and A« is the covariance 
matrix of the noise on the data. 


The Ac matrix is easily demonstrated to be diagonal if 
Ho is treated as a parameter subject to estimation, thus 
making (4) a 4 X 4 matrix equation. Although we have 
no inherent interest in obtaining the value of n©, formu- 
lating the problem in this manner greatly simplifies the 
computations: 

We have 

A-» = A’-A-* A =-^A^A = 


8 


24^ 


, {if/ sin tj/ - cos ^), — [(sin ^(2 + - if/ cos 


Cl> CO 


\j/ + ^ sin ^ cos \f> + 2»i/ cos^ <j/, —2 [sin^ ^ ^ sin f cos V']> ~ 2 [sin ^ cos] 


N 

24 <TAp^ 


^ ^ ^ ^ ^ 2^ sin ^ 


Symmetric 


24 


(5) 


where ^1/ is the half -width of the pass (in degrees of Earth 
rotation), and N is the number of measurements taken. 
Eq. (5) was written assuming that N was large and that 
the indicated summations in A^A could be replaced by 
integrals. This equation is not conveniently inverted 
analytically; the dependence on N and crAp can, however, 
be isolated from and Eq. (4) can be re-written as: 

A = %^^PW (6) 


cTAp was taken as 0.06 m (this is equivalent to the 0.001 
m/sec figure used in SFS 37-38, Vol. III). Fig. 37 is to be 
compared to Fig. 36 in the original article. Sect. H. 
Note that the change to a bounded range error model 
results in an estimated accuracy improvement of approxi- 
mately a factor of 50 for an entire pass of data (^=90 deg). 

As shown, for ^ = 90 deg, the following accuracies are 
obtained: 


= 0.020 m 


Eq. (6) was solved numerically, and the results are = 0.007 m 

plotted in Fig. 37. For this computation, it was assumed 

that data was obtained at the rate of 1 sample/min, and <tc = 0.026 m 
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3 . Inferprefafion of Results 

The numerical results quoted assume a sample rate of 
1/min. Eq. (6) indicates that sampling at higher rates 
will result in even increased accuracy estimates by a 
factor of because <7Ap is assumed not to be a 

function of the sample spacing. As an aside, the original 
article derived a similar error variance equation pro- 
portional to <7.“ /N, but left concealed by this equation 
alone is the factor that is proportional to the sample 
rate when using the assumption of uncorrelated range 
rate noise. Thus the results quoted in the SPS 37-38, Vol. 
Ill, are insensitive to N as long as it is large. 

The formulation here, however, gives rise to the un- 
realistic result that we may drive the error variance to 
zero by simply taking samples at a fast enough rate. 
In order to keep the model realistic, let us assume that 
the data noise is entirely produced by the high-frequency 
drift in the oscillator and draw upon the result of the 
previous article (Sect. H) which shows that the data are 
uncorrelated only if the sample time is greater than the 
round-trip light time to the probe. Thus we must restrict 
N to be 

Nmax = ^ (7) 

0>T 

where r is the round trip light time. In addition, it has 
also been shown that the variance of the incremental 
range error is: 

o-Ap^ — ajf T (8) 

where a/ is a measure of the oscillators performance. 
Substituting (7) and (8) into (6) we obtain 

A = m (9) 

Thus, we see that the performance degrades linearly 
with T, and Fig. 37 is reinterpreted to be the results 
obtainable for a 1-min light time and an oscillator whose 
rms frequency departure A/// is approximately 5 X 10"'^ 
for a 1-sec averaging time. As an example, for a probe 
distance of 1 AU, the numbers given would be increased 
by a factor of 17. 

The use of the estimates of a, h, and c to obtain f, 8, 
and Tg, and the relationship between these estimates and 
errors in station longitude and Universal Time were 
discussed in the reference and will not be repeated here. 

It should be cautioned, however, that the estimation 
procedure implied by this analysis is very sensitive to 


modeling errors or omissions. We have assumed that we 
can measure the range change to 0.06 m, a figure that is 
below our current knowledge of ionospheric eflFects and 
is probably below the perturbations on the spacecraft due 
to random non-gravitational forces. This assumption is to 
be contrasted with that made in the reference which 
implied that at any time during the pass the velocity 
could be measured to 0.001 m/sec accuracy. 

As a final comment, it is noted that in order to achieve 
the accuracies quoted, continuous RF lock must be main- 
tained with the spacecraft throughout the pass. In prac- 
tice, each time lock is lost and then regained, a new rio is 
created and its estimation must be re-initiated. This (even 
assuming the data loss to be negligible) will degrade the 
estimation performance. In the limit, if lock is lost 
between each sample, the situation will revert to that 
characterized by independent range-rate noise, and the 
analysis of the reference applies. 


N67 14431 

J. Orbit Determination of Two 
Spacecraft Through the Use of 
Earth-Based Multiple-Link 
Doppler Tracking Data 

D. W. Curkendatl 

1. Introduction 

The high accuracy of Earth-based doppler tracking 
data and its ability to produce precision orbits for a wide 
spectrum of lunar and planetary missions is well known 
and has recently been discussed in several SPS articles. 
In certain future missions, involving more than one 
spacecraft (or a spacecraft and a capsule), it may be 
useful to extend the basic concept of two-way doppler, 
i.e., tracking station to spacecraft to tracking station, to 
include an additional spacecraft-to-spacecraft link. With 
this scheme, provisions are made to transmit the received 
signal from the first vehicle to the second before its 
return to Earth. Phase coherence with the original trans- 
mission is maintained in all three links and, upon its 
return to the tracking station, the signal would contain 
a doppler shift proportional to the algebraic sum of the 
three relative range rates involved: 

d = P«1 + Pl2 + P2S (1) 
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where s refers to the station, and the numbers identify 
the two spacecraft. 

The qualitative motivation for such a data type is that 
orbital elements which are unobservable from Earth, or 
that are observed so obliquely as to be poorly deter- 
mined, may come solidly into view. Of course, there is 
nothing new in the notion of spacecraft-to-spacecraft 
doppler links; the performance of these measurements 
from an Earth-based system, however, would seem to 
offer the important benefits of increased accuracy with a 
minimum of complexity of the on-board equipment. 
This article explores the utility of obtaining such a data 
type by constructing three navigation problems, each 
involving two spacecraft, and calculating the resultant 
navigation accuracy using the inter-spacecraft link as a 
data type.^ 

2 . Numerical Results: Three Examples 

The numerical results presented in this section were 
obtained from a computer study program specifically 
constructed to analyze the ability to determine the orbits 
of two vehicles simultaneously when data is available 
which relates to the state of both vehicles.^ 

In this program, the capability to obtain five diflFerent 
data types at each of a series of data times was mecha- 
nized; two of these data types have an associated bias. 
The five data types are 

(1) Two-way doppler from Spacecraft No. 1. 

(2) Two-way doppler from Spacecraft No. 2. 

(3) Three-way doppler from Spacecraft No. 1 (this data 
type has a bias associated with it). 

(4) Three-way doppler from Spacecraft No. 2 (with a 
bias). 

(5) Multiplier-link (or stereo) doppler from the Station 
to Spacecraft No. 1, to Spacecraft No. 2, and back 
to the Station. 

All of the physical constants, including the station loca- 
tions, were considered known. Because of this, these 
numerical results may not be taken too literally; the 
comparisons, however, of the relative power of the dif- 
ferent data configurations should be valid. 

'For a description of the mathematical development involved, the 
reader is referred to A AS preprint 66-110, a paper having the same 
title as this article which is to be published in Advances in the 
Astronautical Sciences, 


3. A Rendezvous Problem, Example 1 

The first example treats a rendezvous situation and, for 
illustration, the trajectories used were taken from an 
Apollo mission profile of direct Lunar Excursion Module 
(LEM) ascent to the Command and Service Module 
(CSM). Fig. 38 shows the basic geometry and tracking 
pattern. Ref. 27 was the source for and describes in detail 
the trajectory used, but briefly the LEM rises from 
approximately zero degrees latitude and longitude and 
injects into a slightly elliptical orbit of 1820-km semi- 
major axis. The CSM is in a well-established circular 
lunar orbit of 1886-km semi-major axis. The orbits have 
nearly zero degrees inclination relative to the lunar 
equator and are inclined with respect to one another by 
0,5 deg. From LEM injection, the two vehicles are co- 
visible to one another and the Earth for about 30 min 
before LEM occultation. Rendezvous occurs some 52 min 



EARTH 


Fig. 38. Basic geometry and tracking pattern, including 
inter-spacecraft doppler links for a 
rendezvous situation 
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after LEM injection as a point in space not visible from 
the Earth. For simplicity, it was assumed that no mid- 
course maneuvers, interrupting the tracking and corrupt- 
ing the orbit, were performed. 

For the basic case, data types 1 to 5 were used, and 
then a second case for comparison was run omitting the 
stereo data. The two-way tracking originates from the 
Madrid station, and the Ascension Island station was 
used to obtain three-way data. These were picked to 
eflFect a large latitude difference between the stations, 
which is the direction found by Ref. 27 to be the more 
critical for tracking low-inchnation lunar orbits. The 
noise on the doppler was assumed to be 0.03 m/s for 
1-min samples, a highly conservative number based on 
present experience but used to offset the neglect in the 
uncertainties in the physical constants. The three-way 
data contain a bias arising from an unknown offset in the 
frequency of the oscillators in the master and slave track- 
ing stations. This was taken as 0.03 m/s, representing 
an unknown offset of two parts in 10"^®. The a-priori 
knowledge of the CSM orbit was assumed to be 1 km in 
position and 0.5 m/s in each of the three directions; the 
LEM s initial uncertainty was assumed to be 10 km in 
position and 1 km/s in velocity. 

Fig. 39 plots the results from both of these cases. This 
figure shows the rms uncertainties in the prediction of 
the difference of the rendezvous conditions, both position 
and velocity, plotted as a function of tracking time. 
Although the relative performance between the two data 
configurations varies considerably over the time span 
considered, adding the stereo linkage improves the track- 
ing accuracy significantly, if not spectacularly. After 29 
min of tracking, the rms error in predicting both the 
final position and velocity is approximately reduced by 
a factor of 5. 

The prediction accuracy of the LEM s absolute final 
position and velocity is not shown here, but after the 
full 29 min, the figures are 0.7 km and 0.6 m/s using the 
stereo linkage, and 1.5 km and 1.5 m/s without. Thus, 
a good portion of the factor of 5 shown on the graphs is 
realized, not by better absolute prediction capability, but 
by the ability of the spacecraft-to-spacecraft link to relate 
the two orbits in question to each other. 

A second factor investigated was the effect of removing 
three-way doppler as a data-type from each of the two 
data configurations. The importance of the three-way 
doppler is that, similar to the stereo linkage, it affords 



the observer a second line-of -sight with which to view 
the orbit. The two tracking stations considered are dis- 
placed by 4500 km, which when viewing objects at lunar 
distance, results in pointing directions nearly 0.7 deg 
apart. The effect on the removal of this second line-of- 
sight is dramatically shown in Fig. 40, which plots the 
same information as before, except that no three-way 
doppler data are taken. The accuracy produced by two- 
way doppler alone increases to a final value of 21 km 
and 56 m/s. By way of contrast, the accuracy produced 
by the data configuration which includes the stereo link- 
age degrades but slightly to 0.5 km and 0.4 m/s. The 
implication is clear: at lunar distances, three-way doppler 
data are nearly as powerful as the more difficultly ob- 
tained multiple-link data and the performance gain is 
little more than nominal. At interplanetary distances, 
however, little parallax can be achieved by tracking from 
two Earth'bound stations; in this circumstance, obtaining 
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T, min 

Fig. 40. RMS rendezvous position and velocity errors 
without three-way doppler 

directly or indirectly inter-spacecraft range-rate data can 
easily supplant the function served by the three-way data. 

4. Tracking an Occulted Spacecraft, Example 2 

In our second example we shall use the same multiple- 
link concept to track and obtain navigation data from 
a spacecraft that is normally not visible to Earth-bound 
tracking stations. In this problem, we again borrow from 
the Apollo mission profile and inspect the ability to track 
the CSM/LEM just after de-boost into lunar parking 
orbit by means of a communications satellite in a well- 
established lunar orbit which is visible to both the Earth 
and the CSM. 

As is shown in Fig. 41, the tracking station on the 
Earth broadcasts directly to the communications satellite, 
which in turn sends the received signal to the CSM; this 



Fig. 41 . Tracking geometry for 
behind-the-Moon tracking 


is retransmitted back to the satellite and finally to the 
original station on Earth. The resulting doppler shift on 
the return signal can then be equated to: 

d = + 2pci (2) 

where s is the station, c is the communications satellite 
and I refers to the CSM. This configuration could be 
used to monitor the de-boost maneuver itself, but this 
feature is not analyzed here. It was assumed that track- 
ing commenced immediately after the end of the bum 
and continued for the approximately 30 min that the 
spacecraft was behind the Moon. The spacecraft is in a 
circular orbit lying nearly in the lunar equatorial plane 
and has an altitude of 150 km, as is specified by Ref. 27. 
The communications satellite was placed in an orbit 
inclined 45 deg to the lunar equator with the longitude 
of the descending node at the zero meridian. Its radius 
was two lunar radii, and at the start of tracking was at 
— 90 deg longitude, —45 deg latitude. The relatively 
high inclination was specified in order that the out-of- 
plane components of the spacecraft s orbit would be 
clearly visible from the satellite. In this situation, only 
one data-type is possible, the modification of data-type 5, 
which is specified by Eq. (2). Again, data were taken at 
1-min intervals, and it was assumed the measurements 
were corrupted by 0.03 m/s rms noise. The a-priori 
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knowledge of the initial position and velocity of the 
communications satellite was assumed to be 0.1 km and 
0.05 m/s, respectively, in each of three axes. The corre- 
sponding numbers for the spacecraft were taken to be 
10 km and 1 km/s. 

In this problem, we are not interested in determining 
the relative orbits, as was the case in the rendezvous 
example, but in establishing the absolute orbit of the 
CSM/LEM combination. Fig. 42 plots the rms uncer- 
tainty in the three orbital elements of primary interest: 
(1) the semi-major axis, a, (2) the inclination, i, and (3) the 
radius of closest approach. Tea. 

It is seen that, after the end of 30 min, the uncertain- 
ties in these elements have been reduced to: = 4.5 km, 

cTi = 0.23 deg, and ^rca 17 km. The f^rcc, IS probably 
the most relevant of the three quantities plotted, for it is 
on the basis of the predicted closest approach altitude 
that a hypothetical safe orbit computation would be 
made. The performance shown for this quantity in Fig. 42 
is not as good as perhaps could be desired, but the data 



noise we have assumed is quite conservative and the 
a-priori uncertainty of the spacecraft's position is prob- 
ably needlessly large. Even as computed, an adequate 
safe orbit decision could be made after 15 min of tracking. 

It is interesting, for comparison, to assume for the 
moment that the Moon is transparent and that two-way 
tracking data from the spacecraft can be obtained. If 30 
min of two-way data alone are processed, the resulting 
o-r,ca is 28 km, a factor of 1.6 larger than previously at- 
tained. This degradation occurs in spite of the fact that 
the communications satellite has an orbit nowhere near 
as precisely determined as is the location of the tracking 
station (in fact, it was assumed here that the station loca- 
tion is precisely known). The superior tracking geometry 
afforded by the satellite appears to more than compensate 
for its imprecise location and velocity. 

If we now process both data-types, there is a dramatic 
reduction in <Jrca to 0.25 km, illustrating the eflFectiveness, 
in general, of viewing an orbit from more than one track- 
ing station. This also suggests that in the physically 
realizeable example, the orbit can be quickly determined 
very precisely as soon as the spacecraft becomes visible 
to the Earth. 

5 . An Interplanetary Bus/Capsule-Lander 
Problem, Example 3 

The third and last example deals with a situation that 
might occur on an interplanetary mission, either to Mars 
or Venus. The specific situation analyzed is that the main 
spacecraft (or bus) is travelling on a trajectory that will 
fly by the planet Mars. At 10 days out from encounter, 
a capsule is split off from the bus and is placed on an 
impacting trajectory with the planet. At the time of 
capsule ejection, the orbit of the bus is imperfectly known 
and the capsule's orbit, which has the same position error, 
is further corrupted in velocity by the execution errors of 
the ejection. Both craft are to be tracked until 1 day 
from encounter, at which time definitive orbits of both 
will be determined. 

For illustration, in this example, we shall capitalize on 
a feature of the multiple-link concept that has not been 
mentioned previously. Let us suppose that we are obtain- 
ing the three data types as previously listed, each of 
which is corrupted by noise: 

di = 2f)sji + Til 

do = 2psc + fi2 (3) 

ds — + pBC + Pcs + Tl5 
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where B and C refer to the bus and capsule, respectively. 
If we now form a new data type from a linear combina- 
tion of the first three, we may obtain: 

dn = ds — (di/2) ~ (^ 2 / 2 ) — Pbc + ns - l/2(ni + na) (4) 

An analysis of the origin of the data noise will reveal 
that it is primarily composed of (1) uncorrected atmos- 
pheric eflFects and (2) instability of the oscillator in the 
tracking station. 

It will be recognized that the signals used to obtain 
the three data types listed will traverse the same atmos- 
pheric path and have nearly the same propagation time. 
The major errors, and several others of less significance, 
will be nearly identical in all three sets of data; thus, 
Eq. (4) specifies a derived data-type equivalent to the 
inter-spacecraft range rate that should be unaffected by 
these major error sources. We shall use this technique, in 
this example, and assume that direct measurements of 
Pbc are made which contain smaller errors than either 
of the two-way doppler measurements made to the bus 
and the capsule. 

The bus was placed on a trajectory having a B distance 
of 8150 km;^ the components B • T and B • R were 5480 
and 6016 km, respectively. The components of the cap- 
sule were a B • R of 996 km and a B • T of 5926 km, 
which yields a trajectory inside the capture radius of 
6350 km resulting from the specified hyperbolic excess 
velocity of 3.16 km/s. Starting at 1 day past capsule 
ejection, the data were taken at 1 sample/hr with an rms 
noise of 0.004 m/s (this is effectively the same figure used 
in the two previous examples: the noise is assumed to be 
inversely proportional to the square root of the inte- 
gration time between samples) on the two-way measure- 
ments from the probe and the capsule. A noise % that 
figure was assumed on the psc measurement. 

The a^priority^ assumed on the bus injection conditions, 
was a pessimistic 1000 km in position and 1 m/s velocity, 
1 cr, in each of the three axes. The capsule was given 
identical position uncertainties but positively correlated 
with a coefficient of 0.98. A unity correlation is more 
realistic, but the orbits were uncoupled somewhat in 
their a-priori position so that the subsequent navigation 


^The B is the hypothetical vector from the center of the planet 
to the probe at closest approach assuming the planet has no mass. 
As such, B is perpendicular to the direction of the incoming 
asymptote, S. Normally B is broken into two components, B • R 
and B • T where T is a unit vector perpendicular to S and lying 
in the ecliptic plane and R is the unit vector which completes the 
RST orthogonal system. 


characteristics could better be studied. An additional 
5 m/s velocity uncertainty was added to the capsule s 
a-priority, and this was loosely correlated with the veloc- 
ity of the bus. 

Fig. 43 shows the navigation accuracy results obtained, 
as a function of tracking time, from the two-way mea- 
surements alone and then inclusion of the stereo data. 
Plotted here are the uncertainties in the determination 
of the B vector, 

<tb — ( o’b • R • t) 

for both the bus and the capsule. It is seen that the addi- 
tion of a stereo link improves the knowledge of the target 
parameters of the bus by only a very small margin. 
Initially, the capsule s orbit is improved by the addition 
of Pbc data, in fact, the uncertainty falls by a factor of 3 
and quickly becomes nearly identical to the bus uncer- 
tainty, but by the end of the 8 days of tracking, this gain 
has disappeared. 



Fig. 43. RMS error in prediction of the miss parameters 

Fig. 44 explains this behavior by plotting the rms dif- 
ference errors in the estimation of the miss parameters, 
o-AB. Without the psc data, falls to 300 km at the end 
of the tracking period; the corresponding figure using 
pBc is 5 km, or a factor of 60 improvement. Even with 
only 2 days of tracking, the relative B estimation errors 
are only 70 km. Thus, the conclusion is reached that the 
inter-spacecraft link results in establishing the relative 
orbits of the two vehicles almost immediately even when 
the probes are experiencing rectilinear motion. Referring 
back to Fig. 43, the lower ctb on the capsules absolute 
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orbit is achieved by passing through the well-determined 
bus to capsule orbit, the better knowledge of the orbit of 
the bus. 

This leads us to observe that with an initially well- 
determined bus/capsule orbit, the capsule could be 
ejected and its new orbit quickly redetermined to the 
accuracy of the parent spacecraft. 

6 . Conclusions 

On the basis of the three numerical examples presented, 
it may be concluded that establishing an inter-spacecraft 
doppler link can be very eflFective in determining the 
relative orbits of two vehicles. Even with otherwise very 
weak tracking geometry (Example 3), precise relative 
orbits can be determined with dispatch. In Example 1, 
the stereo data were found to be somewhat more power- 
ful than three-way data at lunar distances. This was true 
even when the relative inclination between the two orbits 
was 0.5 deg; with larger relative inclinations one would 
expect a further improvement in the stereo data's ability 
to determine the out-of -plane components of the orbit. 



Fig. 44. RMS difference in the prediction 
of the miss parameters 


However, it seems equally clear that, at least over short 
data arcs or in weak tracking geometry situations, the 
absolute orbits are not significantly improved with the 
addition of the multiple-link data. Thus, the usefulness 
of this, or related techniques, is only particularly attrac- 
tive when: (1) the orbit of one vehicle is accurately 
known, (2) the objective of the tracking is to establish 
the relative orbits, as in Example 1, or (3) direct Earth- 
based tracking is unobtainable, as in Example 2. 

Acknowledgment. The author is indebted to Dianne 
Wimp for providing the programming and computing 
support needed to obtain the numerical results presented 
in this article. 


,167 

K. Pioneer Project Support 

D. Curkendali 


The Pioneer Project consists of a series of probes whose 
purpose is to gather scientific information concerning the 
deep space environment in the region near the ecliptic 
plane and at a distance from the Sun varying from 0.8 
to 1.2 AU. The DSN is committed to track the Pioneer 
spacecraft, obtain engineering and scientific telemetry, 
determine the orbit of each probe, and provide facilities 
within the SFOF and DSIF to the project for the con- 
duct of the mission and reduction of data. Specific items 
under this task will be selected each reporting period for 
discussion in the SPS. 

For this issue, Motsch, Chaney, and Curkendali have 
contributed an article ( Section L ) analyzing a sinusoidal 
periodicity discovered in the Pioneer VI doppler tracking 
data. It is found that this periodicity arises from a slight 
offset of the spacecraft antenna from the axis of rotation 
This phenomena is used to make a very accurate deter- 
mination of the spacecraft rotation rate and, incidentally, 
to determine magnitude of the offset itself. 

A second article (Section M) by Gallagher presents 
an analysis which relates injection dispersions due to 
launch-vehicle guidance errors to the initial tracking sta- 
tion s acquisition. Two of the parameters most critical to 
early acquisition are selected, and upper bounds on the 
variations of these parameters with injection errors are 
computed. 
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Lastly, Ball presents (in Section N) the second of a 
series of articles analyzing and tracing the efFects of the 
observed Pioneer VI gas leak on the orbit. For this 
analysis, the tracking data were obtained from March 1 
to June 20, 1966. It is concluded that the gas leak, while 
diminishing, is still significantly aflFecting the spacecraft 
orbit determination. 
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L Pioneer VI High-Frequency 
Data Analysis 



■►EARTH 


Fig. 45. Spacecraft geometry 


R. Motsch, D. Chaney^ and D. Curkendall 


1. Introduction 

During the analysis of the Pioneer VI high-resolution 
tracking data, a periodicity was discovered in the data 
noise. This periodicity arose from a slight offset of the 
spacecraft antenna from the spin axis. The resulting 
simple harmonic motion in the station— spacecraft range 
dominated the autocorrelation function of the noise; thus, 
it was necessary to remove the periodicity in order to 
determine the amount of random phase jitter in the dop- 
pler signal and the high-frequency drift characteristics 
of the rubidium oscillators. 

Analysis of the problem showed that, in the process of 
removing the periodicity, one could obtain two param- 
eters of interest: (1) the spacecraft rotation rate, and 
(2) the distance of the spacecraft antenna from the spin 
axis. This analysis, when applied to the Pioneer VI data, 
led to: 

(1) Measurement of the approximately 58-rev/min 
spacecraft rotation rate and the change in this 
rate with a minimal accuracy of 0.01 rev/min. 

(2) Measurement of an apparent antenna offset of 
approximately 1.5 cm. 


2. Spacecraft Rotation and Antenna Offset 
Analysis 

The spacecraft geometry pertinent to this study is 
portrayed in Fig. 45. We have considered the spacecraft 
antenna as a point source located a distance d off the 
spin axis of the spacecraft. 


It can be shown that, for a sampled data system, the 
autocorrelation function of the range-rate error due to 
antenna displacement and spacecraft rotation is given by: 

R(r) = 2d^ sin^ (<o/2fa) cos <ot 

where d is the apparent offset of the antenna when con- 
sidered as a point source, <a is the angular velocity ( rad/ 
sec) with respect to the spacecraft spin axis, and fd is 
the data sampling rate. This expression for R(t) was ob- 
tained under the assumption that, over any one pass, <a 
is essentially constant; i.e., 

2d^ sin2 (o./2/d) = A 

jR(t) = a COS<t>T 

Consider Fig. 46, where R(r) is shown as the dashed 
curve. If the data sampling rate is more than twice the 
frequency of the autocorrelation curve, R(t) can be con- 
structed directly. However, if the data sampling rate is 
less than twice the frequency (rf the autocorrelation curve, 
some pseudo-autocorrelation R(t) will be obtained with 
amplitude of A, but with a period greater than that 



T 


Fig. 46. Amplitude v$ time for the true autocorrelation 
function R(t) and the pseudo-autocorrelation 
function R(t) 
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of R(t). This situation is demonstrated in Fig. 46 by the 
unprimed numbers and the solid curve. 

Let fs be the spacecraft rotation rate and fa, the data 
sampling rate. Then, for there exists a time in- 

terval T such that: 

fJ±l^faT 

Thus, every T sec it will appear that the spacecraft has 
gone through one rotation, as measured by the period 
of R(t) using a data sampling rate fa. From the above 
equation, we obtain: 

Since T is necessarily positive, the above equation can 
be written: 

n==l,2,--- 


It has been shown above that any spacecraft rotation 
rate that is less than one half the data sampling rate 
can be determined uniquely, and any other spacecraft 
rotation rate that is not an integral multiple of the data 
sampling rate can be restricted to several possible choices. 
An example pertaining to these statements is shown in 
Fig. 47. In general, one will need a priori information to 



Fig. 47. Possible spacecraft rotation rates for dafti 
sampling rate of 60 samples/min 
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Fig. 48. Range residuals, normalized autocorrelation function, and normalized filtered autocorrelation function 

for Pioneer VI on January 5, 1966, using Echo DSS 12 data 
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resolve the ambiguity in the pseudo-autocorrelation R{t). 
If such information is available, can be uniquely deter- 
mined, and, since w = 2nrfs, the apparent ofiEset d can be 
obtained from the equation: 

2(P sin2 (o>/2fa) = A 

where A is the autocorrelation amplitude. 

3. Results from Pioneer VI Data Analysis 

The techniques described above were applied to data 
collected from the Pioneer VI probe. The Earth rate was 


removed from the data, a seventh-order polynomial was 
fitted to the resulting numbers, and residuals about the 
fit were formed. From these residuals, the autocorrelation 
function was produced, and a curve of the form A cos {(at) 
was fit to the autocorrelation data. Fig. 48 displays the 
range residuals, the normalized autocorrelation function, 
and the normalized filtered autocorrelation function ob- 
tained using the above processes on data for January 5, 
1966. Information from a spin counter on the spacecraft 
was utilized to uniquely determine /, from <o, and d 
was determined from A. Table 3 lists the results for 
several passes of Pioneer VI, 

Table 3. Results of Pioneer Vi data analysis 
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Fig. 49. Pioneer VI rotation rate vs date 1966 


Date 

utf rad /sec 

A, cps‘xi0 ** 

fi, rev/min 

d, cm 

Computed 

Measured 

1- 5-66 

0.1518 

5.85 

58.555 

58.40 

1.47 

1- 5-66 

0.1513 

7.97 

58,550 

~ 

1.72 

1-30-66 

0.1909 

10.93 

58.177 

58.04 

1.60 

2- 8-66 

0.2003 

11.75 

58.087 

- 

1,58 

3- 4-66 

0.2141 

15,92 

57.956 

57.72 

1.72 

4- 9-66 

0.2292 

17.95 

57.811 

- 

1.71 

4-21-66 

0.2305 

18.37 

57,799 

- 

1,72 


It is apparent that the Pioneer VI rotation rate is de- 
creasing. This rate decrease is plotted in Fig. 49. It is 
known that Pioneer VI experienced a gas leak, and quite 
probably the decreasing rotation rate is associated with 
this leak. 
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M. Pioneer 6‘ Trajectory 
Dispersions and DSN 
Initial Acquisition 

J. F. Gallagher^ 


I . Introduction 

During mission operations, it is important to the orbit- 
determination process that the spacecraft be acquired 
in two-way lock as soon as possible. This is true because 
the slant-range rate from Johannesburg DSS 51 during 

^Pioneer B was launched on August 17, 1966, at which time its 
designation beca me Pioneer VII. ^ 

^The author wishes to express his appreciation to D. Hubiak, W. 
Ross, and F. Borncamp of JPL Technical Section 337 for infor- 
mation furnished on the DSN. 
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the first 10 or 12 min after injection varies rapidly, and the 
doppler data during this early period is therefore power- 
ful. The slant-range-rate curve from DSS 51 for the 
nominal Pioneer B trajectory is shown in Fig. 50. 

For DSS 51 to acquire as early as possible, the station 
must be prepared to track the spacecraft on not only a 
nominal trajectory but also various off-nominal trajec- 
tories which might occur. The off-nominal trajectories 
referred to here are those resulting from such things as 
3a-high or -low thrust performance of the third stage, but 
not ‘‘catastrophes” such as a failure of the third-stage 
motor to ignite. 

There are two important quantities involved in early 
acquisition by DSS 51: rise azimuth, and the time at 
which the slant-range acceleration from DSS 51 first 
falls within the doppler capability of the spacecraft and 
ground station. This time is the earliest time at which 
initial RF acquisition for two-way lock can occur. The 
doppler capability of the spacecraft and ground station 
is rather complicated, depending upon many parameters, 
including bandwidth, sensitivity, and selectivity of the 
ground station and spacecraft receivers. This doppler 
capability also changes with time as improvements are 
made in the equipment. Personnel of JPL Technical Sec- 
tion 337 chose as a working definition of the doppler 
capability for the Pioneer B mission the following: the 



Fig. 50. Slant-range rate from DSS 51 for the 
nominal Pioneer B trajectory 


time at which the slant-range acceleration from DSS 51 
falls below 40 m/ sec^. 

Knowledge of the rise azimuth and its dispersions is 
important because, if these dispersions fall outside the 
beamwidth of the acquisition-aid antenna, a more diffi- 
cult and time-consuming search for the spacecraft may 
be necessary and much early data valuable to orbit de- 
termination may be lost. On the other hand, if the dis- 
persions in rise azimuth fall well within the beamwidth 
of the acquisition-aid antenna, the chances for very quick 
two-way acquisition are greatly improved. Knowledge of 
the time at which the slant-range acceleration from 
DSS 51 falls below 40 m/sec^ and the dispersions on 
this time is important because with this information it 
is possible to work out a procedure for tuning the trans- 
mitter and receiver which will improve the chances of 
quick two-way acquisition and thereby maximize the 
amount of strong doppler data. 

It is sufficient for the DSN to have advance knowledge 
of the nominal trajectory and two off-nominal trajectories 
which are loosely termed 3<T-high and 3<r-low. These 
off-nominal trajectories are those which result from a 
“worst-case” combination of the dispersions in each of 
the six injection conditions. In this case, the term “worst 
case” means worst with respect to launch azimuth and 
the time at which the slant-range acceleration curve falls 
below 40 m/sec^. 

2. Method 

To determine these worst-case trajectories, the follow- 
ing method was used: First a nominal trajectory was 
computed beginning at injection and running out to 2 hr 
past injection. All of the station quantities for DSS 51 
were obtained for this nominal trajectory. Then the 3 <t 
dispersions in each of the six Earth-fixed spherical injec- 
tion conditions were applied, one at a time, to the cor- 
responding nominal injection conditions; the trajectories 
were computed to 2 hr; and the DSS 51 quantities were 
obtained along the trajectory. Trajectories using both 
positive and negative increments in the injection condi- 
tions were calculated to check for linearity. This pro- 
cedure allows one to see the effect of a dispersion in any 
one of the six injection conditions upon any station 
quantity. A worst-case combination of injection condi- 
tions can be gotten simply by picking that set of plus-and- 
minus increments in the six injection conditions which 
results in all-positive or all-negative changes in the sta- 
tion quantity. If linearity holds, then we are guaranteed 
that the combination thus arrived at does, in fact, result 
in the greatest change in the station quantity, and that 
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we have a pair of oflE-nominal trajectories which brackets 
the trajectories we can actually expect to track. 

3. Results 

The Earth-fixed spherical injection conditions and their 
3<r dispersions are shown in Table 4. One at a time, 
these 3<7 dispersions were added to and subtracted from 
their corresponding nominal quantities, and the trajec- 
tories were calculated. The results for the rise azimuth 
at DSS 51 are shown in Table 5. Here we see that the 
smallest change in rise azimuth at DSS 51 is caused by 
the 3<r dispersion in V and the largest change, by the S<r 
dispersion in a, the injection azimuth. The worst-case 
combination of injection conditions to maximize the value 
of the rise-azimuth angle at DSS 51 would be obtained 
by using positive increments in R, <f>, V, and y and a 
negative increment in <7. Increments with opposite signs 
would be used to minimize the rise azimuth, and this 


Table 4. Earth-fixed spherical injection conditions 
and their 3a dispersions 


Condition 


3<r Dispersion 

Geocentric radius R, km 

6791.6385 

22.5 

Geocentric latitude 0 , deg 

-14.052728 

0.258 

Longitude deg 

352.65549 

0.384 

Earth relative velocity V, km /sec 

10.546620 

0.075 

Relative path angle 7 , deg 

2,5642066 

1.05 

Relative azimuth angle <r, deg 

131.08955 

1.05 


pair of trajectories should bracket the trajectory actually 
flown, at least as far as rise azimuth is concerned. Refer- 
ring to Table 5, it appears that a rise-azimuth spread of 
3 deg on either side of the nominal value should cover 
virtually all of the possibilities. Since the beamwidth 
of the acquisition-aid antenna is about 16 deg, the 
rise-azimuth dispersions should cause no trouble in 
acquisition. 



Fig. 51. Slant-range acceleration from DSS 51 
for the nominal Pioneer B trajectory 


Table 5. Effects of injection-conditions dispersions 
upon rise azimuth at DSS 51 


Trajectory 

Resulting rise azimuth, deg 

Difference, deg 

Nominal 

270,47687 



+ AR 

270.87919 

0.40232 

“AR 

270.05910 

-0.41777 

+ A<p 

271.11444 

0.63757 

-A 0 

269.83581 

“0.64106 

+ Ar 

269.65577“ 

-0.82110 

-A0 

269.96299 

“0.51388 

+ AV 

270.50748 

0.03061 

-AV 

270.44529 

-0.03158 

+ Ay 

271.05562 

0.57875 

—Ay 

269.83356 

-0.64331 

+ A<r 

269.60717 

-0.86970 

— A<r 

271.33660 

0.85973 

“The input for the +A0 trajectory was incorrect; however, a check for 
linearity showed that the variation in rise azimuth with injection longi- 
tude is very linear for the small changes in the longitude used. 


Fig. 51 is a plot of the slant-range acceleration of the 
probe as seen from DSS 51. It can be seen that the slope 
of the curve is very steep in the neighborhood where the 
curve falls below 40 m/sec^. The trajectory program does 
not compute the time at which the slant-range accelera- 
tion from a station takes on a given value. However, for 
the purposes of the study, it was sufficient to assume a 
straight-line variation in the slant-range acceleration be- 
tween 7 and 8 min after injection. To find the time of this 
crossing, linear interpolation was done graphically for 
the nominal trajectory and for each of the six pairs of 
off-nominal trajectories. Fig. 52 is a sample of this graphi- 
cal linear interpolation for the dispersions in the crossing 
time caused by dispersions in flight-path angle. Observe 
that the line labelled ‘ nominaF in Fig. 52 is a linearized 
and magnified version of the portion of the curve of 
Fig. 51 which lies between 7 and 8 min after injection. 

The effects of the dispersions in the six injection condi- 
tions upon the time of crossing the 40-m/sec^ level are 
given in Table 6. In Table 6, we see that dispersions in y 
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Fig. 52. Sample plot showing the graphical linear 
interpolation for the dispersions in crossing time 
caused by dispersions in flight path angle 

are most influential here. The worst-case combination of 
injection conditions which causes the crossing time to be 
as early as possible would be a positive increment in 
R, negative in and positive in 0, V, y and <r. This tra- 
jectory and its opposite (i.e,, the one for which the 
crossing time is latest) should bracket the actual trajec- 
tory as far as the crossing time is concerned. From this 


information, the tuning procedure mentioned earlier can 
be worked out. 


Table 6. Effects of injection-conditions dispersions 
upon the time of the slant-range acceleration 
curve falling below 40 m/sec^ 


Trajectory 

Resulting crossing time, 
minrsec after injection 

Difference, sec 

Nominal 

07 : 25.0 

- 

+ AR 

07 : 24.2 

- 0.8 

-AR 

07 : 25.7 

0.7 

+ A<p 

07 : 26.8 

1.8 

-A(p 

07 : 23.0 

- 2.0 

+Ar 

07 : 28 .r 

3.1 

-AS 

07 : 26.9 

1.9 

+ AV 

07 : 23.4 

- 1.6 

~AV 

07 : 26.5 

1.5 

+ Ay 

07 : 19.7 

- 5.3 

-Ay 

07 : 29.7 

4.7 

+ A<r 

07 : 21.5 

- 3.5 

— Aff 

07 : 28.2 

3.2 

«The input for the trajectory was incorrect; however, a check for 

linearity showed that the variation in crossing time with injection longi- 
tude is very linear for the small changes in the longitude used. 


J. E. Ball 

Pioneer VI, launched on December 16, 1965, has been 
in a heliocentric orbit for the last 7 mo. The one anomaly 
that occurred during this time was a gas leak somewhere 
in the pneumatic system of Pioneer VI. Since 90% of the 
gas carried on the spacecraft was depleted 90 days after 
launch (March 16, 1966), JPL’s single-precision orbit- 
determination program (SPODP) attempted to fit a tra- 
jectory over the arc of the spacecraft orbit from March 1 
to June 20, 1966. Fig. 53 is a plot of the difference be- 
tween the observed Echo DSS 12 two-way doppler values 
and the calculated two-way doppler values (residuals) 
for a portion of that arc. These residuals have a similar 
signature to the residuals of the fit over the arc of the 
spacecraft orbit from December 1965 to January 1966 
(SPS 37-39, Vol. Ill, pp. 42-45). However, in this instance, 
the data covers approximately 4 mo. Earlier in the mission. 
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the data arc could not be extended much more than 1 mo 
before excursions in the residuals considerably exceeded 
those shown here. Thus, our relative ability to fit the 
data is improving with time, as would be predicted from 
the fact that the gas flow rate is following a decreasing 
exponential curve. 

Table 7 compares the values of the spacecraft s posi- 
tion and velocity for December 18, 1965, as computed 
from the fit over the arc of the orbit from December to 


January to the values as computed from the fit over the 
arc of the orbit from March through June. The change in 
the injection conditions obtained is larger than the theo- 
retical uncertainties computed, assuming no gas leak to 
be present. This may be seen by comparing the values 
of aQ in Table 7 to the values for <r. The version of the 
SPODP used to fit portions of the Pioneer VI orbit did not 
have the capability of estimating an exponential type of 
force such as that occurring in Pioneer VI from the gas 
leak. 







TRACKING TIME, min 

Fig. 54. Residuals vs tracking time for a portion of the Pioneer Vi trajectory fit for March through June 1966, 
using estimated values for the AU and the deep space stations^ geocentric radius and longitude 
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Table 7. Comparison of Pioneer VI position and velocity as computed from 
two fits of data (Epoch: December 1 8, 1 965) 


Q 

Fit from December 1965 
to January 1966 

1 <r Value 

Fit from March through 
June 1966 

I<r VoliM* 

AQ 

X,km 

430843.04 

78 

430809.56 

82 

33 

Y, km 

-24759.184 

120 

-24011.667 

121 

748 

Z,km 

-40436.293 

240 

-42300.775 

233 

1864 

X, km/sdc 

2.1548440 

0.0001 

2.1566785 

0.00009 

0.00016 

Y, km/ sec 

0.02582071 

0.00002 

0.02623449 

0.00006 

0.00041 

Z, km /sec 

-0.10912526 

0.00007 

-0.11013935 

0.00015 

0.00101 

"These uncertainties are 

the values calculated, assuming no gas teak to be present. 




In addition to the information desired from the prime 
experiments of Pioneer VI, it was hoped that, from analy- 
sis of the Pioneer VI orbit, a better determination of cer- 
tain astronomical constants could result. Therefore, a 
determination of the astronomical unit (AU) was at- 
tempted over that portion of the orbit from March 
through June when the gas leak was small. In addition to 
solving for the AU, the deep space stations' geocentric 
radius and longitude were estimated to help improve 
the fit, as shown in Fig. 54. 

Since the residuals from Fig. 54 are smaller and dis- 
tribute themselves in a linear fashion along the zero 
residual line, as compared to the residuals in Fig. 53, 


improvements to the AU and deep space station locations 
are suggested when the additional parameters of the AU 
and station locations are estimated. However, a com- 
parison of these values to the values determined from 
previous sources (Table 8) indicates that small unmolded 
and unknown forces acting on the spacecraft can pro- 
duce erroneous station location results and other param- 
eter estimation errors. Thus, it is concluded that the 
effects of the Pioneer VI gas leak, while diminished, are 
still significantly affecting the orbit of Pioneer VL Future 
plans include continued processing of the tracking data 
as the gas leak continues to diminish and attempts to 
model the leak in the SPODP to obtain reliable estimates 
of the AU and certain other astronomical constants. 


Table 8. Comparison of values for the AU and tracking station longitudes 
as estimated from Pioneer VI orbit to those estimated previously 


Q 

Estimates from Pioneer Vf orbit 

1<r Value 

Previous estimates 

1(T Value 

AO 

AU, km 

149,598,760 

957 

149,598,500 

500 

266 

Deep space station longitudes 
Echo DSS 12 

243.19530 deg 

_ 

243.19428 deg 

40 m 

102 m 

Tidbinbilla DSS 42 

148.98208 deg 

- 

148.98132 deg 

40 m 

76 m 

Johannesburg DSS 51 

27.686081 deg 

- 

27.685513 deg 

40 m 

57 m 

Robledo DSS 61 

355.75172 deg 

- 

355.75083 deg 

40m 

89 m 
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III. Communications Research and Development 
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A. Wide Band Receiver System 
Frequency Divider 

K. Schreder 

7 * Introduction 

The wide band receiver system (WBRS) is currently 
in the module checkout phase with less than 30% of the' 
modules remaining to be checked. One of the modules 
designed, built, and checked out is a frequency divider 
from 66 MHz. The module is shown in the block diagram 
of the transmitter section of Fig. 39, p. 101, SPS 37-33, 
Vol. III. 

2 . Circuit Description 

The circuit is a locked oscillator divider. Fig. 1 shows 
the circuit configuration. The input circuit is designed 
for a 50% bandwidth to provide a low input voltage 
standing wave ratio (VSWR) over a wide 1-db band- 
width. The first stage, Ql, is an emitter-follower stage 
which is necessary for maintaining a low VSWR for 
greater than a 5% bandwidth. The next stage, Q2, pro- 
vides limiting and insures proper isolation between the 


divider and the emitter follower. The divider stage, Q3, 
is a regenerative divider whose output circuit is tuned 
to the output frequency. Minimum power is taken from 
the divider which allows the divider to remain locked 
over a wide power range. The output stages, Q4 and Q5, 
amplify the signal to an output power greater than +13 
dbm. (The FT-45 transistor was used instead of the more 
common 2N915 or 2N918 because of the higher current 
gain cutoff frequency The FT-45 has a lower current 
gain than the other two transistors, but a gain- 
bandwidth product Ft, about the same as a 2N918 which 
produces the higher F^,) 

3. Performance Characteristic 

For a -|-13 dbm input signal, the output power was 
-|-13.6 dbm for a power drain of 600 mw from a -|-30 
vdc supply. The input VSWR was less than or equal to 
1.10 for a 5-MHz bandwidth and input power variations 
from -\-l to -|-19 dbm. 

The locking bandwidth was greater than 7.4 MHz. The 
percentage of locking bandwidth was 11.2%. (The lock- 
ing bandwidth is defined as the bandwidth in which the 
output signal is synchronous with the input signal.) 


PRECEDING PAGE ELA!-4K NOT FiLVitD. 
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The limiting curve is shown in Fig. 2. When the input 
signal is below —18 dbm, the locked oscillator stops. For 
a ±10% power supply variation, the output power varied 
±0.9 db. Fig. 3 shows the output power variation for 
power supply variations from 40 to 14.5 vdc. For a supply 
voltage of 14.5 vdc or less, the oscillator stopped. The 
total output harmonic distortion was 1.88% with the 9th 
harmonic being the last measurable product. 



INPUT POWER, dbm 


Fig. 2. Limiting curve for h- 3 frequency divider 


. 

ns Contro 




B. Ephemeris 
Oscillators: Modification of 
Mod V Programmed Oscillator 


G. Thompson 


The Mod V programmed oscillator (PO) at the Mars 
DSS has been modified to increase its reliability. Because 
of this modification, the Mod V PO gain constant Ko, is 
no longer unity when used in a closed-loop, phase-locked 
receiver system (Fig. 4). 



TOTAL LOOP GAIN /T = MK^yK^ 


Fig. 4. Phase-locked loop receiver using a PO 


The modification involves the 24-MHz voltage- 
controlled oscillator (VCO) as used in conjunction with 
the Hewlett-Packard (H-P) 5100 A/51 10 A synthesizer. As 
previously reported (SPS 37-32, Vol. Ill, pp. 36-38) the 
24-MHz VCO provides a means of voltage controlling 
the output frequency of the H-P synthesizer at frequen- 
cies from 1 to 50 MHz. The original system design (SFS 
37-32) substituted the 24-MHz VCO for the 24-MHz ref- 
erence in the 100-kHz decade of the H-P synthesizer 
(Fig. 5). This configuration made the VCO-synthesizer 


E 




Fig. 3. Output power variations due to power 


supply changes 


Fig. 5. Original 24-MHz substitution in Mod V PO 
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combination have the sensitivity of the VCO, K^, with 
Ko = 1 (Fig. 4). 

The modification in Fig. 6 substitutes the 24-MHz 
VCO for the 24-MHz reference from the 5110A driver. 
Therefore, voltage control of the frequency can occur in 
the eight least significant digits of the synthesizer out- 
put frequency. The gain of the VCO is multiplied by 
1.1111111, or in other words, Kq = 1.1111111. 

The modification was made because problems had been 
found with the substitution of the VCO in the 100-kHz 
decade. The modification requires a 24-MHz VCO oper- 
ating at 100-mv rms into 50 O. The modification is simple 
to implement and provides a more reliable system. 



Fig. 6. Present 24-MHz substitution in Mod V PO 
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C. Information Systems: Error 
Correction for High-Speed 
Data Transmission 

E. C. Posner and G. Solomon 

1 . Mariner Spacecraft High-Speed Data 
Transmission, E. C. Posner 

a. Introduction. This report describes a coding- 
decoding scheme suitable for use with Mariner 1964, 
1967, and 1969 spacecraft data transmitted over a high- 
speed data line. The method uses a (15,7) double error- 
correcting code and a very simple majority logic decod- 


ing. The probability of failing to decode correctly is 
reduced by a factor of more than 200. 

Data from Mariner 1964, and ultimately from the 
Mariner 1967 and 1969 spacecraft, come to the SFOF 
over a 1200-bps commercial grade high-speed data line. 
The bit rate for Mariner 1964 spacecraft is at most 25 
bps; for Mariner 1967 spacecraft is at most SSVs bps; and 
for Mariner 1969 is at most 66%bps. Of the 1200 bps, 
450 bits are sync bits required by the telephone company 
and are not at the disposal of the user. Thus, 750 bps are 
all that are available to the Ground Communications Sys- 
tem (GCS). Since, at most, 10 data words (a data word has 
7 spacecraft information bits) arrive per second, the chan- 
nel has space for other bits. These other bits are timing 
bits, station AGC bits, etc. There are 150 bps left for the 
spacecraft data. Figuring 10 spacecraft data words per 
minute, each spacecraft data word can be allotted 8 bps 
for error-correcting purposes, yielding coded words of 
length at most 15. 

Fig. 7 shows the present method of error correction. 
An eighth data bit is put into the block — a sync informa- 
tion bit having nothing to do with the spacecraft data. 
These 8 bits are blocked into a 4 X 2 rectangular array. 
A third column is added as even parity on the four rows. 
A fifth row of three bits is now added as vertical parity 
on the three columns. The resulting 5x3 array of 15 
bits is transmitted in a top-to-bottom, left-to-right serial 
scan. 

This (15,8) coding scheme yields a code of minimum 
distance 4, which can correct single errors and still detect 
double errors. The purpose of this note is to propose a 
(15,7) code of minimum distance 5 which can correct 
double errors. The coding of more than one word at a 
time is not considered acceptable at the present time for 
the Mariner missions. 
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fe. The Code, We prox>ose that the special sync bit S 
in Fig. 1 be moved to another portion of that 75-bit 
block which comes in with every word received from the 
spacecraft. For example, there is space in the portion of 
the block reserved for station AGC. Thus, we can add 8 
check bits to the 7 spacecraft data bits to obtain coded 
words of length 15. 

We can now use a (15,7) code instead of a (15,8) code. 
The (15,7) code which suggests itself is the famous 
double-error-correcting quasiperfect Bose-Chaudhuri 
(15,7) code. The definitions need not concern us here; the 
important point is that the code has minimum distance 5, 
and thus can correct all double errors, instead of just all 
single errors. The code is generated from the recursion 
polynomial = x® + x + 1, addition being modulo 2. 
This will be explained in the next section. A (14,7) 
double-error-correcting code does not exist (Ref. 1). 

c. Encoding, Enct)ding would be done on the SDS 920 
TDH computer at the DSIF sites. This high-speed data 
line will not go through the NASCOM communications 
processor at the Goddard Space Flight Center, and so 
no special formatting problems will arise. (This contrasts 
with the case for teletype. See SPS 37-35, Vol. Ill and 
Sect. D following. 

To encode, regard the seven information bits Oo, Ui, (h, 
^ 3 , ^ 5 , as elements of the 2-element field. That is, 

each is 0 or 1, addition being modulo 2. The 8 check 
bits are generated by the rule 

^i+7 — ^i + 3 “h ^4-1 4" 0 ^ ^ 7. (1) 

Extra storage space required in the 920 is minimal (less 
than ten words), and the running time is less than a milli- 
second. 

Because the polynomial x^ + ^ + 1 divides x^® + 

1 over the 2-element field, it can be shown that the bits 
i > 0, if computed by the rule in Eq. (1), would 
just repeat Co through in the same order. This fact is 
needed in the decoding process. 

d. Decoding, We present a ‘majority logic” decoding 
scheme which is especially easy to implement on an IBM 
7044, the computer used as the SFOF communications 
processor. Little extra storage (40 locations at most) is 
required to accommodate the (15,7) decoding scheme, 
and the running time is a few milliseconds. 

The relation in Eq. (1) implies that Ut ^ -\- Ui Uo, 

or better, flo = + Us + Let i = 8 in Eq. (1) to ob- 


tain ai 5 = ail + ^9 + ^8* Since as observed flis = Co, we 
have flo = ^ii + ^9 + Similarly, by letting t = 12 and 
i = 14 in Eq. (1) we obtain Uo = ^*4 + (h 2 +tfi 3 and ao ^ 
+ ^6 + ^ 14 * (III f^ct, these are the only check equa- 
tions involving Uq.) 

Consider the five expressions below, which are all 
equal to the original transmitted bit ao in the absence 
of errors: 


4- Us 4“ u? 

U2 4" ^6 4“ ^14 (2) 

U4 4" ^12 “h ^13 
^ Ug 4” U9 -f- flu . 

Compute the five quantities in Eq. (2) from the received 
word. If two or fewer errors have been made, at most 
two of the five quantities are changed, since no index 
occurs in more than one expression in Eq. (2). Hence, 
decode flo as do, where ^ is 0 if the majority of the five 
expressions of Eq. (2) are zero; similarly, do is 1 if the 
majority are 1. 

By the cyclic property of the (15,7) code, we can cycle 
the received word around one to the left, and decode 
di as ai in the same way. That is, the original fli becomes 
the new Aq after the shift, the original U 2 becomes fli, • • * , 
ai 4 becomes flis, whereas ao becomes flu. By repeating 
this process out to flg, we can decode all the information 
bits. In fact, this shows that the (15,7) code is double- 
error correcting. 

However, we always produce output in this fashion, 
even when more than two errors are made. Since the 
(15,7) code does not always make an error if more than 
two errors are made, we wish to detect whether such 
detectable errors have occurred. To do this, generate an 
fl 7 , • • • , flu by the recursion of the previous section. Count 
the number of disagreements of the word do, • * • , du with 
the original received word. If the number of disagree- 
ments is 2 or less, we assume a double error, single error, 
or no error has been made, and declare the decoded 
information bits do, • • • , do as correct. If, however, there 
are 3 or more disagreements, then we declare that a 
detectable but uncorrectable error has occurred. 

e. Performance improvement. The main improvement 
in performance comes from the fact that double errors, 
which are merely detected by the (15,8) double-parity 
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code, are corrected by the (15,7) cyclic code. Thus, when 
the bit error probability in the channel is 10~®, which is 
hopefully the worst encountered, the (15,8) code yields 
as output a detected but uncorrected error with proba- 
bility about 



that through simple modification, greater and more effi^ 
cient use of the projected coding system can be obtained. 
In particular, most of the data to be discarded can be 
recovered and used. In noisy situations, this will mean 
the difference between having a good channel and hav- 
ing no channel at all. The data is ordinarily spacecraft 
telemetry data from the Deep Space Stations received at 
the SFOF. 


The (15,7) code has for this probability 

15 ^ _i8oJ (10-®)* 2.75 X lO"' 

The 180 represents the number of triple-error patterns 
which are confused with double errors and decoded in- 
correctly. Thus, the mean time between failures to pro- 
duce correct decoded output is increased by a factor of 
220, from 2.8 hr to 26 days. This is a worthwhile im- 
provement. 

The probability of error (probability that the decoder 
outputs a word in error) is low in both systems, the pres- 
ent (15,8) and the proposed (15,7). For example, when 
the input bit error probability is 10“^, and burst errors 
are not occurring, the output word error probabilities are 
less than 2 X 10“'^. Burst errors could be avoided by dis- 
persing the 15 bits of the code word throughout the 
transmitted block of 75 bits. 

A computer program is now being written to incorpo- 
rate this coding-decoding scheme in the Mariner space- 
craft format for the CCS. 
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D. Information Systems: Error 
Detection Equipment 

G. Solomon 


The basic code used by the NASCOM is a shortened 
(1023,990) BCH code with minimum distance 8, giving 
the code a capability of detecting up to 7 errors, correct- 
ing 1 and detecting 6, correcting 2 and detecting 5, or 
correcting 3 errors and detecting 4. As furnished, the 
equipment will detect errors and add to the received 
word of length 597 a 3-bit indicator of error status, i.e., 
it will tell whether the word was correctly received or 
not. The incorrect words would not be decoded. How- 
ever, a computer program can be set up to correct up to 
3 errors in the received word via an algebraic mechaniza- 
tion — a variant of Berlekamp-Solomon decoding proce- 
dures for cyclic codes (1) which is capable of decoding 
the word as fast as the word is received. As a sample 
result, for a probability of bit error of 0.001, which is 
considered very noisy for high-speed data lines, a double 
error-correcting procedure can expect to recover 0.95 of 
such words which are detected in error; a triple error- 
correcting procedure will recover 0.992 of these words. 
This triple error-correcting code will still detect 4 errors; 
the probability of letting an undetected error through in 
the block of 567 information bits would be less than 10 

Without error correction, 0.45 of the received words 
will be in error and would not be used; by using the 
triple error-correcting, quadruple error-detecting proce- 
dure as proposed, only 0.003 of the words would be 
thrown out. The output word error probabilities for bit 
error probabilities of 0.005, 10~^, 10~^ and 10~''’ are 
given in Tables 1 and 2. The estimate of erroneous de- 
coding is very conservative. Work is continuing to obtain 
a more exact figure (the problem is one of finding the 
distribution of words of weight 8 in the shortened BCH 
code). The statistics given make the proposed modifica- 
tion a worthwhile one, considering the economy and 
simplicity of the modifications. 


J. Modification of NASCOM Error-Defection 
Equipment for Error Correction 

a. Introduction, An error detector will be installed at 
the DSN stations as part of the NASCOM high-speed 
data modem to be supplied to the DSN. It is suggested 


b. Error- correction procedures. The (597,564) short- 
ened code used by the NASCOM is encoded using the 
procedure in Ref. (9), To the information bits flo, fli, ^ 2 , 
* • • Osfis) is assigned the polynomial 

563 

i=.o 
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Table 1. p(n), the probability of making n errors 0 < n < 8 







P(3) 

P(4) 



P(7) 

P(8) 


0.0494 

0.1490 

0.2242 

0.2246 

0.1688 

1 

0.1009 

0.0503 




0.5487 

0.3295 

0.0988 

0.0197 

2.944 X 10-® 

3.514 X 10“* 

3.488 X 10“® 





0.0565 

1.693 X 10-* 

3.374 X 10"* 

5.037 X 10-^ 

6.005 X 10-* 

A 

0 




0.9940 

5.964 X 10-® 

1.786 X 10-* 

3.561 X 10"* 

<io-‘» 

<10-** 

• 

0 

V 




Table 2. Output word error probabilities 


Four modes of operotion 

Input 

bit 

error 

probability 

Probability 
of one 
or more 
errors 

Detecting 7 errors 

Correcting 1 
Detecting 6 

Correcting 2 
Detecting 5 

Correcting 3 
Detecting 4 



A 

B 

c 

A 

B 

C 

A 

B 

C 

A 

B 

C 

5 X 10-^ 
10-^ 

10““ 

10“* 

0.95 

0.45 

0.06 

0.006 

0.946 

0.45 

0.06 

0.006 

0 

0 

0 

0 

<0.004 

<10-* 

<10-*'* 

<10-** 

0.80 

0.12 

0.002 

1.8 X 10-* 

0.15 

0.75 

0.987 

0.994 

<0.007 

<10-* 

<10-** 

<10-** 

0.57 

0.023 

3.4 X 10-* 
3.6 X 10“® 

0.42 

0.95 

1 

1 

<0.02 

<10-* 

<10-*® 

<10-*® 

0.17 
0.003 
5 X 10-^ 
<10-^° 

0.62 

0.992 

1 

1 

<0.04 

<10-* 

<10-* 

<10-** 

A = Probability of detected but uncorrectoble error. 

B = Fraction of words with detected errors that are corrected. 

C = Output word error probobiiity (these are estimates and should be regarded as upper bounds). 


This M(x) is then divided by P{x) = /i fs + 1) where 
h(x)=X^^ + X^+l; f^=x^^^ + X^ + X^+X + h fo = 
The remainder 

32 

R(x) = X b,x* 

t ==0 

is formed, and the parity checks (&i) are sent* (In order 
to have a nonzero modulus after division by F(x), certain 
positions in the encoder are complemented by a poly- 
nomial F(x) at the end of encoding. This is done to facili- 
tate synchronization. We take this complementation into 
account by forming the correct M(x) + irom the 
received word by re-adding in the F(x).) From M + R we 
compute S(x), the remainder of the division of M + R 
byF(x). 

Evaluate S(/?^) S,, / = 1 , 3, 5, where ^ is a root of 

fi{x). Note that — 0 for a correctly re- 
ceived word. Relabel Cj = Cm = 0, / < 459. Then 


1022 
t = 0 

Thus, if t errors are made in positions /i, /* 2 , . . it, we have 

Si= £ /?-i. ,S3= ^ X 

t=l i=l i=l 

Letting y — l/p, we have the symmetric power sums in 
the error positions /s y\ . . y^. The program is de- 
scribed in the next section. 

c. Decoding program 

( 1 ) Compute Si, S3, S 5 

(2) Calculate R 3 = S 3 + Si^ 

(3) If R 3 = 0, a single error has been made in the posi- 
tion Si =y^, ( > 459. (Note; NASCOM furnishes in 
effect the fact that the received word has one or two 
errors; if it does not have errors, the entire proce- 
dure is bypassed.) 
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(4) If Rz = 0, check if Ss = (S3R3/S1) + Sj®, If so, then 
2 errors have been made. These are the roots of the 
quadratic equation 

-|- (TlX (T 2 — 0; <72 — R 3 /S 1 . 

The quadratic can be solved and mechanized using 
technique of Ref. 1. 

(5) If S5 (S3R3/S1) + Si®, then 3 or more errors have 
been made. Solve the cubic + SiX^ 4- 0^3 = 0, 
where <72 = (Si^Sa + S5) R3; <73 = (S1S5 + Si^Ss + 
R3 ^)/Rs- The roots of this cubic, which are the 
error positions, can be found either via a search 
routine or by using the method for solving cubics 
in Ref. 1. 

(6) If the above cubic does not have three roots in the 
field generated over the 2-element field by an 
error has been detected which cannot be corrected, 
i.e., 4 or more errors have occurred. 

Work is under way to find the most efficient method 
of implementing this decoding procedure on SFOF 
computers. 
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E. Digital Devices Development: 
Nine-Channel Autocorrelator 
and SDS 900 Series Computer 
Interface Logic 

E. B. Lutz 


I. Introduction 

The following is a description of the logic used in inter- 
facing the nine-channel autocorrelator and the SDS 930 
Computer at the Venus DSS test site. The main purpose 
of interfacing is to optimize the efficiency of the computer 
in controlling external devices (Ref. 3). 

The basic setup is easily explained by means of Fig. 8, 
which shows a 30-ft cable connecting the autocorrelator 
with the SDS 930 through the POT and PIN cards which 
contain the interface logic. These so-called ‘7‘Cards” are 
standard ZX-11 blank cable plugs which are modified to 
contain integrated circuit logic. 



Fig. 8. 30-ft Interface cable from nine-channel 
autocorrelator to SDS 930 


In controlling peripheral units, the computer transmits 
a 24-bit word in parallel, of which a certain number of 
bits is decoded to determine which of the external de- 
vices is to be selected. Other bits are decoded and used 
to generate specific functions inside the unit being se- 
lected. These commands are generated by the computer 
and transmitted through the POT card while data from 
the external device will be entered into the computer 
through the PIN card. 

2. Functions 

Table 3 shows the various functions and the corre- 
sponding octal codes for the nine-channel autocorrelator. 
Bits 3 through 11 describe the EOM for external units 
(octal 023), and bits 0 through 2 and 18 through 20 spe- 
cify the class to which the nine-channel autocorrelator 
belongs (octal 06). In addition, the unit designation (bits 
21 through 23) distinguishes between units belonging to 
the same class and is set to octal 0 for the nine-channel 
autocorrelator. The various functions such as load, pre- 
cess, etc., are uniquely specified by bits 12 through 17. 
Whenever data is to be sent to the computer (PIN Oper- 
ation), the function code assumes the octal value of 40. 

3. Description 

Figs. 9 and 10 show the logic diagrams of the POT and 
PIN cards, respectively. Both cards utilize Signetics inte- 
grated circuits to a large extent, in combination with 
passive discrete components, such as diodes, resistors, and 
capacitors. Those lines which carry information from the 
computer's C-register are designated by Cq, Ci, C2, etc., 
and are in one-to-one correspondence with bits 0 through 
23 in Table 3. In addition, such lines as Ssc, Qq2, etc., 
perform the necessary control and timing functions which 
are thoroughly explained in the SDS manuals (Ref. 4). 
The signal lines LLl through LLS and PTUE transmit 
data to the correlator, and their functions can be found 
in Table 3. 
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Load D/ A 

Stop count — start output 
Process output 
Restart count 
Shift D/ A, insert zero 
Shift D/ A, insert one 
Lower X-V pen 
Raise X-Y pen 
PIN setup 
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Table 3. 9-Channel autocorrelator function commands 


C-Register 


EOM EOM Function 

1st octal Istoctol 2nd octal 1st octal 2nd octal 


0 1 2 I 3 4 5 6 7 8 9 10 11 12 13 14 I 15 16 17 


Class 


2nd octal 


18 19 20 I 21 22 23 


TO NINE-CHANNEL 
CORRELATOR 


N2NI6.3 1 

+6.5 V 

-TO V 


i\_yj 


5|oa 





'IN754A ; 

HN746A 

► GROUND 

LOGIC EQUATIONS 

X ~ CQ'Ci'Cg* C|Q* C||* C|Q*C|^ • c 

SET FFI*i.r-C,2-Sys 

RESET FFI» Rti +St + .T-Cig-Sys 

oROUND 





Fig. 9. POT 

card logic 
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When an EOM 40 (PIN setup) is executed, FFl is 
set, enabling the gates on the PIN card which then will 
recognize data sent from the correlator through the chan- 
nel lines CHI, CH2, etc. Those lines then enter the C- 
register through Cdl5, Cdl6, etc. The level shifters (LS) 
serve the purpose of translating negative logic from the 
correlator into the appropriate positive levels which are 
standard to the SDS 930 computer. FFl, when set, also' 
signals the computer through Rt that it must exit the 
wait phase and start receiving data. In addition, an inter- 
rupt line (INT 201) is derived from the PIN card logic 
and fed to a corresponding interrupt jack-mounted on a 
card which plugs into the same connector cage. 


various logical conditions. Final tests at the Venus DSS 
were found to be satisfactory, and the logic is currently 
installed and operative. 




U441 

Microwcrv© Maser Development: 
Extension of Klystron 
Tuning Range 


R. Clauss 


In order to insure reliable operation, the interface logic 
circuits are powered by 6.5-v supplies which are derived 
from the standard SDS -f- 8 v. This was done primarily 
to avoid exceeding the 8-v breakdown ratings on the 
integrated circuits. 

4. Conclusion 

Initial tests on the cable plus interface cards were per- 
formed at JPL on an SDS 930 computer, simulating the 


I. Introduction 

Two tunable traveling wave masers have been de- 
livered to Goldstone (February and June, 1966). Descrip- 
tive information and data for these amplifiers can be 
found in SPS 37-37 through 37-40, Vol. Ill under the 
heading ‘‘Microwave Maser Development." The basic 
difference between the new TWM’s and masers previ- 
ously used at the Goldstone Deep Space Stations is tun- 
ability. The new slow-wave comb structure has been 
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/ 


designed and loaded to give a tuning range which covers 
both 2295 and 2388 Me. When the signal frequency is 
changed, the pump klystron frequency must also be 
changed. An external over-coupled cavity is used with 
the pump klystron to extend its electronic tuning range. 
No mechanical adjustments at the maser package are 
necessary when the signal frequency is tuned to either 
2295 or 2388 Me. This frequency change may be accom- 
plished entirely from the maser instrumentation rack. 

2. Dual Frequency Klystron Package 



VA 246E and VA 246R klystrons have been used to 
pump S-band masers for several years. This field experi- 
ence has shown that these tubes have adequate power, 
good stability, and reliability. These klystrons, however, 
do not have sufficient tuning range to cover maser opera- 
tion at both 2295 and 2388 Me. A brief survey of avail- 
able klystrons has not yielded a tube suitable for the 
required tuning range. A change to other pump sources 
(such as a backward wave oscillator) might present new 
stability, reliability, or compatibility problems and has 
not been considered at this time. 

Laboratory tests of a VA 246R klystron with an ex- 
ternal overcoupled cavity have resulted in a split mode 
of operation with frequencies spread much more than 
the normal tuning range. The external cavity is formed 
by the use of irises and a one-half wavelength section of 
guide through the klystron heat sink. Fig. 11 shows the 
assembly sequence for the klystron and external cavity. 
Figs. 12 and 13 show the irises which fit a standard wave- 
guide flange. Capacitive tuning is accomplished with 
a tuning screw in the heat sink. Klystron tuning range 
is determined by the coupling and Q of the external 
cavity. The iris sizes determine the loaded Q and 
coupling between the klystron and external cavity. 



DIMENSIONS IN INCHES MATERIAL: ALUMINUM 

Fig. 12. iris 1 of Fig. 11 



DIMENSIONS IN INCHES 


MATERIAL: ALUMINUM 


Fig. 13. Iris 2 of Fig. 11 


Mechanical tuning adjustments are required after 
assembly of the dual frequency klystron package prior 
to installation on the maser package. A block diagram 
of the klystron package is shown in Fig. 14. Test equip- 
ment to be used during the tuning procedure is shown 
in Fig. 15. By frequency modulating the klystron a 
display of frequency versus power may be viewed on an 
oscilloscope (for initial adjustments). During normal 
operation with the maser the klystron is in a CW mode. 
The klystron frequency is determined by the reflector 
voltage. In Fig. 16 (a graph of frequency versus reflector 
voltage) it is interesting to note that a double-valued 
region exists. In this region (shaded on Fig. 16) the 
klystron frequency is dependent upon the direction of 
reflector voltage change. For example, if a setting of 
325 V is approached from 350 v, the klystron will oscillate 
at 12.872 kMc. If the setting of 325 v is approached from 
300 V, the klystron will oscillate at 12.724 kMc. Even 
in this double-valued region stability has been excellent, 
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FROM 

POWER 

SUPPLY 


TO PUMP 
POWER 
MONITOR 


TOTWM 


Fig. 14. Dual-frequency klystron package block diagram 



Fig. 15. Test equipment block diagram 
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both during laboratory testing and field operation. The 
TWM requires approximately 100 mw of pump power. 
Fig. 17 shows pump power output versus frequency. The 
data was taken during preparation of a dual-frequency 
klystron package which is now in use at the Mars DSS. 


The phenomenon of "frequency pulling” due to a 
cavity external to the klystron is well known and is 
usually avoided. Here it is used to extend the tuning 
range of a fix-tuned klystron. 



12.66 12.70 12.72 12.74 12^6 1236 12.90 12.92 

FREQUENCY. kMc 


Fig. 17. Klystron power versus frequency 


G. Frequency 
Control: Atomic Hydrogen 
Frequency Standard Evaluation 


W. Higa 


A simple technique has been used to evaluate the 
stability characteristics of a traveling wave tube (TWT). 
In addition to providing pertinent information regard- 
ing the TWT as a possible follow-up amplifier to the 
atomic hydrogen maser, the experiment itself is important 
in establishing a method for future evaluation of the 
atomic frequency standard. 

As shown schematically in Fig. 18 a phase bridge is 
used to measure the relative phase shift between two 



Fig. 18. Block diagram for phase stability measurement 
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signal paths with the test device (unknown) in one of 
the paths. A pair of matched diodes was used as syn- 
chronous detectors and provided adequate sensitivity. 
In order to minimize differential phase errors due to 
frequency excursions in the signal generator, it was 
necessary to equalize the electrical lengths of the two 
paths. This was achieved by means of a precision S-band 
delay line in one of the paths. 

The calibrated delay line also provided the means for 
calibrating the sensitivity of the phase bridge. A reso- 
lution of the order of 0.05 deg was readily achieved. 
Fig. 19 shows a typical recording with passive elements 
(coaxial cable) in both paths of the phase bridge. Short- 
term jitter is of the order of ±0.05 deg, and the drift 
over the period of an hour is of the same order of 
magnitude. 

Fig. 20 shows the stability test with the TWT (Watkins- 
Johnson Model WJ 268-2). The short- and long-term 
phase variation is seen to be increased slightly to around 
±0.25 deg. The origin of the short-term jitter needs to 
be investigated; its regularity suggests the possibility 
of inadequate filtering in the built-in power supply and 
may be easily resolved. 

The dual-channel traveling wave maser (TWM) for 
the hydrogen line frequency (1420 MHz) is being com- 
pleted. Stability tests similar to that described here will 
be made before using the TWM to evaluate the short- 
term stability of the hydrogen maser. 



Fig. 19. Calibration test-run for system of Fig. 18 
with passive elements 
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Fig. 20. Test results for TWT (Note reduced sensitivity 
as compared with Fig. 2) 
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H. Improved RF Calibration 
Techniques: Noise Temperature 
Calibrations 

7, Y. Ofoshi and C, 7. Stelzried 

1 . R&D Cone 

Calibrations of the Venus DSS R&D cone were per- 
formed during May and June, 1966 at JPL on the roof of 
Bldg. 238. The calibration work basically consisted of 
measurements of: (1) transmission line insertion losses 
and reflection coefBcients, (2) antenna noise temperature, 
(3) receiver noise temperature, and (4) high and low gas 
tube temperatures. 

Descriptions of the RF cone instrumentation, the RF 
calibration system, and a tabulation of recently measured 
insertion losses and reflection coefficients were presented 
in SPS 37-40, Vol. Ill, pp. 93-100. 

Noise temperature measurements were performed by 
the Y-factor method, which involves the use of two 
thermal reference noise standards and switching between 
the unknown and the thermal noise standards. For the 
cone calibrations, the two thermal standards used were 
a cryogenically cooled waveguide termination and an 
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Fig. 21. R&D cone antenna noise temperature defined at the antenna output reference flange 


ambient waveguide termination. The YTactor equations 
as well as probable errors for this measurement technique 
are discussed in SPS 37-39, VoL IV, pp. 86—91. 

Fig. 21 shows the results of antenna noise temperature 
calibrations made at 2388 Me for the antenna in the 
circular polarization mode. The antenna noise tempera- 
ture in this figure is defined at the antenna output refer- 
ence flange which (for this cone) is the input to the 
circular waveguide rotary joint preceding the turnstile 
junction. Each data point represents the average of 10 
sets of Y-factor measurements. The probable error asso- 
ciated with each data point includes the probable errors 
of the effective cryogenic load temperature, ambient 
load temperature, Y-factor ratios, and the measurement 
precision index. The probable error, however, does not 
include the errors due to mismatch, which are treated 
separately. 

A calibration run was made each hour over approxi- 
mately 24-hr periods to permit study of antenna tempera- 
ture variations due to weather and radio sources. The 
first 24-hr sequence used a liquid nitrogen-cooled wave- 
guide load (MMC SR 8135 SN 004) as the cryogenic 
reference standard, and the second 24-hr sequence used 
a liquid helium-cooled waveguide load. Comparisons of 
the results from the two 24-hr runs provide a check on 


the consistency of the antenna temperature measurements 
as well as a cross check on the stability and accuracy 
of the cryogenic loads. In Fig. 21 it is of interest to note 
the periodic rise in the antenna temperature caused by 
the Sun passing overhead during the noon hours. 

Figs. 22, 23, and 24 show plots of the 2388-Mc. line- 
loss temperature contribution, receiver temperature, and 
gas tube temperatures, respectively. The line-loss tem- 
perature contribution is the difference between antenna 
temperature defined at the receiver reference flange and 
the antenna temperature defined at the antenna output 
reference flange ( previously defined ) . The receiver refer- 
ence flange for this cone is defined to be the output of 
the 26-db directional coupler (SPS 37-39, Vol.III). 

The overall average noise temperatures for the R&D 
cone at 2388 Me are summarized in Table 4. The overall 
averages were computed from the data points shown 
in the figures, with the exception of omitting those an- 
tenna temperature data points obtained when the Sun 
was overhead. The probable error shown with each 
average is the measurement precision index of the overall 
average. The tabulated uncertainties due to mismatch 
errors are the present best estimates and are based upon 
the noise temperature mismatch error analysis for two 
port networks presented in SPS 37-40, Vol. IV, pp. 202-208. 
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Fig, 22. R&D cone noise temperature contributed by antenna line lots via Arm A 



Fig. 23. R&D cone maser gain and the receiver and follow-up noise temperatures 
defined at the receiver reference flange 
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Fig. 24. R&D cone high and low gas tube noise temperatures defined at the receiver reference flange 


Table 4. Summary of average noise temperatures 
for R&D cone at 2388 Me 


Noise source 

Reference 

flange 

Average noise 
temperofure, 

Peak uncer- 
tainties due 
to mismatch. 

Antenna vio 
Arm A 

Antenna output 

5.0 ± 0.05 fpe;« 

±1.0 

Transmission tine 
loss contribu’ 
tion via Arm A 


8.4 ±0.01 {p9)u 


Receiver 

Receiver 

7.2 ± 0.12 (pejM 1 

±0.2 

System via Arm A 

Receiver 

20.6 


Followup receiver 
contribution 

Receiver 

1.6 ± 0.05 rpe^jf 


High gas tube 

Receiver 

42.4 ± 0.05 


Low gas tube 

Receiver 

4.1 ±0.02 ip9)u 


denotes the measurement precision index of the average. 


2. 85-Ff O/sh/R&O Cone 

After calibrations were completed at JPL, the R&D 
cone was shipped to the Venus DSS and installed on the 
85-ft antenna system. The results of two weekly noise 
temperature calibrations of this system at 2388 Me are 
shown tabulated in Table 5. These preliminary system 


temperature measurements of the 85-ft antenna system 
give approximately 22.0® K (with 0,4® K for the followup 
receiver) as compared to a 19.4® K system temperature 
measured for the R&D cone on the ground if we assume 
the same followup contribution for an equal-basis com- 
parison. 

The difference in system temperatures for the R&D 
cone on the ground and on the 85-ft dish should be due 
mainly to the difference in the antenna temperatures. If 
the system temperature difference of 2,6® K is added to 
the antenna temperature of 5.0 ®K measured for the R&D 
cone on the ground, then the antenna temperature meas- 
ured for the R&D cone on the 85-ft dish should be 7.6 ®K. 
However, examination of the preliminary data in Table 5 
reveals that the antenna temperature measured for the 
R&D cone on the 85-ft antenna is approximately 6.6®K. 
The discrepancy of 1®K is believed to be caused by un- 
certainties in the liquid-nitrogen— cooled load used in the 
measurement. Further examination of the data in Table 5 
shows that the receiver temperature measured during 
the weekly calibrations is hotter than those measured 
when the cone was on the ground. A difference between 
the calibrated and actual value for the cryogenic load 
temperature can result in the measured receiver temper- 
ature being too high and the antenna temperature being 
too low. 
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Table 5. Summary of 85-ft dish/R&D cone weekly calibrations at 2388 Me 


Noise source 

Reference flange 

Turnstile 
junction arm 

Noise temperatures^ ‘’K 

Peak uncertainties due 
to mismatch, °K 

GMT day 198 

GMT day 205 



A 

6.6 ± 0.48 (p.e.) 

6.3 ± 0.48 (p.e.) 

±1.2 

Antenna 

Antenna output 







B 

6.8 ± 0.48 (p.e.) 

6.6 ± 0.48 (p.e.) 

±1.1 

Transmission line 


A 

8.6 

8.4 


loss contribution 


B 

8.3 

8.1 


Receiver 

Receiver 


6.8 ± 0.57 (p.e.) 

7.6 ± 0.57 (p.e.) 

±0.3 



A 

21.9- 

22,3 


System 

Receiver 

B 

21.9 

22.2- 


Follow-up receiver 

Receiver 


0.4 



contribution 






High gas tube 

Receiver 


44.8 

44.3 


Low gas tube 

Receiver 


4.4 

4.2 


‘Difference between given system temperature and that obtained by adding contributing temperatures is due to round-off. 


A study is presently being made to determine the 
causes of discrepancies in the cryogenic load tempera- 
ture. One possible source of error being investigated is 
the pressure differential between the atmosphere and 
the gas of the cryogenic fluid inside the dewar. 


Evaluation of the RF characteristics of a number of 
radome materials was made from data received on the 
Venus DSS cone. Radiometric measurements were made 
with and without samples of the radome material over 
the horn (Fig. 25), and four separate computer programs 
were used to reduce and evaluate the data taken. 
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I. RF Techniques: Evaluation of 
Radome Materials for the 
Ascension Island DSS Antenna 

6. Seidel and C. Steizried 


I. Introduction 

RF characteristics of radome materials have been eval- 
uated by the use of radiometric measurements. Computer 
techniques were used for the evaluation and error anal- 
ysis data handling. 

Unfortunately, the time available for this experiment 
did not allow for repeatability data on the various sam- 
ples. Therefore, in some cases the tolerances on the 
results are of the same order of magnitude as the results 
themselves. It is felt, however, that the results of the 
experiment with a probable error of approximately 0.01 
db are still useful. 



Fig. 25. Venus DSS cone with sample of radome 
material over the horn 
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2 . Theory 

A block diagram of the radiometer system used for 
the radome material evaluation is shown in Fig. 26. The 
loss of the material L/j at an ambient temperature To 
is determined from system effective noise-temperature 
measurements looking at the sky ( T'^ )g and (T'sa)i, with 
and without the material placed between the horn pointed 
at zenith and the sky. The effective noise temperatures 
of the sky T' and T 2 with and without the material in 
place, are related by Ref. 5. 

T'-T2=(To-T,)(^1-^^ ( 1 ) 


The noise-temperature measurements are most con- 
veniently referred to the maser input ( Ref. 6, IBM 1620 
computer program CTS/20B) through the transmission 
line loss L 2 

= (3) 

Manipulating Eqs. ( 2 ) and ( 3 ) 

= (Ts^),-(Tsa), 



The probable error of Li? is given by Ref. 7 


Hie measured change in system temperature is therefore 

(T's^)2- ( 2 ) 


♦ 

rt 


RADOME MATERIAL LOSS 


Tsa 




HORN 



and 


PEl 2 = probable error of the transmission line 
loss, I >2 


TRANSMISSON LINE LOSS L2 


Tsa 


LOW NOISE 
PREAMPLIFIER 


Fig. 26. Block diagram of radiometer system for 
radome material evaluation 


PEtq — probable error of the ambient tempera- 
ture, To 

PEt 2 — probable error of the sky temperature, T 2 

PE{t^^ )2 — probable error of the system temperature 
measurement with radome material in 
place due solely to measurement scatter 

PE{t^^)^ = same as except with radome 

material removed 


=r probable error of system temperature 
measurements due to instrument bias 


errors 
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The system temperature bias errors PE are corre- 

lated by the measurements, whereas the errors due to 
system temperature measurement scatter are uncorre- 
lated and add in the usual squared manner. Eqs. (4) 
and (5) have been programmed for the IBM computer 
and designated CTS/27A. 

Fig. 26 can also be used to derive the effect on system 
noise performance of the evaluated radome material for 
any receiving system. The increase in effective system 
temperature defined at the input of the low noise pre- 
amplifier is obtained by manipulating Eq. (4) 


(Ti 


^Tai — ‘ 




( 6 ) 


where 



AT, 


8.1 


^TgA 

To-Tj 

^Tba 


I 


Ft) 


and the probable errors are as previously defined. Eqs. 
(6) and (7) have been programmed for the IBM 1620 
computer and designated CTS/27B. 


The parameters in Eq. (6) pertain to the receiving sys- 
tem under evaluation. The probable error is 




= ^PE 


3AT 


8A 


+ PE 


'^*1 BL , 

' BATsa 


+ {PEr^ + PEl^) 




BU 


IT 


BATs 


BTo 


( 7 ) 


3 . Results 

Table 6 lists the pertinent results of this evaluation. 
Some of the samples were unaltered, some had seams, 
and some were soaked in brine (several days) before 
measurement, as indicated. 

The results of the experiment indicate that, from an 
RF standpoint, Sample V-25306 should be selected for 


Table 6. Results of tests for various radome materials 


Sample 

Ln 

db 

db 

•K 

•K 

AT«ai 

•K 

•K 

Test Cenditiont 

15187 

0.0282 

0.0095 

1.52 

0.511 

1.54 

0.520 

No seom, brine-soaked 

15187 

0.0310 

0.0079 

1.67 

0.425 

1.70 

0.432 

Seam, brine-soaked 

15878 

0.0585 

0.0090 

3.14 

0.482 

3.19 

0.491 

Seam, brine-soaked 

15878 

0.0650 

0.0092 

3.49 

0.493 

3.55 

0.501 

No seam, brine-soaked 

V-25284 

0.0265 

0.0093 

1.43 

0.500 

1.45 

0.509 

Seam, brine-sooked 

V-25306 

0.0162 

0.0090 

0.87 

0.485 

0.89 

0.493 

Seam, brine-soaked 

15878 

0.0072 

0.0088 

0.38 

0.475 

0.39 

0.483 

No seam 

V-25284 

0.0071 

0.0084 

0.38 

0.454 

0.39 

0.461 

No seam 

15187 

0.0079 

0.0086 

0.42 

0.464 

0.43 

0.472 

Seam 

V-25306 

0.0128 

0.0097 

0.69 

0.523 

0.70 

0.523 

No seam 


Lm 

^Tsai 


= measured insertion loss of the sample 
= probable error of the insertion loss measurement 

= expected Increase in the system temperature of the Ascen- 
sion Island DlSS receiving system when the system is at an 
ambient temperature of 25 

AJba 2 — tome os except at an ambient temperature of 30 "C 

= probable error of 
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the radome material. However, since none of the mate- 
rials will vastly degrade the system performance, other 
considerations (i.e,, structural strength) may prove to 
be the deciding criteria on the material. 


■ JV67 14445 

J. cTTicient Antenna Systems: 
Aperture Blockage and 
Surface Tolerance Loss 
Calculations for Non- 
Uniform Illumination 
and Error Distribution 

A. Ludwig 


If an antenna with a reasonably well-defined aperture 
is divided into several segments (concentric rings for ex- 
ample), the efficiency of the whole antenna is related to 
the efficiency of each segment by 


( 1 ) 


r N -12 

= (-nPn^n)^J 

where 

rf = overall efficiency 
= efficiency of the nth segment 
N = number of segments 

Pn = fraction of total radiated power intercepted by 
the nth segment 

On = fractional aperture area in the nth segment 
All efficiencies are defined by 

G G 


v = ^ 


Go 


( 2 ) 


47tA/a2 

where in the case of segments 

G = the gain of the nth segment, measured in the di- 
rection of maximum gain of the whole antenna 
(not necessarily the maximum gain of the seg- 
ment) assuming that the power incident on the 
segment is total radiated power (no spillover) 

A r= physical aperture area of the nth segment 


Eq. ( 1 ) simply expresses the summation of the field con- 
tributions of each segment, assuming that the contribu- 
tions from all segments are in phase and normalized by 
total area and total radiated power. For overall efficiency, 
spillover is accounted for by the fact that the sum of Pn 
is in general less than 1. Eq. ( 1 ) may be easily extended 
for the case of phase differences in the contributions of 
individual segments. 

One application of Eq. (1) lies in the evaluation ot 
aperture blockage loss. Areas of the aperture that are 
shadowed by a subreflector, feed, or support structure, 
may be considered as one or more segments with zero 
efficiency, or with an estimated efficiency in the case of 
partial RF transparency of the obstacle. The effect of a 
single totally blocked segment is found to be 

v' = [ (v)^ - (anPnVn)’^ ]^ = 1 - ' (3) 

where 

— overall efficiency with the nth segment totally 
blocked 

rj = overall efficiency in the absence of blockage 

rfn = efficiency of the nth segment if it were not blocked 

This formulation has been previously used by Gray 
(Ref. 8 ). Typically since illumination is more 

likely to be uniform over a segment, and since spillover is 
included in rj but not in For the special case of iy = iyn> 
and an = Pn (percentage of intercepted energy equals 
percentage of area) Eq, (3) reduces to 

77 ' = 17 [1 — ttn]^ (4) 

which is a form frequently used for aperture blockage 
loss. As an example of the error introduced by using 
Eq. (4) instead of Eq. ( 3 ), in a recent problem subreflec- 
tor blockage represented 1 . 1 % of the aperture area, yield- 
ing a blockage loss factor of 0.978, using Eq. (4). The 
fractional energy in the blocked region was 2.54%, con- 
siderably higher than the fractional area, since the illumi- 
nation was at a maximum in this region. The unblocked 
overall 1 ; was 0.739. Assuming that = 1 (the illumina- 
tion is virtually uniform in the blocked region), Eq. (3) 
yields a blockage loss factor of 0.962 (which agrees well 
with a more sophisticated computer calculated value 
using an integral technique, Ref. 9, similar to that given 
in Ref. 10). In this case the difference in the two compu- 
tations is only 0.07 db, but for careful efficiency calcula- 
tions this is not insignificant. Also, for higher blockage the 
difference is amplified. 
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For 




< < 1 this becomes 


A second possible application of Eq. ( 1 ) is in determin- 
ing surface tolerance loss. When surface errors are not 
uniformly distributed over the aperture, segments over 
which errors are nearly uniformly distributed may be 
evaluated for RMS error and the loss applied to the indi- 
vidual Tfn values. If each segment is individually best-fitted 
( as it should be ) , a modified version of Eq. ( 1 ) , including 
the phase difiFerences of the contributions of each seg- 
ment, should be used. For the special case of all contri- 
butions adding in phase, efficiency would be given by 

-ClI— 

where an is the RMS error of the nth segment.^ It is of 
interest to determine an equivalent overall a that would 
yield the same loss given by Eq. (5). This may be deter- 
mined by 

exp - ^ ^ J " 

r ^ / 47 ro - 2\%“]2 

= 2 (o^nPnqn CXp y J (6) 

Defining a function 



(7) 

It is easily found that 

i /'$> (anPnrin)^ 

( 8 ) 

n = 1 


'The technique of dividing an aperture into segments is independent 
of the method used to calculate losses for each segment; a discus- 
sion of the validity of the form used above when applied to seg- 
ments of various sizes is beyond the scope of this article. However, 
it should be noted that one must be careful that the criteria for the 
validity and accuracy of a given formula are satisfied over each 
segment. 


(9) 

Z (o.nPnVn)'^ 

n = l 

Therefore, for low surface loss, the error should be 
weighted by the square root of segment area, power and 
efficiency. For high loss, Eqs. (8) or (5) should be used 
to determine the equivalent a. A weighting of error by 
area and power has been suggested by Mezger (Refs. 11 
and 12 ) , but he does not give a derivation, and therefore 
the reason for the disagreement with the above derivation 
is not known. It is worth noting that for uniform illumina- 
tion, Tjn — 1.0, a„ = Pn, the Weighting reduces to the 
square root of area squared, or simply area, which is 
again a conventional method of weighting error. 


It has been determined that a beamshaping flange on 
a subreflector for a Cassegrainian antenna can appreciably 
improve antenna efficiency ( Ref. 13 ) . This article is pri- 
marily concerned with the effect of operating at a fre- 
quency significantly above the design frequency of the 
flange. 

2 . Evaluation of Frequency Dependence 

The subreflectors for the 85-ft- and 210-ft-D DSIF 
antennas have flanges designed for S-band operation. To 
evaluate the performance of these subreflectors at several 
frequencies, the subreflector configurations and illumina- 
tion were input to a JPL-developed computer program 
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K. Efficient Antenna Systems: 
Frequency and Size Dependence 
of a Subreflector Beam shaping 
Flange 

A. Ludwig 

1. Introduction 
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Fig. 27. Subreflector configurations (o) 210-ft-D antenna (b) 85-ft-D antenna 


which calculates scattering from an arbitrary surface of 
revolution (Ref. 14). The computed scattered patterns 
were then input to another computer program which 
evaluates antenna efficiency (Ref. 15). Both subreflectors 
were evaluated at S, C, X, and Ku-band. The subreflector 
configurations are shown in Fig. 27. The flanges for both 
subreflectors are 12 in. wide, since it has been determined 
experimentally that the optimum flange is about two 
wavelengths wide, independent of the subreflector size 
(for large subreflector sizes in wavelengths). The sub- 
reflector illumination, which was assumed to arise from a 
dual mode feedhorn (Ref. 16) scaled for each frequency 
(such that the pattern would remain constant) is shown 
in Fig. 28. Several of the resulting computed scattered 
patterns are shown in Fig. 29. It is seen from Fig. 29 that 
operating above the design frequency causes substantial 
distortion at the edges of the patterns. The phase patterns, 
which are not shown, are also perturbed in this region. 
Scattered patterns from the subreflectors with the flange 
removed (pure hyperboloids) are shown in Fig. 30. 


The resulting antenna efficiencies for the illumination 
patterns shown in Figs. 29 and 30, and for the other cases 
evaluated, are shown in Fig. 31. Efficiency includes losses 
due to spillover past subreflector and main reflector, am- 
plitude and phase illumination losses, and cross- 
polarization losses. It is seen that the flange results in 
substantial improvement in performance at the design 
frequency, and has very little effect on efficiency, either 
positively or negatively, at the higher frequencies. How- 
ever, it should be noted that the subreflectors without 
flanges are smaller than optimum, and the efficiencies for 
these cases would be 1 to 2% higher if the hyperboloids 
were somewhat larger. 

The most interesting fact illustrated by Fig. 29 is that 
the flange is much more effective on the smaller subreflec- 
tor. The illumination efficiency of the smaller antenna is 
actually greater than that of the larger antenna — a com- 
plete reversal of the normal situation as illustrated by the 
no-flange cases. 
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Fig. 29. Scattered patterns of subreflectors with flange (a) 210-ft-D antenna S-Band, 
(b) 210-ft-D antenna X-Band (c) 85-ft-D antenna S-Band (d) 85-ft-D antenna X-Band 
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This curious situation is the result of the fact that, on 
the smaller subreflector, the flange subtends a larger por- 
tion of the feed pattern, resulting in 3% less forward 
spillover. This increase more than makes up for the larger 
diffraction effects of the smaller subreflector. Making the 
flange larger for the large subreflector would result in 
distorted patterns similar to the higher frequency pat- 
terns, and making the hyperboloidal section larger would 
cause more loss due to increased illumination taper than 
could be gained by decreased spillover. Thus, it appears 
that for the hyperboloid/flange combination, smaller 
subreflectors are better than larger subreflectors. The dif- 
ference in efficiency is small, but it indicates an area of 
potential improvement in the design of large subreflectors. 

Hb7 Ikkk? 

“t.^AdvancecTDevelopment of 
Microwave Antenna Subsystems: 
100-kw SCM Cone Tests 

R. W. Hartop 

As part of the continuing effort to maintain state-of-the- 
art capability in the microwave antenna subsystem of the 
DSIF, advanced development has begun on the design 
of waveguide components capable of handling from 500- 
to 2000-kw CW transmitter power. SPS 37-40, Vol. Ill, 
p. 80, described some special tests which have been per- 
formed on present Goldstone Duplicate Standard (GSDS) 
components in order to determine their limitations under 
very high power. 

During this reporting period, these tests were extended 
with the testing of a complete GSDS S-band Cassegrain 
monopulse ( SCM ) cone assembly at power levels up to 
100 kw CW. These tests were performed at the Venus 
DSS, using an available 100-kw klystron. The SCM cone 
was mounted on one of the elevators so that the RF power 
radiated through the roof port, and an R&D traveling 
wave maser was installed in the cone with complete in- 
strumentation to monitor system noise temperature and 
to detect noise bursts due to corona discharges or con- 
taminants. In addition, a 12-channel recorder was used 
to monitor thermal temperatures of critical components 
within the cone by means of calibrated thermocouples. 

Two complete power runs were made at various power 
levels from 20 through 100 kw. The first test used the RCP 
bridge circuit of the SCM feed and was completely suc- 
cessful. After shutting down and switching to left-handed 
circular polarization (LCP), 100 kw was again trans- 


mitted through the feed, but after 3 min at full power 
( preceded by 5 min at 75 kw ) a sudden noise burst was 
observed on the spectrum analyzer and the chart recorder. 
The transmitter was manually turned off, and the power 
was then gradually increased in 10- to 15-kw increments 
as in previous runs. This time, however, noise bursts 
would occur every time at about 40 kw, and the tests 
were concluded at this point. 

It had been anticipated that one of the components that 
might give trouble during very high-power testing was 
the waveguide load that is used to terminate several 
unused ports in the system (such as the cross-difference 
ports of the SCM feed). An attempt was made during the 
tests to monitor the temperature of these loads to prevent 
them from reaching their 350 ®F melting point. However, 
because of the low thermal conductivity of the load 
material, combined with the high conductivity of the 
aluminum waveguide to which the thermocouples were 
attached, it was difficult to estimate what the actual load 
material temperature was at a given time. 

Thus, when the noise burst problem arose during the 
LCP test, it was assumed that one of the cross-difference 
loads had overheated and introduced contaminants into 
the system. However, later investigation showed this 
assumption to be incorrect; a load on the receiver arm of 
the diplexer had slightly overheated and given off enough 
vapor to cause the noise bursts. Fortunately, no contami- 
nants were introduced permanently into the system, and, 
after replacing the load, the SCM cone passed a standard 
diplexing test at the Echo DSS antenna test range. 

The fact that a load which is theoretically isolated by 
90 db from the transmitter signal could overheat has 
prompted several further investigations. In one test, sim- 
ilar loads were power- tested to determine their capacity. 
Briefly, it was found that the load could tolerate about 
30 w indefinitely with a load temperature reaching 215® F. 
After 24 min at 50 w, however, the load material began to 
smoke. This suggests that the load within the cone must 
have been exposed to a power level in excess of 50 w, 
which effectively proves that the signal must have been 
a harmonic of 2110 MHz in order to pass through the 
2200 MHz cutoff sections of the diplexer. Allowing for 
poor impedance match of the waveguide components at 
the harmonic frequencies, it appears that the klystron was 
radiating a total harmonic power of at least 100 w. Fur- 
ther tests will be performed to substantiate these findings. 

An anomalous result of the tests was the peculiar be- 
havior of the system noise temperature during diplexing 
at high power levels. The phenomenon is being investi- 
gated in the laboratory as maser availability permits and 
will be discussed further in the next reporting period. 
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Despite the difiBculty with hannonic power, the 100-kw 
SCM cone tests were highly successful and indicate that 


the design goal of 2000-kw CW microwave components 
may well be achieved in the relatively near future. 


References 

i 

I 

1 . Berlekamp, E., Rumsey, H., Solomon, G., ''Solutions of Algebraic Equations Over 

I Fields of Characteristic 2/' SPS 37-39, Vol. IV, Jet Propulsion Laboratory, Pasa- 

dena, California, June 30, 1966, pp. 219—226. 

2. Peterson, W. W., "Error-Correcting Codes," John Wiley & Sons, New York, N.Y. 

1961. 

3. Martin, W. L., "Digital Communication Tracking: Interfacing with SDS 900 Series 
Computers," SPS No. 37-39, Vol. Ill, Jet Propulsion Laboratory, Pasadena, Cali- 
fornia, May 31, 1966, p. 52. 

4. "SDS 925/930/9300 Computers INPUT-OUTPUT," Scientific Data Systems, Santa 
Monica, California, January 1 965. 

5. Steizried, C. T., "Temperature Calibration of Microwave Thermal Noise Sources," 
IEEE Transactions on Microwave Theory and Techniques, Vol. MTT-13, January 
1965, p. 129. 

6. Steizried, C. T., Otoshi, T. Y., "R&D Listening Cone — Advanced Antenna System 
Daily System Noise Temperature Calibrations," SPS 37-39, Vol. Ill, pp. 86—103. 

7. Worthing, A. G., Gaffner, J., "Treatment of Experimental Data," John Wiley and 
Sons, Inc., New York, N.Y., 9th Printing, 1 960, p. 209. 

8. Gray, C. L., "Estimating the Effect of Feed Support Member Blocking on Antenna 
Gain and Sidelobe Level," Microwave Journal, pp. 88—91, March 1964. 

9. Ludwig, A., et al., "Computer Programs for Antenna Feed System Design and 
Analysis," Sect. IX, Technical Report 32-979, Jet Propulsion Laboratory, Pasa- 
dena, California (to be published). 

10. Wested, J. H., "Shadow and Diffraction Effect of Spars in a Cassegrainian Sys- 
tem," p. 2118, Microwave Laboratory, Danish Academy of Technical Sciences, 
Copenhagen, Denmark, March 1966. 

1 1 . Mezger, P. G., "An Experimental Check of Antenna Toleronce Theory Using the 
NRAO 85-Foot and 300-Foot Telescopes," 1964 PTGAP Internotional Symposium 
Program and Digest, pp. 181—185, September 1964. 

12. Mezger, P. G., "An Experimental Check of Antenno Tolerance Theory Using Radio 
Astronomical Techniques," NRAO Internal Report, April 1965. 

1 3. Potter, P. D., "A Simple Beamshaping Device for Cassegrainian Antennas," Tech- 
nical Report No. 32-214, Jet Propulsion Laboratory, Pasadena, California, January 

1962. 

1 4. Rusch, W. V. T., "Scattering of a Spherical Wave by an Arbitrary Truncated Sur- 
face of Revolution," Technical Report 32-434, Jet Propulsion Laboratory, Pasa- 
dena, California, May 1963. 


L 


97 



JPL SPACE PROGRAMS SUMMARY NO. 37-41, VOL. Ill 


References (Confd) 

15. Ludwig, A., "Antenna Feed Efficiency," SPS 37-26, Vol. Ill, April 1964, pp. 200- 
208. 

1 6. Potter, P. D., "A New Horn Antenna With Suppressed Sidelobes and Equal Beam- 
widths," Microwave Journal, June 1963, pp. 71— 78. 


98 


JPL SPACE PROGRAMS SUMMARY NO. 37-41, VOL. HI 


IV. Communications Development Engineering 
N67 14448 


A. Step Recovery Diode Multiplier 

E, R, Caro 

Development of the step-recovery diode has made 
possible efficient high-order frequency multiplication 
without the attendant circuit complexities of Varactor 
frequency multipliers. 

A frequency multiplier has been developed for use in 
the planned Block IIIC modification of the DSIF track- 
ing system local oscillator and transmitter exciter multi- 
plier chains. In order to achieve the modification with a 
minimum amount of change, the physical design of the 
unit was made compatible with the available mounting 
space in the system. The existing multiplier chains in the 
system include a X32 module. In the modified chain, 
this multiplication is accomplished in two steps: a solid 
state X4 and amplifier module, and a step-recovery 
diode X8 unit. A multiplication factor of 8 was chosen 
in order to achieve the power output required consistent 


with the input power capability of presently available 
diodes. This also places the input frequency in the VHF 
region of 260 to 280 Me where sufficiently high levels of 
power can be generated with transistor circuits. An added 
advantage of a X8 unit is the fact that the nearest side- 
band at the output is 260 to 280 Me away from the center 
of the output bandpass. This allows sufficient suppression 
of undesired sidebands to tolerable levels with a mini- 
mum number of output filtering stages. Considerable 
effort was exerted towards making the design as simple 
as possible without compromising performance so that 
the unit would be easily reproducible and require a 
minimum amount of testing time. 

Fig. 1 shows a photograph of the multiplier, and Fig. 2 
shows a schematic diagram. The multiplier consists of a 
step-recovery diode mounted in a two-cavity output 
bandpass filter. The diode is series-fed through a low- 
pass filter FL-1 which isolates the input network from 
the output network. LI, Cl, C2 and C3 form a matching 
section to transform the diode impedance to 50 ohms. R1 
is a bias resistor. 
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Fig. 1. Photograph of step-recovery diode X 8 
frequency multiplier 



J2 

OUTPUT 


Fig. 2. Schematic diagram of X 8 frequency multiplier 


The low-pass filter FL-1 is a 3-element Chebyshev 
filter designed for a 350-Mc cutoff. The filter uses dis- 
tributed coaxial capacitive elements and a helical center 
conductor for the inductive element. 

The diode taps directly into the center conductor of 
the first output filter cavity. The tapping point into the 
cavity was determined experimentally for optimum match 
between the diode and the output load. Tapping directly 
into the center conductor also serves as a convenient heat 
sink for the diode. 


It was discovered during temperature testing (0 to 
50° C) that variations in the length of the center conduc- 
tors due to a change in temperature were epough to 
detune the output cavities and cause a shift in the band- 
pass response of the unit. To correct the condition, the 
adjustable portion of the center conductors was made 
out of Invar. 

Fig. 3 shows the bandpass response of the unit at three 
different temperatures, 0, 25 and 50° C, after incorporat- 
ing the Invar center conductors. As can be seen, very 
little detuning ( < zb 0.25 Me) occurred over this tem- 
perature range. 



I I I I I 

2090 tllO 2130 


FREQUENCY, MHz 

Fig. 3. Bandpass response of X 8 frequency multiplier 
at 0, 25, and 50^C 


Performance of the unit is presented in Table 1. Fig. 4 
shows input versus output curves for four typical diodes 
tested. Note that the output increases with the input 
until saturation occurs at about 4 W of drive power. 


The two output cavities are coupled through an iris 
in the common wall between them. Output coupling is 
accomplished with a loop. Alignment of the unit is 
straightforward. An input signal is applied to the unit, 
and the tuning controls Cl, C2 and L2, L3 are adjusted 
for maximum output. Since a bandpass symmetry require- 
ment has been imposed on the unit, optimum alignment 
of the output cavities is more easily accomplished with a 
swept input signal. Tuning of the output cavities is ac- 
complished by varying the length of the center conductors. 


Tests made on 20 diodes of the same type purchased 
in three different batches indicate that the uniformity of 
the diodes is very good. Very little retuning was required 
after mounting each diode into the module. For the 
exciter X8, where 500 MW of output power is required, 
with an input drive level of 4 W, the yield of usable 
diodes was 75%. However, for the lower level X8 used 
in the local oscillator, where only 200 MW of output is 
required, all of the diodes tested were acceptable. 
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Table 1. Performance of X 8 frequency multiplier" 



Exciter X8 
fo = 2110 Me 

LOXB 

f, = 223S Me 

Input power, dBm 

-1-36 

+33 

Output power, dBm 

+27 

+24 

Bandwidth (-1 dB), Me 

40 

40 

Input impedance, ^ 

50 

50 

Spurious signal level 



below desired signal 



6th harmonic, dB 

>40 

>40 

7th harmonic, dB 

>40 

>40 

9th harmonic, dB 

>40 

>40 

1 0th harmonic, dB 

>40 

>40 

o Input dynamic rang*, s«e Fig. 4. 
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B. MSFN Transmitters - 

R, L Leu 

Implementation of the DSN/MSFN dual transmitter 
was started at the Pioneer Station. Testing of the first 
20 kW transmitter was interrupted by the Pioneer B 
tracking commitment after approximately 90% of the 
tests were completed. The tests will be completed on a 
non-interference basis with the tracking commitments. 
With the exception of incidental amplitude modulation 
(lAM), the transmitter meets or exceeds the design speci- 
fication for all tests which have been accomplished. The 
lAM tests were inconclusive. 



Fig. 4. Input versus output curves of X 8 frequency 
multiplier for four different diodes 


The tuning range of the multiplier was made suffi- 
ciently wide (1900 to 2300 Me) to allow the use of the 
same unit for both exciter and local oscillator applica- 
tion by applying the appropriate input signal and tuning 
the output cavities to the corresponding output signal. 

Conversion efficiency of the X8 multiplier averages 
12.5%, with some diodes exhibiting as high as 25% effi- 
ciency at below saturation levels. 

This same basic unit was also operated as a X 16 mul- 
tipher (140 Me to 2240 Me) by redesigning the low-pass 
filter for a lower cutoff frequency. Performance of the 
unit as a X 16 was quite satisfactory with a corresponding 
decrease in efficiency because of the higher multiplication 
factor. 


The measurement of lAM showed the Klystron com- 
plied with the specification (— 60 dB below the carrier) 
from 2 kHz on up. From 200 Hz to 2 kHz the test equip- 
ment had a dynamic range of only 53 dB. A better test 
setup will be needed to measure — 60 dB below the 
carrier from 200 Hz to 2 kHz. 

The ripple on the output of the high-voltage power 
supply (20 kV, 3 A) was measiured at 25 V peak-to-peak. 
Using the measured value, the incidental phase modula- 
tion (IPM) contributed by the Klystron was calculated at 
1.5 deg, peak-to-peak. The measured value, using the 
test setup of Fig. 5, was 1.2 deg, peak-to-peak. The mea- 
sured results are shown in Fig. 6. 



Fig. 5. Incidentol phase modulation test setup 


With the exception of the heat exchanger, all equip- 
ment for the second transmitter has been delivered and 
is being installed. The second transmitter is to be tested 
on a non-interference basis with the Pioneer B mission. 
The installation and testing of the dual transmitter con- 
figuration is scheduled for completion in September 1966. 
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^ ~ • Fig. 6. Incidental phase modulation across 5K70SG Klystron 

W67 14450 

■ - -C. Krystrort- Protective Device 

£. J. Finnegan 


I. Introduction 

The mercury pool Ignitron crowbar reported on in 
SPS 37-40, Vol. Ill, was tested during this reporting 


2. Ignitron Crowbar 

A mercury pool Ignitron is a cold cathode tube. Mer- 
cury is vaporized by a high current pulse applied to the 


S| 



*Rt IS RESISTANCE OF KLYSTRON HV 

WHEN IT ARCS (APPROX 0-1 II) 

Fig. 7. Test setup for measuring crowbar response 
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ignitor immersed in the mercury pool, thus making it a 
■ mercury gas tube. This type of device is ideal for use as 

I a crowbar due to its capability of handling very high 

I current pulses (thousands of amperes) for short durations 

without any damage. 

The circuit used to test the Ignitron is shown in Fig. 7. 
In the first test, the Ignitron was connected in parallel 
with the spark gap crowbar to test its effectiveness. The 
spark gap alone gives a loud sharp discharge. With the 
\ Ignitron operating in parallel with it, one sees only a 

. very faint arc, thus indicating the efiEectiveness of the 

j Ignitron. 

Fig. 8(a) shows the current discharge wave form as 
viewed when discharging the spark gap only. The spark 
gap takes approximately 3000 fxs to dissipate the energy 
from the 1 /^F capacitor in the power supply. 

Fig. 8(b) shows the current wave form as viewed when 
the Ignitron is used in place of the spark gap crowbar. 



1000 M«/cm 
(a) SPARK GAP CROWBAR 



10 ^«/cm 


(b) IGNITRON CROWBAR 

Fig. 8. Spark gap crowbar and Ignitron crowbar 


By comparing Fig. 8(a) with Fig. 8(b), it can be seen that 
the surge current of the Ignitron is approximately 100 
times greater and the discharge time is approximately 
100 times faster than the spark gap. The Ignitron clears 
the line of dangerously high currents considerably faster 
than the spark gap. This should give greater protection 
to the Klystron, since the high current is shunted from 
the arcing Klystron at a faster rate. The energy that 
reaches the tube may be expressed by 



(see Fig. 7); Rt is approximately 0.1 O when the Klystron 
arcs. The energy reaching the tube with a crowbar and 
operating at 20 kV as it was in the test would be about 
1.4 Joules; this would be adequate protection for the 
Klystron. 

Further testing and modifications are necessary before 
final conclusions may be drawn; however, the test data 
obtained seems to indicate that there is an advantage to 
the use of an Ignitron tube as a protective device. 

N67 14451 

D. Advanced Transmitter 
Development 

R. E. Arnold 

1. Introduction 

During this reporting period, tests were conducted on 
the 1-MW high-voltage power supply at the Goldstone 
Venus Station, and a linear exciter was implemented and 
tested for amplitude-modulating the Mariner 100-kW 
transmitters. This report summarizes the results of the 
tests. 

2. The 1 -MW High-Voltage Power Supply Tests 

Tests were performed on the 1-MW power supply 
installed at the Goldstone Venus Station in order to sub- 
stantiate the conclusions resulting from previous tests 
described in SPS 37-38, Vol. Ill, pp. 56-61. The conclu- 
sion from those tests was that the high-voltage ripple 
and noise could be reduced by replacing the Amplidyne 
exciting the field of the 400-Hz generator with a faster 
responding programmable solid-state power supply. The 
substantiating tests consisted of first breadboarding a 
simple 400-Hz, 3-phase full-wave bridge silicon control 
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Fig. 9. Simplified SCR system test diagram 


rectifier (SCR) power supply at JPL and then tempo- 
rarily installing it in the 1-MW power supply system at 
Goldstone. Fig. 9 shows the test setup in which the SCR 
power supply replaced the Amplidyne and a solid-state 
diflFerential amplifier replaced the vacuum tube regulator 
amplifier normally used in the system. 

Fig. 10(a) depicts a recording of the high-voltage ripple 
during normal operation at Goldstone before the tem- 
porary test setup was installed. The test was run at a 
nominal 28-A generator field current, producing a Kly- 
stron beam voltage of 30 kV at 7.4 A. 

Fig. 10(b) is a recording of the high-voltage ripple 
using the SCR power supply and differential amplifier. 
The test was run at a nominal 23-A generator field cur- 
rent (the near-maximum output of the breadboard SCR 
power supply) for a Klystron beam voltage of 26.5 kV at 
6.4 A. The low-frequency noise and ripple has been al- 
most entirely eliminated. The recording also shows a 
marked increase in 60-Hz ripple, which is thought to be 


due to the instrumentation. The test did fully demon- 
strate that the beam voltage's low-frequency ripple char- 
acteristics can be significantly improved by replacing the 
Amplidyne with a faster responding SCR power supply. 

Because of the encouraging results of the above tests 
and those described in SPS 37-35, Vol. Ill, pp. 56-61, a 
specification for a programmable 175- V, 100-A solid-state 
power supply was written. Subsequently, a contract was 
awarded to build two units; one unit to excite the 1200 
kW 400-Hz generator, and the other unit to excite the 
1750 HP drive motor. They are currently being fabri- 
cated by a vendor, and delivery of the first unit is ex- 
pected in September 1966. 

3. Mariner Transmitter Noise Modulation 

On 8/5/66, as part of the regular once-a-month 100-kW 
Mariner transmitter support, the 100-kW transmitter was 
noise amplitude modulated during radio contact with 
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Fig. 11. Mariner noise modulation experiment, block diagram 


the Mariner IV spacecraft. The intent of the experiment 
was to determine the ability of the spacecraft s trans- 
ponder to return the noise; no commands were sent. 

The amplitude modulation was obtained by modulat- 
ing a 1-W traveling wave tube amplifier (TWT) placed 
between the 3-W buflFer amplifier of the Mariner exciter 
multiplier chain and the Klystron transmitter. The DC- 
coupled modulation input jack, a feature of the particular 
TWT used, was employed to obtain the required modu- 
lation. The block diagram of the experimental setup is 
shown in Fig. 11. The combined Klystron/TWT output 
power characteristics versus TWT modulation input volt- 
age are shown in Fig. 12. The non-linearities at low output 
powers are due to the TWT modulation characteristics. 
The beam voltage of the Klystron was raised slightly 
during the experiment to obtain peak powers of 100 kW. 

The low-frequency modulating signal was obtained 
from a 0- to 35-Hz white-noise generator located in the 



TWT MODULATION VOLTAGE, V 


Fig. 12. Transmitter output power versus TWT 
modulation voltage 
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30 20 10 0 -10 -20 -30 -40 

RELATIVE POWER 

Fig. 13. Transmitter total spectrum 


Venus Station control room. The experiment require- 
ments necessitated that the noise bandwidth be reduced 
by the filter shown in Fig. 11. 

Fig. 13 shows the transmitter total spectrum as re- 
corded by the Venus Station s autocorrelator. The TWT 
was biased for 50-kW transmitter average power, the 
filter RC time constant was 0.5 s, and the amplitude 
modulation was 100%. 

The next Mariner transmitter support will require sup- 
pressed-carrier double-sideband transmitter modulation. 
Preparations are currently being made to provide such a 
modulator. 


N67 14452 

E. S-Band Planetary Radar 
Receiver Development 

C. F. Foster 

1. Introduction 

During this reporting period, the construction and in- 
stallation of the 2388-MHz bistatic radar receiver was 
completed. This receiver will be used for the Mercury 
and Venus planetary radar experiments. 

2. System Description and Implementation 

The 2388-MHz bistatic radar receiver (see Fig. 14) 
consists of several subsystems: (1) The 2388-MHz to 


30-MHz converter, or first converter, (2) the 30-MHz to 
455-kHz converter, or second converter, (3) the 2388- 
MHz test transmitter, (4) the Mars Station programmed 
oscillator, described in SPS 37-38, Vol. Ill, pp. 54-63, 
(5) the Mars-to-Venus Station microwave link, described 
in SPS 37-39, Vol. Ill, pp. 106-108, (6) a portion of the 
Mod IV receiver, described in SPS 37-38, Vol. Ill, pp. 
66-71, and (7) the Venus Station correlator, described in 
SPS 37-25, Vol. Ill, pp. 33-38. The first converter (see 
Fig. 15) is contained in a standard antenna box mounted 
in the 3A section of the Mars Station cone (see Fig. 16). 
It has the following specifications: Input and output im- 
pedance is 50 n, the single-sideband noise figure is 9.25 
dB, the conversion gain is +59.6 dB, the mixer balance 
is 21 dB, and the 1-dB bandwidth is 10 MHz. The second 
converter is also contained in a standard antenna box 
mounted in the 3A section of the Mars Station cone (see 
Fig. 17). It has the following specifications: Input im- 
pedance is 50 Q, output impedance is 75 Q, the conversion 
gain is selected by means of a manual gain control sys- 
tem (SPS 37-40, Vol. Ill), and has a maximum gain of 
+ 55 dB. The converter bandwidth is 400 kHz. This band- 
width was established by the use of a 30-MHz center 
frequency crystal filter having a 400-kHz 3-dB band- 
width. The 30-MHz signal is mixed down to 4554cHz and 
amplified by a video amplifier. A plot of the final system 
bandwidth is shown in Fig. 18. The output of the second 
converter is sent via the microwave link to the Venus 
Station Mod IV receiver where it is divided between the 
455-kHz doppler loop (SPS 37-38, Vol. Ill, pp. 54r-63) and 
the AM spectrum channel. The 400-kHz bandwidth is 
reduced to 27.3 kHz in the Mod IV AM channel for the 
Mercury experiments. The output of the AM channel is 
then sent to the correlator for signal processing. 
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Fig. 14. Bistatic radar block diagram 
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Fig. 15. 2388-MHz to 30-MHz converter block diagram 
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Fig. 16. Photo of 2388-MHz to 30-MHz converter 
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Fig. 18. Bistatic radar converter passband plot 



Fig. 17. Photo of 30-MHz to 455-kHz converter 


Monitoring of the Mars Station converter condition 
(local oscillator drive and power supply voltage) is ac- 
complished by using preset high- and low-level detectors 
adjusted to sense nominal system levels. Their outputs 
are transmitted to the Venus Station via intersite phone 
lines and displayed at the Mod IV receiver (see Fig. 19). 

The programmed oscillator, the maser, and the refer- 
ence frequency generators at the Mars Station have to be 
shared by the 2388-MHz bistatic radar receiver and the 
Mars Station 2295-MHz experimental receiver. Because 



Fig. 19. 2388-MHz bistatic subsystem monitoring 
diagram 


of this fact, and the remoteness of the various subsys- 
tems at the Mars Station, the several cable and frequency 
changes necessary when changing from one system to the 
other are accomplished by coax switches operated from 
a single front panel selector on the Mars Station 2295- 
MHz experimental receiver (see Fig. 20). This provides 
instant change-over from one receiver to the other. 

The signal of the 2388-MHz coherent test transmitter, 
which in past systems had been injected into the maser 
by way of the maser instrumentation box, was provided 
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Fig. 20. 2388/2295-MHz receiver change-over switch 


The 2388-MHz bistatic radar receiver installation is 
successful and is being used to support the Mercury 
radar experiment. 

N67 l'i',53 

F. ResisTive'libop Directional 
Coupler Calibration 

R. C. Chernoff 


7 . Introduction 

The resistive loop type directional coupler presently 
used as a reflected power monitor on high-power trans- 
mitters was calibrated by two different methods. The 
results agree to within experimental error. The frequency 
variation of directivity was also measured and was found 
to be small (1.5%) but adequate for the application. 


with another injection point (see Fig. 21). This was done 
because, in the past, whenever a noise tube or a coax 
switch was changed in the maser instrumentation box, 
the path loss between the maser and the test transmitter 
required recalibration. The 2388-MHz and the 2295-MHz 
test transmitters are both located in the alidade room to 
allow receiver calibrations to be made by a single oper- 
ator. The test transmitter was calibrated using the CW 
power method (SPS 37-38, Vol. Ill, pp. 41-51). 


2 , General 

Resistive loop waveguide directional couplers are used 
as reflected power monitors on all JPL high-power 
S-band transmitters currently in DSN or planetary radar 
service. Although this type of coupler is generally inferior 
to the multi-hole type in directivity and bandwidth, it 
does have the advantage of compactness. This type of 
coupler consists of a loop which protrudes into the wave- 
guide through the E-wall and which is terminated by 
coaxial lines on each end (see Fig. 22). One of the lines 
is terminated in a detector, while the other is terminated 


DIREXmONAL COUPLER 



Fig. 21. 2388/2295-MHz test transmitter maser 


in a tuner and a coaxial load. Directivity in this device 
results essentially from the “right hand rule” for the 



input path 


Fig. 22. Resistive loop waveguide directional coupler 
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propagation of electromagnetic waves. According to this 
rule, the orientation of the magnetic vector with respect 
to the electric vector in the reflected wave is opposite to 
that in the incident wave. The device simultaneously 
couples to both the magnetic field (via the loop) and 
the electric field (via the coaxial termination); and by 
judicious adjustment of the magnetic coupling, high 
directivity can be achieved. (See Refs. 1 and 2 for quan- 
titative discussions of the device.) 

For reflected power monitoring in high-power trans- 
mitters, the directivity required will depend upon the 
maximum allowable reflected power. This is typically 
about 1% ( — 20 dB) of the forward (incident) power, so 
that the required directivity is at least 30 dB if the full- 
scale error in indicated reflected power is to be held to 
less than 10%. The resistive loop coupler is capable of 
this directivity, although over quite small bandwidths 
(see Fig. 23), if carefully adjusted. If, as is usually the 
case, the transmitter is required to operate only in a 
narrow frequency band, the frequency sensitivity is not 
a serious handicap. However, the initial adjustment is 
tedious, and the device, at least in its present form, is 
easily put out of adjustment by even moderately careless 
handling. Moreover, there is some question as to the 
temperature stability of the device. Temperature insta- 
bility, if appreciable, would seriously affect the perform- 
ance at high CW power levels. 

3. Calibration Procedures 

A resistive loop coupler in use as a reflected power 
monitor in the 100-kW transmitter test bed was recently 
removed and recalibrated in an attempt to improve its 



Fig. 23. Directivity of resistive loop type directional 
coupler 


somewhat erratic performance. The coupler was adjusted 
for approximately 35 dB directivity at 2113 Me. Direc- 
tivity was optimized first by rotating the loop and then 
by adjusting the double stub tuner. The purpose of the 
tuner is not, as might be supposed by analogy with 
directional couplers of the coupled transmission line type, 
to match out the residual VSWR of the termination. 
Its primary purpose in this device is to cancel the self- 
inductance of the loop. Directivity was measured by two 
methods. In the first, a sliding short terminated the wave- 
guide system as shown in Fig. 24. Assuming that the 
coupler is inserted in the waveguide system in the 
‘normal” manner, i.e., as if the sliding short were the 
antenna (or dummy load) and the slotted section were 
the Klystron, the RF voltage at the detector (SWR meter) 
would be 

\En\ = \DCE^ + CExe^^^l 

= |CE,j • |D + (1) 

where Ei — amplitude of the incident wave, 

C = voltage coupling coefficient, 

= 27t/Aj, 

I = distance from coupling loop to short, 

D = directivity, 

Xg = guide wavelength. 



Fig. 24. Setup for directivity measurements 
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We assume, because the correction is small, that jp| = 1 
for the sliding short. From Eq. (1) we see that if D = 0 
(infinite directivity), |£d] = |C£i|. As the short is moved 
toward or away from the coupler, the phase of changes 

while the phases of D, C, and £i do not. Therefore, the 
resultant, Ed, is the vector sum of a rotating and a con- 
stant vector. Optimizing the directivity, therefore, consists 
of attaining the smallest possible fluctuation of Ed as the 
short is moved in and out. Moreover, from the observed 
fluctuation in Ed we can calculate D in the following 
manner: First, we assume that our detector is square-law, 
as would be the case for the crystal detector used in 
Fig. 24. The power absorbed by the detector (which we 
assume to be well-matched) is 

Pd = I Ed 

where is the detector impedance. 

Thus, Max Pn (Max | En |)“Gl 

- [[CE, |•(1 + |D|)^G,] 
and Min Pd = (Min | Ed |)^G/. 

= [|C£,1(1-|D1)=G,] 

The “wiggle,” W, observed on the dB scale of the square- 
law detector, is 

= 20 log (1^) (2) 

Solving for the directivity in dB, 

/ 1 01^/20 _ 1 \ 

Directivity = 20 log D = 20 log (3) 

The above derivation assumed that Ei remained con- 
stant as the short is moved in and out. This will be true 
only if the system of Fig. 24 is perfectly matched looking 
from the coupler toward the generator. The waveguide 
attenuator serves as a pad to isolate the system from the 
coax-to-waveguide transition mismatch. The slotted sec- 
tion is used to check the constancy of | Ei | as the short is 
moved. While this arrangement worked reasonably well, 
the match could have been improved by substituting a 


waveguide tuner for the slotted section, which would in 
turn be inserted on the load side of the coupler in order 
to adjust the tuner for the best attainable match. 

Another assumption used in the above derivation, viz., 
that the coupler is connected in the waveguide system in 
the “normal” manner, is by no means essential; the method 
works equally well if the coupler is reversed. Eq. (1) 
becomes 

I Ed I = I CE, I • 1 1 + I (4) 

but all the subsequent equations stay the same. 

While probably the simplest and quickest method of 
measuring directivity, the above method is not the most 
accurate. For very small W (and D), one can derive the 
simple error formula. 

8 7 dW 

d [Directivity (dB)] - 

where W is in dB. Thus, if W = 0.3 dB, corresponding to 
Directivity — 35 dB, an error of 0.05 dB in W causes a 
1.5 dB error in directivity. In order to improve upon this, 
it was necessary to substitute a sliding load for the sliding 
short in the setup shown in Fig. 24. In this measurement, 
the coupler must first be connected in the usual manner, 
as if the sliding load were the antenna. The load is then 
moved toward or away from the coupler until at least 
one pair of maxima-minima are observed. The dB differ- 
ence between them, W, and the relative power level 
(square-law detector reading) of the minimum (or maxi- 
mum) are recorded. The coupler is now reversed and the 
new power level, which of course is higher than the first 
by the order of magnitude of the directivity, is recorded; 
or, preferably, an accurately calibrated waveguide atten- 
uator is used to reduce the incident power level to yield 
the same relative power indication as that of the minimum 
(or maximum) previously observed. Suppose 

Max Pli ~ the maximum relative power indication ob- 
served as the sliding load is moved with the 
coupler normally connected, 

Fd 2 = the relative power indication with the coupler 
reversed, 

.„d^=10,og(-^). 
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Then 


+ 201«g(i±f^) ,fo,|D| > l.l 

Directivity - 20 log D = / 

(lOI«g(-^) + 201og(l^:|:ril) .(orlD I < IPI 


[See Sect. 5 for the derivation of Eq. (5).] The above 
described method avoids the ‘perfect” termination re- 
quirement of some methods of directivity measurement, 
but it does, unlike the first (sliding short) method, require 
either reversing the coupler, or interchanging the source 
and detector, in order to measure the coupling coeflScient 
(C). The second (sliding load) method has the additional 
disadvantage that one must know in advance whether 
|D|>|p|or|D|<|p|. A more refined method of 
directivity measurement (Ref. 3) uses a sliding load with 
adjustable p so that | p | is made equal to | D |. (The Min 
Pli is a null, W = oo, and the second term in Eq. (5) is 
simply 6 dB.) This device, however, is not commercially 
available. A similar procedure is described in Ref. 4, the 
only difference being that the p of a conventional low- 
VSWR sliding load is adjusted by home-made methods 
(paper clips or adhesive foil) to maximize W. 

4. Calibration Results 

By the sliding short method the measured directivity 
at 2113 Me is 32.5 dB, and by the sliding load method it is 
32.7 dB. Each of these is the average of three readings. 
The frequency variation of directivity was measured by 
the sliding short method and is shown in Fig. 23. The band- 
width for 30 dB or greater directivity is 1.5% which is of 
of the same order of magnitude as that reported in the 
literature (Ref. 1). The smaller directivity maximum at 
2470 Me was not expected, and it is not known at present 
if there is any way in which directivity can be maximized 
at two or more selected frequencies. 

The high-power performance of the coupler will be 
evaluated in future work. This will consist of (1) simul- 
taneous low-power measurements of the waveguide plus 
termination (water load or feed horn) reflection coefiicient 
by both the coupler and a slotted line, and (2) remeasure- 
ment of the coupler directivity at high temperature (60 to 
100° C) to simulate high-power operation. 


5 . Derivation of Eq. (5) 

The voltage detected by the coupler when terminated 
by a sliding short of reflection coefficient p is 

= (CEO (D -h 

using the same notation as in Eq. (1). Therefore, as the 
load is moved in and out we observe 

MaxlE^I = |C£,|(|D1 + jp|) 

Min|E^,| = |C£,|-(|D| - |p|). 

Since our detector is square law, we actually measure 
MaxP,, = |CEip(|D| + |p|)^G^ 

MinP,, = |C£,p(|D|-|p|)^G^ 

IfiD|>|p|, 

I DCE, I (GO^ = % [(MaxPz,0^ + (Min P^O'^]- 

(6) 

Let W = 10 log . Then 

® MinPi,! 

(Min Pi,0^ = (Max (7) 

With the coupler reversed, the power input to the detector 
is 

Pl2 = I CE, I^G^. (8) 

Substituting from (7) into (6) and dividing (6) by (8) gives 

Du-ectivity = 10 log — = h 20 log ^ 

Pl2 ^ 

(9) 

The result for | D | < j p j follows almost immediately by 
changing (6) to 

I DCE, I (G,)^ = % [(MaxP^i)^ - (MinP^O^j. 

(10) 
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V. Deep Space Station Engineering and Operations 


N67 14454 

A. flight Project Support 

J. Orbison and W. J, Kinder 

1 . Lunar Orbiter I Mission 

The Lunar Orbiter 1 spacecraft was launched from 
Cape Kennedy on Wednesday, 10 August 1966. After the 
launch and upon confirmation of lunar injection, Echo 
DSS secured the station for preparation of the first view 
period, 11 August. Commands for the first orbital engine 
de-boost firing were transmitted from the Echo DSS dur- 
ing the early portion of the fourth view period, Sunday, 
14 August. 

Successful reception and video photographic process- 
ing of the spacecraft picture readouts were accomplished 
for the original high orbit and the final low orbit picture 
taking sequences, and continued through 31 August. In 
addition to the successful mission function. Time Synchro- 
nization Correlations were performed between stations 
Echo, Madrid, and Woomera using the Mark I Ranging 
Subsystem. The latter was augmented by a team from 
the National Bureau of Standards at Boulder, Colorado, 


who brought a Cesium Atomic Time Standard for a time 
comparison with the Echo DSS Rubidium Atomic Fre- 
quency Standard. The NBS time correlations were per- 
formed during the Earth-lunar cruise mode; preliminary 
data indicates an accuracy with WWV of 156.4 /as. 

2. Pioneer VII Mission, j . orbison 

The Pioneer VII spacecraft was launched on Wednes- 
day, 17 August 1966. The Pioneer DSS, acting as backup 
to the Echo DSS for Pioneer VI, was designated the 
Prime Station for the launch and early tracking of Pioneer 
VIL The overlap of the Lunar Orbiter and Pioneer VII 
view periods were such that the Echo DSS could only 
have tracked a portion of either spacecraft’s passes for 
any given day. The nature of the early passes of both 
spacecraft made it imperative that full view periods be 
observed. 

Prior to the launch of Pioneer VII, NASA/AMES per- 
sonnel arrived at Goldstone to direct a series of spacecraft 
Type II Orientation Tests. From the Pioneer Ground 
Operational Equipment (SFS 37-35, Vol. Ill, pp. 15-16) 
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located in the Echo DSS control building, commands are 
first entered and verified through the TCP computer, and 
stored in the command encoder awaiting transmit time. 
Actuation of the transmit function, microwaved the com- 
mand from the Echo DSS to the Pioneer DSS transmitter 
and from there to the spacecraft. The spacecraft's telem- 
etry is received at the Pioneer DSS and microwaved to 
the Echo DSS where it is processed through the Pioneer 
GOE demodulator/synchronizer and into the TCP com- 
puter for teletype routing to the Space Flight Operations 
Facility in Pasadena. 

During the Pioneer DSS second pass of Pioneer VII, 
Thursday, 18 August, the Type II Orientation aligned 
the spacecraft high-gain antenna with the Earth. A daily 
tracking schedule was thereafter maintained through the 
end of August. 

Two TCP computers, designated Alpha and Beta, are 
the basic units for the Telemetry and Command Processing 
(TCP-II) subsystem (SPS 37-38, Vol. Ill, p. 102). Included 
are the transfer racks allowing rapid interchange be- 
tween Lunar Orbiter and Pioneer mission functions, the 
NASCOM High-Speed Data Line equipment and the 
associated teletype printers. Normal operation provides 
for Lunar Orbiter data processing through the Alpha 
computer and Pioneer data processing through the Beta 
computer. In the event of a computer failure during over- 
lap tracking, the spacecraft priority at that time would 
determine use of the operational computer. Transfer 
between the computers is accomplished nominally within 
5 min, resulting in minimal loss of data, 

3. Pioneer VI Mission, j. Orbuon 

The Pioneer VI spacecraft is in the 9th month of its 
flight. The Mars DSS assumed the prime tracking of the 
Pioneer VI beginning with Pass 196, 29 June 1966, with 
telemetry and command processing being accomplished 
at the Echo DSS via microwave in a manner similar to that 
described for the Echo/Pioneer DSS arrangement. 

Until Pass 234, 6 August, the Echo DSS also provided 
a transmitter and command function for the Mars DSS. 
Because of prime power difficulties with the diesel gen- 
erators, the Mars DSS delayed testing the 20-kw trans- 
mitter installation on the 210-ft antenna until additional 
prime power was provided. Beginning with Pass 235, 
7 August, the Mars DSS assumed the command trans- 
mission to the Pioneer VI, freeing the Echo DSS for the 
approaching Lunar Orbiter Operation. Telemetry and 
command processing and prime magnetic tape recording 
continues to be performed at the Echo DSS. 


4. Mariner IV Mission, j . Orbison 

Twenty-two months after its launch on 28 November 
1964, the Mariner IV spacecraft is being tracked for sev- 
eral hours on a one-day-a-month schedule. During the 
last scheduled tracking period on 5 August 1966, two-way 
lockup was established using the 85-ft antenna and the 
lOO-kw transmitter at the Venus DSS and the 210-ft 
antenna of the Mars DSS for receiving. Commands were 
transmitted to the spacecraft low'gain omni-antenna, and 
resultant spacecraft telemetry was recorded. 

5. Surveyor B Mission, J. orbison 

The Surveyor 1 Mission was officially terminated in 
July, and the Pioneer DSS began station preparation for 
the Surveyor B Mission scheduled for a September launch 
from Cape Kennedy, Florida. In general, the Surveyor 
equipment and the S-band system have remained un- 
changed. Except for equipment relocation in the S-band 
control room to permit installation of the Multiple Mission 
Support Area (MMSA),only minor assembly modifications 
have been performed. Inter-area compatibility tests were 
performed, and the equipment has been maintained in 
operational readiness for the launch. 

6. MSFN Apollo Implementation at DSS 1 1 , 

42 and 61, W.J. Kinder 

Should a failure or serious degradation of performance 
occur at the Apollo Prime station, mission requirements 
commit the co-located DSN Apollo Backup station to be 
fully capable of substituting for the Prime station within 
a short interval of time. The co-located DSN stations 
(Ref. 1) normally will participate in the Apollo lunar 
missions and their simulations on a full-time basis from 
launch minus 2 weeks to the end of the mission and at 
other times on a scheduled basis. 

The 85-ft antenna Backup station will have a dual 
capability, thereby permitting communication with two 
Apollo spacecraft simultaneously (Ref. 2). Different fre- 
quencies will be assigned to the dual functions; however, 
only a single DSN 85-ft antenna will be used. Dual 
tracking, including data formatting, is accomplished at 
the DSN site and will employ separate NASCOM cir- 
cuits to the computers at MSC's Mission Control Center 
(MCC-H) providing full redundancy to the Prime station. 
Two-way voice is shared between the Prime and Backup 
sites in that while the transmitter and receiver are 
located at the Backup site, the “downlink” voice de- 
modulators are located at the Prime station, with the 
latter s NASCOM circuits used to relay voice to and from 
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NOTES; 

I. BOTH KEY OPERATED MISSION SWITCHES MUST BE ENGAGED, MODE 
MAINTAINED AFTER SELECTION 


LEGEND 

(A) JPL SUPPLIED 

(B) JPL INSTALLED 

(C) 6SFC SUPPLIED 

(D) GSFC INSTALLED 


2, EXISTING FACILITIES, POWER PLANT AND DISTRIBUTION TO BE ENLARGED 

TO COMPENSATE FOR ADDITIONAL MSFN SITE REQUIREMENTS 

3, J- BOXES REQUIRING SWITCHING ARE LABELED- ACTIVATED BY 

DSN/ MSFN MISSION SWITCH 

4, ACQUISITION AID SUBSYSTEM (DSS 42) TEST CONTROLS IN DSN WINGS ONLY 


Fig. 1. MSFN/DSS Apollo Backup station system interface diagram 
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Fig. 2. MSFN/DSS ApoHo Backup station conversion diagram; DSN equipments 
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MCC-H. Two-way data (commands, telemetry, TV, etc.) 
is handled similar to two-way voice. However, for both of 
these functions, the DSN Apollo stations will be imple- 
mented with subcarrier modulators. Subcarrier demodu- 
lation and decommutation will be performed at the Prime 
sites using composite spectra from the DSN Apollo sta- 
tion s carrier phase demodulators and/or the 50 MHz IF 
(FM) output. Selected spacecraft telemetry (e.g., trans- 
ponder phase error and AGC) is routed back to the DSN 
site after processing by the Prime station, A microwave 
carrier two-way communication system links the co- 
located Apollo Prime and Backup stations. 

Selected DSN stations to be modified for Apollo support 
are as follows: 

Goldstone-Pioneer MSFN/DSS 11 

Canberra-Tidbinbilla MSFN/DSS 42 

Madrid-Robledo MSFN/DSS 61 

Table 1 summarizes the major Apollo functional require- 
ments. 

Table 1. Functions of the MSFN Apollo Backup 
station system 


(1) Acquire and track spacecraft signals through an integral 
acquisition subsystem. 

(2) Transmit voice, data, and ranging signals to spacecraft. 

(3) Transfer track from the acquisition subsystem to the main 
receiving subsystem. 

(4) Determine range and range rate parameters of spacecraft 
and transfer this information to recorders and/or the NASA 
Communications Network. 

{ 5 ) Receive video, voice and telemetry information from the 
spacecraft, and transfer this information to recorders and/or 
the NASA Communications Network. 

(6) Accept pointing information from specified sources (drive 

tapes, manual) and utilize this information to drive the 
antenna. 

(7) Maneuver or position the antenna according to any manual 
drive inputs within the specified maximum rate and ac> 
celeration. 

(8) Operate in programmed search modes. 

(9) Follow a second antenna in slaved mode. 

(10) Provide slaving signals for driving other antennas. 

(11) Provide visual readout of antenna pointing angle, record 
this information, and/or transfer this information to the 
NASA Communications Network. 


A Master Mission control switch automatically changes 
the DSS equipment configuration to either MSFN or 
DSN mission support. Except for the site intercom/ 
communication manual patches and transmitter/maser re- 
tuning, the DSN/MSFN station will be in one configura- 
tion or the other, by the activation of the Master Mission 
control switch featuring the following singularities: 

(1) Single activation; requiring the simultaneous opera- 
tion of two identical keys — one in each control room. 

(2) Dual location; DSN and MSFN control room. 

(3) Mission mode maintained after selection. 

After key operation, four indicator lamps will signify 
switching center closures, i.e., the microwave subsystem; 
the angle encoder interface; and servo interface; and the 
Optical Tracking Assembly. A remote wind velocity and 
direction display will be furnished to the MSFN servo 
console from the common site moimted sensors. 

Initial considerations in the JPL approach to the mission 
requirements were to perturb the present DSN S-band 
installation as little as possible and to assess or utilize 
existing GSFC Prime site engineering. In Fig. 1, the 
MSFN/DSS Apollo Backup station system interface dia- 
gram delineates the responsibilities in both equipment 
and installation shared between JPL and GSFC. Only 
the JPL equipment engineering conversion will be dis- 
cussed and is shown in Fie. 2. 

The 85-ft antenna microwave subsystem equipment 
includes two traveling wave maser (TWM) low-noise am- 
plifiers and associated closed-cycle refrigeration; S-band 
Cassegrain monopulse cone; test signal control assembly; 
monitor receiver assembly; and required microwave hard- 
ware, cables, and waveguide. 

A microwave switch junction box, operated by the 
master mission switch and installed in the antenna elec- 
tronics room, permits either the MSFN or DSN use of 
the common antenna-mounted microwave equipments. 
Similarly, located in the DSN control plenum, a status 
signal switch junction box routes DSN microwave con- 
figuration status signals and receiver IF attenuator com- 
mands to either network. Accidental radiation from the 
MSFN collimation tower equipment will be prevented 
by an interlock controlled by the status signal switch. 
Either the MSFN or DSN wing can be in a total opera- 
tional system configuration simultaneously with the other 
networks system in a partial or independent test configu- 
ration mode. The common test signal control assembly 
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will permit test signals from either of the two test 
sources, i.e., transponders or exciter translators, to be 
injected into the MCFN crystal mixer input circuits of the 
multiple receiver reference and angle channels. 

For simultaneous transmission of two separate carrier 
signals, two 20-kW transmitters will provide MSFN 
S-band frequency power amplification (SFS 37-38, Vol. Ill, 
p. 72). Combining two signals at high level, each 20-kW 
transmitter will dissipate 50% of its power in a 3-dB 
hybrid combiner waterload with the remaining power 
being radiated. The two transmitters will be controlled 
from a new transmitter rack placed in the DSN control 
room adjacent to the microwave subsystem. Status signals, 
indicative of transmitter operation, showing beam voltage 
applied to the power amplifiers are routed to the MSFN 
control room for monitoring purposes. 

GSFC will use the existing GSDS Optical Tracking Aid 
(OTA) consisting of cameras, filters, lenses, /-boxes, cables, 
and control panels with camera controls and aperture 
door operation available to either network monitor via 
the mission switch-operated optical tracking and switch- 
ing junction box located in the DSN control room. This 
arrangement assures only one optical reference system to 
align the mechanical axes. The Optical Acquisition Aid 
(OAA) assembly will remain unchanged. 

Structural modifications to the antenna electronics room 
to accommodate the additional MSFN equipments con- 
sist of a new catwalk and guard rail to facilitate servicing 
of the antenna junction box and installation of new floor 
beams. Additional antenna cable trays and wrapup drums 
are also required. Three spare SCM cones complete with 
microwave components (less maser) will be furnished. 
One cone will be stored at each individual DSN Apollo 
site in its original shipping container under an environ- 
ment of pressurized dry nitrogen. 

The GSDS S-Band Acquisition Aid Subsystem (SAA) 
requires modifications to accommodate both MSFN and 
DSN frequencies. This will be accomplished by broad- 
banding present waveguide components. Modifications 
will be accomplished at Goldstone and Tidbinbilla only, 
since no Acquisition Aid is available at the Madrid 
station. 

An angle encoder switch junction box routes antenna 
angle data digital signals to either the DSN Tracking Data 
Handling Subsystem or GSFC s Antenna Position Pro- 
grammer. Analog angle data signals for the GSFC servo 


system ball tracker are circuited through isolation ampli- 
fiers. This switching center will be located in the present 
TDH rack in the Hydro-mechanical building. 

To permit the antenna drive signals originating in either 
control room to control the common antenna hydraulic 
drive assembly, a servo switch junction box, located in 
the Hydro-mechanical building pump room is required. 
All three antenna sites will be retrofitted with new DC 
tachometers. 

A new SCM collimation tower 50-pair cable will be 
installed adjacent to the existing DSIF collimation tower 
cable for MSFN use. 

Low-speed TTY data telephone pair cables connect the 
GSFC Tracking Data Processor and the DSS communi- 
cation room. The JPL Tactical Intercom and Communica- 
tion Assembly (TIC) will also be interfaced with the 
GSFC 112- A intercom subsystem via a cable to the latter s 
vertical distribution frame and by appropriate pre-mission 
manual telephone patches. This arrangement will allow 
GSFC to use the presently installed DSS telephone instru- 
ments. 

J PL-provided facilities will include floor space for the 
MSFN configuration at the selected sites. Utilities, com- 
prising power, air conditioning and water in amounts 
required by GSFC for proper operation of the MSFN 
equipment, will also be provided. 


AI67 14455 

V 

B. Facility Construction and 
Equipment Installation 

J. Orbison and H, L. Richter 

1. Mars Deep Space Station , J.Orbison 

The installation of the S-band system continues con- 
currently with the Mars DSS operation. Originally, the 
operation was performed using the JPL R&D receiver 
and maser, and the Station Digital Instrumentation Sub- 
system. The transmitting to the tracked spacecraft was 
performed by the Venus DSS and the Echo DSS when 
required. Both of the latter stations also provided the 
necessary recording and processing via microwave from 
the Mars DSS. 
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Essentially, the Mars DSS has received the total S-band 
system and is currently completing the installation. The 
major exception is the use of the JPL R&D maser with 
the DSIF receiver for current tracking of the Pioneer VI 
spacecraft. The 20-kW transmitter was placed in opera- 
tion on 7 August. The high-voltage power supply for the 
transmitter currently is located in the large service area 
on the first floor of the pedestal (Fig. 3). Eventually it will 
be relocated to its permanent position in the main equip- 
ment room on the alidade. 

The master equatorial assembly was positioned in the 
room atop the central instrument tower early in June. 
Currently it is in the installation and checkout phase of 
leveling and alignment. 


The diesel generators, which provide the present power 
to the Mars DSS, experienced excessive torque effects in 
the rotors and required disassembly for rework. An 
1126-kW diesel generator unit (Fig. 4) was temporarily 
installed near the generator building to provide station 
power while the three generators were undergoing repairs. 
The No. 2 generator is in operational standby for emer- 
gencies and will itself be reworked when the other two 
become operational. 

Construction of a commercial power line from the re- 
cently completed Southern California Edison substation, 
5 miles east of the Venus DSS, to the Mars DSS was 
begun in late July. Scheduled for a fall completion, it 
will provide adequate power requirements for planned 


Fig. 3. Mars DSS transmitter high-voltage power supply 
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Fig. 4. Temporary 1 126-kW diesel generator at Mars DSS 


station operations, and will also serve as backup to the 
Station’s diesel power which is used during mission- 
critical periods. The pole line will terminate approxi- 
mately 1 mile south of the 210-ft antenna, and the cabling 
will proceed underground to the substation site near the 
diesel generator building. 

2 . Pioneer Deep Space Station, j Orbison 

The Pioneer DSS requirements to maintain support to 
the normal DSN missions, and also help in the implemen- 
tation of the MSFF, required extensive station changes. 
Initially, these dual requirements began with an enlarge- 
ment of the Control Building (SPS 37-27, Vol. Ill, p. 5; 
37-28, Vol. Ill, pp. 13-15; 37-29, Vol. Ill, pp. 10-11). Two 
opposed identical wings were constructed with equal 
tunnel access to the hydro-mechanical building, and the 
85-ft antenna. 

Installation of the MSFF S-band system began in 
July 1966 and is in progress in the west wing of the con- 
trol building (Fig. 5). Equipment installation to date 
includes the receiver and transmitter exciter, antenna 
servo electronics, antenna position programmer, tracking 
data processor, and portions of the telemetry equipment. 


Subsystem power-on tests, preliminary interfacing and 
basic operational tests are in progress. 

The additional weight of the second transmitter and 
the maser assemblies requirements of the MSFF on the 
antenna required the addition of 4900 lb of counter- 
weights to the west side of the hour angle gear. Equip- 
ment relocations, new cable trays, and other changes to 
the declination room and the antenna were also made. 
A cable wrap-up, patterned after the one developed on 
the Echo DSS 85-ft antenna (SPS 37-35, Vol. Ill, p. 17) 
has been installed and tested. 

The original 10-kW Klystron power amplifier was 
changed to a 20-kW unit, necessitating modifications to 
the high-voltage power supply located in the hydro- 
mechanical building to provide the additional power. 
The MSFF 20- kW Klystron power amplifier and associ- 
ated equipment were also installed in the declination 
room, while the high voltage power supply was installed 
adjacent to the original unit. Near the base of the north- 
east leg of the antenna are the two heat exchangers for 
the transmitter power amplifiers, with a newly fabricated 
header permitting either unit to be used with either 
transmitter. 
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Fig. 5. Manned space flight facility control room installation 


I 



Fig. 6. Dual 20-kW transmitter simplified waveguide switching 


A simplified block diagram of the dual transmitter Cassegrain Tracking Feed or the Acquisition Aid an- 
installation is illustrated in Fig. 6. The switching arrange- tenna. Through the use of the power combiner, both 

ment permits either transmitter to be coupled to the transmitters can feed the antenna simultaneously. 
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The parametric amplifier was removed to provide space 
for mounting the MSFF maser in the declination shaft 
room. By using a locally fabricated DSN/MSFF RF 
switching assembly, complete flexibility between the two 
systems is possible. Either system can control the antenna 
drive, use either or both masers with the respective 
receivers, transmit individually or simultaneously, and 
perform the system transfer quickly, when needed. 

Equipment relocations were performed in the DSN 
control room to provide space for the Multiple Mission 
Support Area, and the dual capability needed for the 
DSN/MSFF function. Relocated were the Analog Instru- 
mentation subsystem and the Antenna Microwave ground 
control racks. An additional Analog Instrumentation sub- 
system, recorders, and the TCP computers for the Telem- 
etry and Command Processing subsystem will be installed 
at a later date. 

3. Echo Deep Space Station, J.orbison 

Lunar Orbit er and Pioneer VII equipment freezes 
delayed the installation of the second Analog Instrumen- 
tation Subsystem delivered just prior to the final testing 
and launch of the Lunar Orhiter. The equipment racks 
were placed in position in the control room and readied 
for interrack and interface cabling. Two Ampex FR-1400 
magnetic tape recorders for the second system arrived in 


late June and were installed and placed in operation for 
Pioneer spacecraft recording, a function they are cur- 
rently performing for Pioneers VI and VII. 

4. Venus Deep Space Station, j. orbison 

Construction of the Laboratory and Office Building, 
G-60, is in progress (Fig. 7). The outside walls, and roof 
are in the final stage of construction, and the interior 
features installations are in process. 

5. DC Tachometer Installation on DSIF 85-ff 
Polar Antennas, h. l. Richter 

In order to improve the interface of rate feedback 
signals with MSFN servo electronics and to upgrade the 
rate loops of the standard DSIF 85-ft diameter polar 
tracking antennas, it was decided to replace the presently 
used 400-Hz AC tachometers with DC tachometers. The 
DC tachometer offers a number of advantages over the 
AC machine when used in a basically DC servo system 
as the DSIF employs. These advantages manifest them- 
selves in the fact that the need for field excitation and 
demodulation is eliminated thereby making the system 
simpler. By elimination of the above, the noise level of the 
feedback signal is decreased to that of pickup and any 
phase shifting because of noise is eliminated. The DC 
tachometer also provides better feedback sensitivity at 



Fig. 7. Laboratory and office building, G-60, at Venus DSS 
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low rates because the signal is not degraded by the 
demodulator. Finally the DC tachometer, in the size used, 
is much less expensive, costing approximately 25% as 
much as the comparable AC machine. 

The DC tachometer chosen was the Servo-Tek model 
SB-740B-1. This machine is a standard DC permanent 
magnet tachometer generator with built-in temperature 
compensation and a nominal gain of 20.8 V/1000 rpm. 
This tachometer was chosen on the recommendation of 
Collins Radio Company, who supplied the MSFN system 
and who have gained experience on a number of an- 
tennas. The tachometer was laboratory tested for gain 
and life characteristics. The machine tested was found to 
be linear at approximately 21 V/1000 rpm from 0.5 to 
3000 rpm with a signal-to-noise ratio of 0.5 at 1 rpm and 
a ratio of 1 at 25 rpm, increasing steadily upward from 
there. The tachometer was life-tested for 1000 hr at a 
nominal rate of 2000 rpm and showed no visible loss of 
gain or increase of noise. Catalog life on the tachometer 
brushes is given as 100,000 hr, and bearings are assumed 
to exceed this. 

The testing and trial installation of the tachometers 
was performed on the DSS-12 antenna declination drive. 
The declination axis was chosen because it could be 
modified without jeopardizing any tracking commitments 
of that station since the declination axis does not move to 
any appreciable extent when tracking deep space targets 
and because the theoretical models of the hour angle and 
declination rate loops are nearly identical. The only dif- 
ference between the two is that the hydraulic poles are 
slightly difiFerent in location, this being due to different 
line lengths between the servo valve and the motors. 
The difference can be considered negligible, however. 

The modification to the servo group electronics involves 
removal of all demodulators, which demodulated the 
400-Hz rate signal, and installation of new wiring for the 
additional tachometers added to motors not so equipped 
in the old configuration and designing the proper com- 
pensation for the loops to obtain maximum gain and 
bandwidth. 

The low-speed rate loop was compensated with a 9:1 
lead-lag compensator with lead comer at 40 rad/s and a 
10:1 lag-lead compensator with lag comer at 0.63 rad/s 
for additional phase margin and DC gain, respectively. 
The final Bode diagrams for the low-speed rate loop are 
shown in Fig. 8. The final open-loop DC gain constant is 
800. For the high-speed rate loop, it was necessary to 



Fig. 8. Declination low-speed closed-rate loop 


close the pressure loop around the servo valve in order 
to obtain additional phase margin and hydraulic damp- 
ing. The high-speed pressure loop was compensated 
with a 12:1 lead-lag compensator, the lead comer at 
28 rad/s, resulting in the closed pressure loop transfer 
function shown in Fig. 9. After closing the pressure loop 
as indicated, the high-speed rate loop was closed using 



Fig. 9. Decimation high-speed closed-pressure loop 
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a 10:1 lag-lead compensator with the lag comer at 0.63 
rad/s. This resulted in the Bode diagram of Fig. 10. The 
final open rate loop DC gain is 230. 

Response curves were obtained from magnitude and 
phase data measured on a Boonshaft and Fuchs transfer 
function analyzer. 

The physical change to the servo amplifier chassis was 
planned so that no mechanical modifications to the chas- 
sis or front panel were required. Sockets formerly used 
for demodulators were rewired to become sockets for the 
compensation and summing networks for the new and 
existing feedback signals. Planning the modification in 
this manner allowed the change to be made without any 
visible change to the equipment. Detailed installation 
procedures are set forth in JPL Procedure DZA-1026-IP. 

The DC tachometer modification is now being readied 
for installation on all DSIF 85-ft diameter antennas, 

C. DSIF Station Control and Data 
Equipment 

R. Flanders, G, Jenkins, A. Burke, H. Baugh, 
and W. Frey 


J, Introduction 

The DSIF Station Control and Data Equipment 
(SCDE) is to be implemented to provide increased sup- 
port for the flight projects as a part of the DSIF station 
equipment. In the performance of this objective, the 
equipment is to provide the capability for processing 
spacecraft telemetry and command data in direct sup- 
port of the mission, and for station monitoring, time and 
frequency tagging, and antenna pointing as a part of the 
mission independent functions. The specific subsystems 
involved with these tasks are the Antenna Pointing Sub- 
system (APS), Digital Instrumentation Subsystem (DIS), 
Station Monitor and Control Console (SMC), Telemetry 
and Command Processor (TCP), and the Frequency and 
Timing Subsystem (FTS). The function of each of these 
subsystems, together with the current status, is described 
in the following paragraphs. 

2 , Antenna Pointing Subsystem 

The Antenna Pointing Subsystem (APS) provides the 
DSIF stations with the capability to point antennas under 



Fig. 10. Declination high-speed closed-rate loop 
transfer function 


computer control in addition to the existing automatic 
tracking and manual control modes of operation. The 
primary function of the subsystem is to point an antenna 
toward a spacecraft or celestial body ephemeris that is 
either derived internally based on parameter inputs or 
supplied by external data contained on precomputed 
drive tapes. The subsystem consists of an SDS 910 com- 
puter with associated peripheral equipment, interface 
equipment and a subsystem control panel. The APS oper- 
ates in conjunction with the servo assembly of the An- 
tenna Mechanical Subsystem and the antenna angle 
readouts of the Tracking Data Handling Subsystem. The 
APS will add the following capabilities to the DSIF: 
(1) Rapid acquisition time. (2) Tracking when there is 
no RF signal or the received signal strength is too low 
for automatic tracking operations. (3) More flexibility and 
precision in control of antenna movement. 

The Antenna Pointing Subsystem is being implemented 
in two steps. The first step, the Interim Antenna Pointing 
Subsystem, meets the flight project antenna pointing re- 
quirements on an interim capability basis. The second 
implementation step will provide a full capability An- 
tenna Pointing Subsystem (APS 1). The interim APS 
utilizes existing equipment to perform its function, while 
the APS 1 requires the addition of new computer inter- 
face equipment to operate. Both the interim and full 
capability APS are described in detail in SPS 37 - 38 , 
Vol. III. 
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The Interim APS has been implemented in DSIF Sta- 
tions 12, 14, 41, 42, and 61. Station 11 will be implemented 
in September 1966, and Stations 62 and 72 will be in- 
stalled in the last quarter of this year. The Interim APS 
is presently providing antenna pointing support for the 
Lunar Orhiter Flight Project at Stations 12, 41, and 61. 

The first implementation of the APS 1 equipment 
will be during October 1966 at Station 14. Stations 11, 
12, 41, 42, 51, 61, 62, and 72 will be implemented in 
1 -month intervals beginning in May 1967. At the present 
time, hardware evaluation and program checkout are 
being conducted on the prototype APS 1 which is located 
at the Laboratory. 

3. Digital Instrumentation Subsystem, Phase II 

The Digital Instrumentation Subsystem Phase II (DIS 
II), together with the Station Monitor and Control Con- 
sole Phase II (SMC II) provides the facility for system 
monitoring at each station in the DSIF. The DIS/SMC II 
performs the following functions: (1) Station performance 
and status monitoring. (2) Alarm monitoring of critical 
station parameters. (3) Preparation of a permanent record 
of station performance and failure data for real time and 
post-mission analysis. The DIS II provides for the acqui- 
sition, processing, recording, and communication of sta- 
tion data. 

The DIS II is implemented with the major elements of 
the DIS I configuration (SPS 37-16, -20, -21, and -34, 
Vol. Ill), with the exception of the Communications 
BuflFer which has been assigned to the Telemetry and 
Command Processor (TCP), and the SDS 910 computer 
and peripheral equipment which is being utilized as part 
of the Antenna Pointing Subsystem (APS). The DIS II 
utilizes the SDS 920 computer, peripheral equipment, 
magnetic tape transports, and data gathering equipment 
from the DIS I. To this is added the following elements: 
(1) Y-BuflFer interlace. (2) Memory expansion of 8,192 
words of magnetic core memory to increase memory 
capacity to 16,384 words. (3) Monitor Page Printer and 
Coupler. (4) Line Printer and Coupler. (5) Cartridge 
Loading Magnetic Tape Unit. (6) Communications Buffer 
with full duplex teletype and high-speed data line com- 
munications capability. 

A prototype of the DIS II is being procured and will 
be utilized for: (1) evaluation of the configuration and 
performance prior to procurement of production quanti- 
ties, (2) verification of the subsystem documentation, 
(3) checkout of the operational programs and diagnostics. 


and (4) establishing interface compatibility with the mat- 
ing subsystems in the stations. Completion of the proto- 
type evaluation is scheduled for the fourth quarter of 
1966. Subsequent to the completion of the evaluation, a 
production quantity of the DIS II will be implemented. 
This will require the retrofit of the existing DIS I at each 
of the stations involved. The initial installations will be 
performed at DSS 11, 12, 61, 62, and 71 and are scheduled 
for operation during the first quarter of 1967. The re- 
maining stations are scheduled to be updated to the 
DIS II configuration during the latter part of 1967. Sup- 
porting documentation is to be prepared, for completion 
simultaneously with the hardware availability. 

The program for the DSIF System Monitor function is 
being prepared to replace the Interim Monitor Program 
currently used at the stations. The Interim Monitor Pro- 
gram is presently operating with the DIS I and performs 
limited monitoring functions. The permanent System 
Monitor Program will be available for use at selected 
stations when the DIS II becomes operational, and will 
provide a full monitoring capability, including perform- 
ance and alarm monitoring, data gathering and proces- 
sing, and station report preparation and communication. 

4. Station Monitor and Control Console 

The primary objectives of the Station Monitor and 
Control Console Phase I (SMC I) are: (1) to provide the 
facility for the display and control of station performance 
at a central location, for use by the station manager at 
each site, and (2) to display information to indicate non- 
standard system or subsystem operation. SMC Phase II 
will add the capability of: (1) displaying selected station 
parameters, (2) an English text printout of station per- 
formance and status, (3) system alarm monitoring under 
control of the DIS to indicate station failures, and (4) a 
teletype page printer. SMC I is essentially complete with 
exception of the countdown clock. Procurement of the 
clock is in process and should be complete by May 1967. 

Procurement of a prototype Program Alarm and Con- 
trol Panel has been initiated and is scheduled for com- 
pletion by September 1966. An order for six additional 
panels will follow upon completion of the prototype 
evaluation. The prototype line printer has been procured 
and will be shipped to DSS 11 for testing and evaluation. 
Thereafter, five stations will be equipped with line 
printers as part of SMC II. 

A graphical recorder will be mounted in the station 
monitor console to provide the facility for displaying 
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actual data versus predicted performance data on an 
X-Y plot. Procurement of a prototype graphical recorder 
has been initiated, for evaluation of the recorder prior 
to procurement of an additional quantity. The X-Y re- 
corder will be evaluated at JPL and at DSS 11 to verify 
proper performance, establish interface compatibility with 
companion equipment, and assure the validity of the 
supporting documentation. Prototype evaluation is to be 
completed during the fourth quarter of 1966, with oper- 
ation at five stations scheduled for the first quarter of 
1967. 

5. Telemetry and Command Processors, Phase II 

The Telemetry and Command Data Processors Phase II 
(TCP II) are the major assemblies in the Telemetry and 
Command Data Handling Subsystem Phase II (TCD II). 
A description of the TCP II has been presented in SPS 
37-38, Vol. Ill, Sect. VI-C. The major changes in TCP II 
during this period are as follows: 

(1) The TCP II-A configuration has been eliminated 
and the stations formerly configured; namely DSS 
11, 42, and 62, will now be implemented with the 
standard TCP II configuration. Table 2 shows the 
modified TCP II configuration. 


Table 2. DSIF TCP II configuration 


DSIF 

Idtntificafien 

TCP 

configuration 

Operational 

date 

n 

11 

December 1966 

12 

II 

February 1966 

41 

n 

March 1966 

42 

II 

November 1966 

51 

ll-B 

October 1965 

61 

II 

March 1966 

62 

11 

December 1 966 

71 

n 

March 1 966 

72 

ti 

April 1 966 


(2) The TCP II-A scheduled for installation into DSIF 
11 in December has temporarily been installed in 
the Spacecraft Assembly Facility (SAF) at JPL for 
support of the Mariner 67 Project. This unit will 
be returned to the DSIF in January 1967. 

(3) TCP II has been participating in Flight Project 
Support at four stations. DSS 12, 41, 61, and 71 
have operational TCP II assemblies, and these 
stations have been utilizing the TCP II in provid- 


ing support to the Lunar Orbiter and Pioneer flight 
projects. At the present time, not enough data has 
been processed by the TCP II for a quantitative 
analysis. 

6. Frequency and Timing Subsystem, Phase II 

a. Background. The present Frequency and Timing 
Subsystems (FTS I) has some serious deficiencies in the 
areas of reliability, resolution and accuracy, and drive 
capability. Increasingly stringent mission requirements 
have indicated the need for a new subsystem, FTS II, 
which can correct these difficulties. 

b. Reliability, In order for the FTS to provide correct 
time information, all elements must operate essentially 
error-free for long periods of time. The digital elements 
in the present subsystem can not do this satisfactorily, 
since the modules are of an older design and the logical 
structure makes no provision for error detection or cor- 
rection. 

FTS II will be implemented with redundant assemblies 
using Hi-Rel Modules so that the probability of a module 
failure is greatly reduced and the subsystem can continue 
to operate correctly even if failures do occur. An example 
of this type of structure is the divider chain which pro- 
duces 1 pps from a 1 MHz input. This assembly will 
consist of three identical divide-by 10^ counters with 
triplicate majority logic voters operating on the 1 pps 
signals. If two of the signals agree and one disagrees, 
the good outputs will be used, the dissenting signal 
suppressed, and an alarm indication given. The offending 
element may be “locked out,” requiring operator inter- 
vention, or it may be automatically synchronized to the 
two and re-connected. 

c. Resolution and accuracy. The present FTS digital 
clock has provision for advancing or retarding to set 
the clock to the proper time in relation to WWV, but the 
changes can only be made in increments of 100 /xs, so 
the time can not be set to closer than ±50 fxs. The tim- 
ing pulses out of the unit have considerable phase jitter 
in view of present requirements, so there is a need for 
improvement in this area. 

In FTS-II, it is planned that the clock-setting resolu- 
tion be 1 /xs, as there are an increasing number of track- 
ing operations where this order of resolution is important. 

The phase jitter of the timing pulses from FTS-I is a 
consequence of the use of slow-speed digital modules 


128 




JPL SPACE PROGRAMS SUMMARY NO. 37-41, VOL. Ill 


and unclocked logic for the generation of the key pulses. 
FTS-II modules will have much shorter delay times, and 
all key pulses will be directly clocked by the basic fre- 
quency source so that the phase jitter should be less than 
a few tens of nanoseconds. 

The absolute accuracy of the timing information is a 
function of the basic input frequency. FTS-I uses ru- 
bidium frequency standards, which can nominally pro- 
vide the necessary accuracy, but they have presented 
some reliability problems. FTS-II will use similar stand- 
ards, initially, although some work is being done by the 
vendors to improve the reliability of the units. 

d. Capability. As the tracking stations grow in com- 
plexity, there is a continual demand for additional 
services and additional driving capability on existing sig- 
nals. This generally leads to overloaded drivers and little 
or no isolation between subsystems. Overloading in some 
cases has caused malfunctions in the station clock. 

The digital dividers in FTS-II will be completely iso- 
lated from all output circuits by short-circuit-proof ampli- 
fiers, so that no external activity can upset the basic 
timing chains. Similarly, the digital clocks will be isolated 
from their customers by amplifiers which also provide 
the level-shifts required by the difiEerent systems of logical 
elements used in other subsystems. Key timing pulses and 
frequencies will have individual distribution amplifiers 
for each customer, to aid in preventing interaction be- 
tween subsystems. 

An additional capability to be provided by FTS-II is 
that of flexibility. In the checkout and simulation area 
there are requirements for switching the sources of tim- 
ing signals so that one project may be checking its 
equipment with tape playback signals while another 
project may be tracking with real-time signals. All com- 
binations of timing sources and customers should be 
available by suitable switching, without requiring that 
the FTS-II digital clock itself be interfered with in any 
way. 

e. Status. At the present writing, all the functional 
requirements for FTS-II have been established, and the 
basic design philosophy has been defined. Procurements 
have been initiated for the digital modules for the proto- 
type unit, and for the documentation. Investigations have 
been started in the area of long-term stability of fre- 
quency standards, and in the area of techniques for 
synchronizing overseas stations to a central timing sys- 


tem to be established at Goldstone as a portion of FTS-IIA 
and FTS-IIB. 

7, Summary 

A limited quantity of the Station Control and Data 
Equipment has been operational in support of recent 
flight projects. Use of the equipment during these mis- 
sions has provided experience for the station personnel, 
as well as an indication of expected performance when 
implementation of the final configuration is complete. As 
a result, some improvements have been incorporated into 
the development. For the most part, the performance has 
been favorable. 


N67 U457 

D. Venus Station Operations 

M. A. Gregg, E. B. Jackson, and A. L. Price 

1 . Experimental Activities 

During the period 14 June to 30 June, the Venus Sta- 
tion was not in an operational configuration; instead, it 
was conducting maintenance, modifications, upgrading, 
and station cleanup. The monthly cooperative (with the 
Mars Station) experimental track of the Mariner IV space- 
craft took place on 1 July, and full operational status 
with a return to 24-hr/day operation started on 5 July. 

The cooperative, with the Mars Station, tracks of the 
Mariner IV spacecraft (on 1 July and 8 August) returned 
one- and two-way telemetry data (recorded on magnetic 
tape at the Mars Station), spectograms (via microwave 
to the Venus Station), and some two-way doppler data 
to be utilized in orbit analysis. The spacecraft transmitted 
power and frequency in the two-way mode appear to be 
nominal; however, the spacecraft in the one-way mode 
is being received at a frequency of approximately 600 Hz 
(at carrier frequency) different from predictions. The 
exact amount of difference and possible reasons are still 
under study. 

During the period of 12 July to 5 August, a bistatic 
planetary radar experiment with Mercury as the target 
was conducted at the Venus Station. The transmitting 
antenna was the Venus 85-ft Az-El, and the receiving 
antenna was the Mars 210-ft Az-El. The signal was fed 
to Venus at 455 kHz via the intersite microwave link, 
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with both ranging and total spectrum data being col- 
lected. The round-trip range to Mercury was found to 
be initially some 5000 /is less than predicted. However, it 
moved into closer agreement with prediction until, at 
the end of the experiment, the range error was only 
2200 fis less than predicted. This range error must be 
corrected by additional system delay between the Mars 
and Venus Stations via the microwave link. This delay is 
still undergoing measurement but is on the order of 
150 /Lts. Thus, at the end of the experiment, the actual 
round-trip range was less than the predicted round trip 
range by approximately 2350 /xs, or approximately 438 
miles. 

The Mercury bistatic planetary radar experiment was 
complicated by the present errors associated with the 
210-ft antenna pointing system which, at a maximum, 
approximate one-half of the 3-dB beamwidth in both 
elevation and azimuth axes. A day of experimental time 
was utilized to obtain oflFset data by which the nominal 
pointing data could be modified. This data was obtained 


by taking signal strength measurements, integrated over 
a 5-min period, from Mercury. These measurements were 
taken in groups of three arbitrary oflFsets in both axes 
and then a parabola was fitted through the three points 
thus obtained. The resulting oflFsets versus indicated 
azimuth and elevation are illustrated in Fig. 11. 

2. Subsystem Performance 

a. 100-kW transmitting systems. The R&D transmitter 
for this reporting period has been operated only during 
the month of July in conjunction with the Mercury bi- 
static experiments. During this time the transmitter has 
operated a total of 358.5 filament hours and 192.6 beam 
hours. Some internal arcing and outgassing of the tube 
were experienced after the shut-down period. 

For this period, the Mariner 100-kW transmitter has 
operated a total of 61.8 filament hours and 38.4 beam 
hours. However, during testing of the system the Mariner 
Klystron developed a short from the cathode to body, 
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and on 20 June S/N 5 was removed and replaced with 
S/N 4. During these hours of tracking and testing, only 
one other failure occurred, which was the waveguide 
switching at the start of a Mariner track. A total of 2 min 
of tracking time was lost. During this period the system 
has also been used in evaluation of the design for the 
solid-state regulator and exciter for the 1 MW synchron- 
ous motor and generator. 

The testing of the new X-3060 Eimac 100-kW Klystrons 
has been in progress during this period. A leak developed 
in S/N H5-30 beyond the normal pumpdown capability 
after 77.3 filament hours and 16.1 beam hours, and the 
tube was returned to Eimac. 

Support was given to testing the standard DSN diplex- 
ing cone at 100 kW. 

b. Receiving systems. The X-band system 8448 MHz 
was used in a lunar radar experiment. One module 
failure, a radio frequency keyer module, occurred. This 
was repaired at the Laboratory and re-installed in the 
receiver. 

The Venus Station Mariner receiver was not used dur- 
ing this period. Mariner spacecraft reception was accom- 
plished by feeding the 455 kHz signal from the Mars 
Station Mariner receiver, via the microwave link, into 
the AM channel of the Venus Station Mod IV receiver. 

Reception of signals during the Mercury bistatic ex- 
periment was handled similarly. The 456 kHz signal from 
the Mars Station 2388 MHz receiver was fed, via the 
microwave link, into the AM channel and the main loop 
channel of the Mod IV receiver. 

The 30 MHz portion of the Mod IV receiver was used 
in obtaining spectrograms of the Venus Station Mariner 
transmitter transmitted signal. This was accomplished 
by feeding 30 MHz from the RWV receiver into the 
Mod IV receiver. The Mod IV receiver s third conversion 
oscillator, a frequency synthesizer, was used to center 


the spectrum within the passband of the digital processing 
equipment. 

A failure in the programmed local oscillator, resulting 
in large doppler errors in the digital subsystems, occurred. 
It was found that by substituting an external frequency 
synthesizer for the one used in the PLO, the trouble was 
corrected. A new synthesizer was installed in the PLO, 
and the operation is now normal. The old synthesizer was 
returned to the Laboratory for repair. Some tape reader 
problems have occurred. Corrective work is underway. 

3. System Improvements 

a. Transmitting systems. The Owens Lab 500- V power 
supplies which were purchased for the new R&D exciter 
antenna-mounted equipment have both failed in less 
than 48 hr of operating time. 

New control circuitry associated with the operation 
and protection of the pumps for the 400 kW system was 
designed, installed, and documented. 

b. Receiving systems. The installation of all the Venus 
Station receiving systems has been completed during this 
period. 

The X-band digital phase modulator has been modi- 
fied to facilitate its proper adjustment. The modification 
involved the removal of one of the external adjustable 
controls, thus simplifying the setup procedure. 

The 30 MHz leakage problem in the Mod IV receiver 
has been corrected. This was accomplished by rerouting 
the AM channel 30 MHz mixer AGC cable. 

The system checkout of the new central frequency 
synthesizer (CFS), by the Laboratory Project Engineer, 
has been completed. The unit has now been turned over 
to the station for operation and maintenance. The old 
interim CFS has been returned to the Laboratory. 
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